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LIMIT DESIGN REINFORCED CONCRETE BEAMS 


SYNOPSIS 


The paper presents complete analyses the inelastic behavior two- 
span continuous reinforced concrete beam, for both simple plastic and strain- 
hardening, moment-rotation relationships. The analyses predict rotation re- 
quirements that may compared with limiting rotation capacities design. 


INTRODUCTION 


Inelastic concepts are becoming increasingly important structural engi- 
neering. has been evident research workers for some time that inelastic 
theories are essential adequate understanding the actual behavior 
structures. More recently, designers have discovered that plastic and ultimate 
strength design methods happily combine the advantages rationality and 
economy. 

The plastic theory structures is, course, generally applicable any 
structural material with suitable elastic-ductile, load-deformation character- 
istics. Both structural steel and reinforced concrete can shown, with cer- 
tain limitations, possess the requisite characteristics. These limitations— 
essentially the tendency instability steel, and limited ductility rein- 
forced concrete—have stimulated flood research that opening the way 
general use plastic concepts structural design. 


Note.—Discussion open until December 1960. extend the closing date one 


month, written request must filed with the Executive Secretary, ASCE. This paper 
part the copyrighted Journal the Structural Division, Proceedings, the Ameri- 
can Society Civil Engineers, Vol. 86, No. July, 1960. 


and Dept. Civ. Engrg., Univ. Waterloo, Waterloo, Ontario, 


Lecturer, Dept. Civ. Engrg., Univ. Manitoba, Winnipeg, Canada. 
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This paper presents complete analysis the elastic and inelastic be- 
havior simple, continuous, reinforced concrete beam. Plastic deforma- 
tions required for the realization ultimate load capacity are determined for 
both simple-plastic and strain-hardening behavior, and inelastic 
deformation reinforced concrete beams are established. Comparison de- 
amount plastic moment redistribution possible. The results number 
experiments intended test the theory are also given. Although the case 
studied, propped cantilever beam, has only limited direct usefulness, the re- 
sults may given application the design continuous reinforced 
concrete beams. 

Notation.—The letter symbols adopted for paper are defined where 
appear and are arranged alphabetically, for convenience reference, 
the Appendix. 


ANALYSIS BEHAVIOR FIRST YIELD 


Consider propped cantilever span shown Fig. single con- 
centrated load, applied point distant from the fixed end. The mo- 


ment diagram shown Fig. and completely described and M2, 
the moments sections and 
For elastic behavior, My; 


and: 
Let 


the ratio the yield capacities the two critical sections, and consider only 
<1.0 


Then for 
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CASE 


FIRST YIELD 


SECTIONS 


REDISTRIBUTION 


ULTIMATE 


FIG, 1.—LOAD-MOMENT RELATIONSHIP (IDEAL) 
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first yield will occur section For yield will occur simul- 


taneously and and collapse mechanism will develop immediately. 
Fig. may used locate the section first yield when and are 
known. Fig. the ratio loads for first yield section and sec- 


tion 


plotted against for various values and 


ANALYSIS BEHAVIOR AFTER FIRST YIELD FOR IDEAL ELASTIC- 
PLASTIC RELATIONSHIP 


Moments and already been pointed out, yield may first occur 
These alternatives are termed “Case and “Mechanism” respec- 
tively. The behavior first yield is, course, entirely elastic. The be- 
havior after first yield may seen considering Case first yield the 


(See Fig. 1). the load increased into the range 
there plastic rotation section 1(under constant moment while 
increasing towards order that plastic collapse mechanism can form 


essential that the plastic rotation capacity section adequate per- 
mit reach When section yields, the ultimate capacity 


the structure, and further deformation will occur under constant load until 
one the “hinges” destroyed, crushing the concrete. Case be- 
havior will similar Case the case termed Mechanism both critical 
sections yield simultaneously and plastic deformation required since the 

show the variation the moments the critical sections with increas- 
ing load, they are expressed ratios This yield moment has been 


used define andis, definition, the largest yield moment capacity 
the beam. Loads will expressed ratios Using such dimension- 
less ratios facilitates the analyses and leaves the expressions functions 
two parameters, and 
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FIG, 3.—FIRST YIELD LOAD RATIO VERSUS LOAD POSITION 
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FIG, 2.—FIRST YIELD LOAD VERSUS POSITION LOAD 
Ry 
| 
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Case I.—From considerations statical equilibrium the ultimate load for 
ideal plastic behavior 


From Eqs. and 


(7) 
P.. K(3-K)(m+1-K)° © @ @ re 2 


The moment when yield just reached 


The loads the range are given 


The nature plastic redistribution can seen from Fig. where values 


for the critical sections are shown against for 0.4 and various 
Pult 

values The dotted lines Fig. show first yield moment ratios, and 

versus (Eq. 7), and indicate how these 

approaches unity—when the critical sections yield simultaneously. More gen- 


eral variations are shown Figs. and which and are given fora 
range values ult and 


Case analysis similar that just described, may 
shown that 


and the ratio moments for section when yield has just been reached 
section 


eid! 


LIMIT DESIGN 
The load ratio the redistribution range is: 


Curvatures and Rotations 
Case I.—The curvature due moment first yield 


which the flexural rigidity the section. 
After first yield the left hand part the beam, AB, undergoes rigid body 
movement with constant moment its right end, and additional load must 


1.0 


FIG, 4.—MOMENT RATIO VERSUS LOAD RATIO 


“carried” the right hand part the beam, The plastic hinge 
separates the two elastic parts the beam, and spread over some finite 
distance called the hinge length, Considering free body diagrams the 


0.4 
08 
a2 0.4 0.6 0-8 1.0 


for) 
> 
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two portions the beam, Fig. can seen that now behaving 
cantilever with the differential load (Pp Py,) its free end. The incre- 


ment curvature the hinge, over that first yield, will the total angle 
rotation the angle discontinuity between the two elastic parts the 
beam spread over the hinge length. 


From Fig. the increment curvature section beyond that first 
yield 


which the total rotation the hinge, and and are the respec- 
tive angles rotation the left and right portions the beam. 
Considering only the incremental load (Pult acting the cantilever 


and taking new datum, the angle between the tangent the canti- 
lever and the new datum, 


The deflection the end the cantilever, the incremental deflection 
section over that first yield, 


the angle through which rotates undergoes rigid body movement 


Combining these gives 


(18) 


The curvature ultimate load equal the curvature first yield (Eq. 
13) plus the incremental curvature ultimate (Eq. 18): 
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Rewriting terms and gives the rotation factor: 


Case curvature due moment first yield 


After the first hinge has formed section the beam acts simply 
supported, Fig. Proceeding before with new datum for incremental 
loads only, the angle discontinuity the hinge the angle between the 
tangent beam and the new datum: 


The increment curvature the hinge 


which Hy, the hinge length section for the propped cantilever. The 
curvature ultimate load equal the curvature first yield (Eq. 20) plus 
the incremental curvature ultimate (Eq. 22): 


The rotation factor then 


Fig. very useful design tool. The left half this diagram represents 
Case failures. The rotation factor, given Eq. 19, plotted against the 
ratio first yield ultimate load, given Eq. for common values and 


Similarly the right half represents Case failures, with Eq. plotted 
against Eq. 10, inverted. 
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FIG, RELATIONSHIP 


RELATIONSHIP 


IDEAL 
CASE 


Case I.—At first yield, elastic theory, the deflection section 


The total deflection under the load ultimate found, from Eqs. and 24to 


Case first yield, deflection section 


(4-K) 


From Fig. the incremental deflection for the simply supported beam ulti- 
mate load is: 


The total deflection under the load ultimate found from Eqs. and 


2,3 


ANALYSIS BEHAVIOR AFTER FIRST YIELD FOR ELASTIC- 
STRAIN HARDENING RELATIONSHIP 


stated previously, the expressions for loading conditions beyond first 
yield are dependent upon assumptions made regarding moment-curvature re- 
lationships. The group design equations and charts just developed are based 
upon ideally simple elastic-plastic moment-curvature relationships. more 
precise, though more complex, treatment will now given, based upon 
consideration elastic-strain hardening moment-curvature characteristics 
(Fig. 7). 

Behavior first yield elastic and the expressions previously derived 
will apply. the loading range, the load the formation 


collapse mechanism, there yielding tension reinforcement section 


LIMIT DESIGN 


Increasing deformations will produce slightly higher resisting moments. Sec- 
tion behaves elastically with increasing load until reached, when sec- 
tion reaches its yield moment capacity. Behavior may such that section 
reaches ultimate resistance first yield occurs section which case 
The load which the mechanism forms then taken the ul- 
timate load capacity the structure. Subsequent behavior may such that 
section has additional rotation capacity and reaches its ultimate moment 
capacity, while section may may not reach its ultimate moment capacity 
ultimate load (Fig. 10). For Case there are parallel stages behavior. 


FIRST YIELD AT SECTION | FIRST vieto AT SECTION @ 


9.—ROTATION FACTOR VERSUS FIRST YIELD LOAD RATIO 


Moments and Loads.—It will found convenient express the moments 
and loads ratios My, and Py, respectively, the expressions de- 
veloped. 


Case I.—Up first yield the usual elastic equations apply. the forma- 
tion the mechanism 


ultimate load, ultimate failure section assumed. However, there are 
three possible conditions moment From Fig. 
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CASE 


FIRST YIELD SECTIONS 


ULTIMATE 


FIG, RELATIONSHIP (NON-IDEAL) 


The following parameters are then defined: 


ult, ult, 


For (i), where ultimate moment section develops section reached 
first yield, the load ratio is: 


ult 


For (ii), section has greater rotation capacity before ultimate failure, per- 
mitting the moment section reach the value MEP, 


The load ratio first yield then 


(33) 


and the load ratio ultimate is: 


For (iii), both sections reach ultimate moment capacity together, and the load 
ratio first yield given Eq. 33. The load ratio ultimate 


Considering now the moment ratios, Eqs. and apply first yield 
the case simple elastic-plastic behavior. the formation the mechanism 
the moment ratio section proportion from Fig. 11: 


Case first yield the elementary elastic equations apply. the for- 
mation the mechanism derivation similar that for Case gives: 
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FIG, VERSUS LOAD RATIO (NON-IDEAL) 


CASE 
1.0 
| | 
FIG, 
CASE 
1.0 
Ms 
My, 
P 1.0 


ultimate load, ultimate failure section assumed. 


There are three 
conditions moment section 


The various load ratios are: 


(i) 


= “eee . . (39) 


(ii) first yield: 


and ultimate: 


ose eee eee eee 41) 
m+m, 


and for (iii), first yield Eq. applies, and ultimate: 


1-K 


The moment ratio first yield found from Eq. 11. the formation 
the mechanism, the moment ratio section (from Fig. 12) 


Curvature.—The curvature relationships developed are for the criti- 
cal hinge, which are the first form, and are expressed terms the curv- 
ature first yield. 

Case I.—At first yield Eq. applies. For the simple case the loading 

range Py, the moment section remained constant Now 


there increase the moment from My, some intermediate-value 


The moment section increases value when the mechanism 
first formed. 


Consider the beam action Fig. made two parts: 
one, the pure cantilever, for the ideal case, with the incremental load 
its free end, and the other, cantilever beam with the incre- 
mental moment (MEP, applied its free end. Since the com- 
pletely elastic part the beam the right the two actions can 
superimposed achieve the actual behavior 

The effect this incremental moment decrease the angle discon- 
tinuity well decrease the deflection the elastic part 
the beam the left the hinge, subjected increased curvature the 
moment section increases. The hinge length increases with increasing 


load, and its value can determined from study the bending moment dia- 
gram. 


the formation the mechanism, proceeding the simple case,the 
expressions for angles and deflections are: 


(44) 


(46) 


and the incremental curvature 


The curvature associated with the strain-hardening portion 


the non-ideal, moment-curvature relationship Fig. can thought 


and purely plastic hinge curvature. 


The curvature ratio the formation the mechanism, expressed terms 
the parameters, is: 


LIMIT DESIGN 


ultimate load, the curvature ratio becomes: 


which 1.0 <mg 


Case first yield Eq. applies. range <P< there 


increase the resisting moment the first hinge. Considering only the 
increment load beyond that causing first yield, the beam behavior can 
separated into two parts. First, for the simple plastic case, the simple 
beam with the differential load this must added the effect 


the change moment the hinge, given consideration simple beam 
with the differential moment atits right support. The effectis 
course reduce the curvature and deflection the beam. 


the formation the mechanism, the expression for angle, curvature and 
curvature ratios are follows: 


3EI 


ultimate load the curvature ratio is: 


Case I.—At first yield Eq. applies. the formation the mechanism the 


total deflection under the load found combining Eq. and Eq. 45: 


(54) 
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ultimate load, the total deflection under the load for which 
expression may obtained substituting and MEP, Mult; 


Eq. 54. 

Case first yield Eq. applies. the formation the mechanism, 
the incremental deflection under the load the the simple beam 
with load less the deflection the simple beam due the moment 


The total deflection under the load then: 


4 - - - - - 


ultimate load, the deflection under the load for which ex- 
pression may obtained from Eq. substituting and 


MOMENT AND ROTATION CAPACITIES REINFORCED CONCRETE 


Analysis moment and rotation capacities are well established for rein- 
forced concrete.3,4,5,6,7 The following treatment based the usual as- 
sumptions: 


(1) that concrete stress distribution may described the parameters 
shown Figs. and 14, 

(2) that concrete has tensile strength, 

(3) that strains vary linearly across section, 

(4) that the steel reinforcement perfectly bonded the concrete, 

(5) that the stress-strain relationship for the steel known, 

(6) that the concrete will crush when strain reached, and the dis- 
cussion course limited under-reinforced sections. 


“The Ultimate Load-Theory Applied the Design Reinforced and Prestressed 
Concrete Frames,” Baker, Concrete Publications, 1956. 

Ultimate Strength and Deformation Plastic Hinges Reinforced Concrete 
Frameworks,” Chan, Magazine Concrete Research, Vol. No. 21, Nov- 
ember, 

“Plastic Hinging the Intersection Beams and Columns,” Ernst, ACI 
Journal, Vol. 28, No. 12, June, 1957, pp. 1119-1144. 

“Ultimate Slopes and Deflections—A Brief for Limit Design,” G.C. Ernst, Pro- 
ceedings, ASCE, Vol. 81, No. 583, January, 1955. 

“Elastic-Plastic Design Single Span Beams and Frames,” H.A. Sawyer, Pro- 

ceedings, ASCE, Vol. 81, No. 851, December, 
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Tension Reinforcement the usual elementary elastic theory 


which 


TRANSFORMED STRESS 
AREA DISTRIBUTION DISTRIBUTION DISTRIBUTION 
(a) 


FIG, 13.—TENSION REINFORCEMENT 


Cy: ASS Ce 


Te As fa 
TRANSFORMED STRESS STRAIN STRESS 
AREA DISTRIBUTION DISTRIBUTION DISTRIBUTION 
(a) 


FIG, 14.—TENSION AND COMPRESSION REINFORCEMENT 


and the other quantities are defined Fig. 13. 
limit 


The curvature the elastic 


parameter describing the effective percentage steel, 
nls | J | 


Analyses have shown that the maximum ultimate re- 
sistance reinforced concrete beam can written 


which Kj, Kg, and are parameters used describe the concrete stress- 
strain curve and the steel stress failure (Fig. 13). The value fsy 
ultimate depends largely upon the stress-strain relationship for the rein- 
forcing steel. For mild steel with stress-strain diagram shown Fig. 


developed Gaston. denotes the strain. 


For stresses the strain-hardening region, when and 
the case strain hardening material, the stress is: 


which Egy the modulus elasticity for the steel the strain hardening 
region and indicated Fig. 15. The parameters these equations 
were studied Gaston® and later Hognestad, Hanson and 
The most complete results indicate that 


3900 0.35 


The curvature ultimate 


Eate 


“An Investigation the Load-Deformation Characteristics Reinforced Concrete 
Beams the Point Failure,” J.R. Gaston, Siess, and Newmark, 
Studies, Univ. No. 40. 

Stress Distribution Ultimate Strength Design,” Hognestad, 
Hanson, and McHenry, ACI Journal, Vol. 27, No. December, 1955, pp. 455-479. 


INDEX 


REINFORCEMENT 
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MILD STEEL WITH 
STRAIN HARDENING REGION 


TYPICAL RELATIONSHIP 
MILD STEEL 


TENSION STEEL STRAIN. 


FIG, 15.—TYPICAL REINFORCEMENT CURVES 


IDEAL STRAIN HARDENING 
MATERIAL 


where and are strains the steel and concrete, respectively. This, 
together with the previous expressions would enable and their ratio 
determined, though perhaps not readily. For use plastic design, 
simpler approach essential. 

From the assumption linear strain distribution 


and for equilibrium forces the section (Fig. 13b) 


that 


Introducing this into the expression for gives more tractable expres- 
sion: 


and the rotation capacity then 


which may very considerably simplified noting that 


thus giving 


MAX 


Pf, 
A a 
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For design, Eq. may approximated 


Combined Tension and Compression Reinforcement.—In this case 

has developed expressions for moments resistance, stresses and strains 

which, though similar those just described, are rather more complex and 
need not reproduced detail. The moments resistance are (Fig. 14) 


and p’, refer compression reinforcement. Using the same pro- 


cedure the case simple tension reinforcing, approximate expres- 
sion for the curvature ratio may established: 


TESTS 


Outline.—A test program limited scope was undertaken supplement 
the analytical parts the study. tests were made all. Two simply- 
supported single-span beams and six symmetrical two-span continuous beams 
were tested. 

The beams all had in. in., and the span lengths 
were 4ft. Properties the various are given indetail 
Intermediate grade hi-bond reinforcing steel was used throughout. Tension 


tests the steel showed characteristic stress-strain relationship without 
yield plateau, but rather with tendency commence strain-hardening im- 
mediately after yield. Fig. shows typical stress-strain curve for the re- 

inforcing steel together with average characteristics for the different sizes 


used. High early strength Portland cement was used fabricating the beams. 
Stress-strain characteristics for typical concrete speciman are also shown 
Fig. 16. Web reinforcement was designed resist full shear ultimate 


load and consisted simple vertical stirrups made from 1/4-in.- diameter 
plain mild steel bars. 


and 
where 
pf - 


July, 1960 


TABLE 1.—PROPERTIES OF BEAM CROSS SECTIONS 


0.0068 


0.0140 
0.0141 


0.0278} 0.0208 


44.6 
0.00193 
78.9 


0.007 


FIG, 16.—STRESS-STRAIN CURVES FOR REINFORCING STEEL AND CONCRETE MIX 
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The test beams and associated concrete cylinders were moist-cured until 
the day their test. Stress-strain the concrete cylinders were ob- 
tained with extensometers. Inthe tests the continuous beams particular care 
had taken avoid the introduction secondary stresses due uneven 
supports. the beam tests the total applied load was measured, and the cen- 
tral reaction was weighed with very stiff load cell. Strains were measured 
with gauges the concrete, and exposed portions the reinforcing steel. 

Deflections were measured with direct-reading dial gauges. 

Results.—A detailed study was made the behavior and mode failure 
each beam. Analyses were made show the variation with load strain, 
curvature, deflection, and position neutral axis. Figs. and show the 
results such analyses for beam number and conjunction with the char- 
acteristics the materials shown Fig. permit precise determination 
the mode failure. Figs. present comparisons theoretically 
derived and observed moment and load ratios for the continuous beams all 
stages failure. provides additional comparative experimental and 
theoretical information. 

particular problem the study was the determination the flexural 
rigidity the various cross sections. Computed values curvature ratios 
and deflections are, course, largely dependent the accuracy the de- 
termination. The flexural rigidity may computed from either the product 
the transformed section moment inertia and the initial tangent modulus, 
the gross effective moment inertia d3/12). Comparing these values with 
experimental determinations obtained from initial slopes moment-curvature 
diagrams showed considerable scatter between computed and observed results. 
The analyses given are based entirely the simple gross moment inertia 
combined with the initial tangent modulus elasticity for the concrete. Al- 
though less than entirely satisfactory this simple approach for the deter- 
mination curvature ratios and deflections, and appears good any 
method available. 

Examination the results Figs. and Table show satisfactory 
agreement between derived and observed results. appears that the theory 
based upon elastic-strain hardening behavior more accurate than that based 
upon ideal elastic-plastic behavior. Results from both theories appear 
conservative, and the simpler theory would appear entirely satisfactory 
basis for design. The ability estimate deflections first yield seems 
fairly good, though there was tendency underestimate ultimate deflections. 
This due part tothe ultimate larger than anticipated. The use 
estimating required curvature capacities. 


APPLICATION THEORY 


The analyses presented for the propped cantilever with single 
load and the design expressions and charts developed using the ideal elastic- 
plastic, moment-curvature relationship represent the complete solution one 
part the general problem reinforced beam. Similar 
design expressions and charts can evolved for other continuous loading 


uniformly distributed varying, for partially loaded spans and for other com- 
binations loadings. 
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VARIATION WITH LOAD 


FIG. 17.—CONCRETE AND STEEL STRAIN, CURVATURE 
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FIG, 24.—MOMENT RATIO VERSUS LOAD RATIO 
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FIG, 25.—TANGENT MODULUS, VERSUS CONCRETE STRENGTH 
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FIG, 26.—CALCULATED VERSUS EXPERIMENTAL FLEXURAL RIGIDITY 
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continuous beam can considered made propped cantilevers 
(the exterior spans) and fixed partially fixed beams (the interior spans). 
parallel set design expressions and charts forthe most im- 
portant loading conditions encountered the fully fixed and partially fixed con- 
tinuous beam segments. The design expressions for the fixed beam with single 
concentrated load are presented elsewhere.10 Ineach particular case load- 
ing, the design expressions charts would present complete solution forall 
the variables and for the particular continuous beam segment considered. 

expressions for each beam segment and loading type considered 
include moment-load relationships that can written from considerations 
elasticity and the collapse bending moment diagram. Aconsideration ductility 
required ensure that there sufficient rotation capacity each the 
early forming hinges order that the collapse mechanism canbe formed. This 
requires aknowledge the order which are formed. For the case 
the propped cantilever with single concentrated load this simplifies the 
critical section which first yield will occur. The rotation capacity required 
function the geometry span length and position load well 
the parameter the ratio yield moments the critical sections. addi- 
tional consideration apparent importance reinforced concrete safety 
factor against first yield critical section. Flexural rigidity decreases 
quite rapidly yielding the tension steel takes place. From the point 
view service working conditions, working loads sufficiently be- 
low the yield value prevent any sudden reduction flexural rigidity. There 
should also some margin safety against first yield occurring critical 
sections under working loads. 

The spans, load positions and values working ordesign loads are normal- 
known continuous beam design problem. The ultimate loads that the 
structure will designed for can arrived the use suitable load 
factors. The value the parameter can chosen asto give any desired 
modes failure. Thus, making use Figs. and either Case Case 
Mechanism failure can produced for the propped cantilever that was ana- 
lyzed. 

From consideration Fig. similar diagrams, the required values 
the rotation factor and the ratio first yield ultimate load are deter- 
mined. The first yield load can then computed well the safety factor 
against first yield. The rotation capacity required can computed from the 
rotation factor. For design the hinge length can conservatively taken 
being equal the effective depth the beam. The moment and rotation ca- 
pacities for singly and doubly reinforced concrete cross sections have been 
presented and are used proportion the critical sections for adequate yield- 
moment capacity, and ductility rotation capacity. 


APPENDIX—NOTATION 


depth transformed area. 


width transformed area. 


depth from compression face beam center longitudinal 
tensile reinforcement. 


“Application Plastic Design Theory Reinforced Concrete Continuous Beams,” 
Berwanger, Sc. Thesis, Queen’s Univ., 1957. 


=, 
7” 


w 
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1960 
modulus elasticity. 


modulus elasticity strain hardening region. 
steel stress failure. 


unit compressive stress concrete. 
hinge length. 


moment inertia. 


increment length parameter. 
parameters used describe the concrete stress-strain curve. 
span length. 
moment force. 


ratio yield capacities critical sections. 
point force. 


load the formation collapse mechanism. 
parameter describing the effective percentage steel. 


deflection. 
deformation. 
total rotation hinges. 
angle rotation left portion beam. 


angle rotation right portion beam. 
curvature. 
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SECONDARY STRESSES PARALLEL WIRE SUSPENSION CABLES 


Thomas Alan Wyatt! 


SYNOPSIS 


The working stress parallel wire cable normally based the aver- 
age stress across section the cable. This paper analyzes the variation 
stress across the section, caused deflection the cable. The influences 
clamping bands and tensioned wire wrapping are considered. Numerical 


examples show that secondary stresses may important part the ef- 
fects live loading. 


INTRODUCTION 


The advantages the suspension bridge with parallel-wire cables for long 
spans are well known, and relatively short time considerable advances 
have been made the methods design and analysis. other fields this 
has been accompanied increase the permissible stresses. One the 
characteristics this type bridge that the live load has relatively little 
effect the tension the cables which implies that the range stress 
service low and also that adequate load factor the live load may 


Note.—Discussion open until December 1960. extend the closing date one 
month, written request must filed with the Executive Secretary, ASCE. This paper 
part the copyrighted Journal the Structural Division, Proceedings the Ameri- 
can Society Civil Engineers, Vol. 86, No. July, 1960. 

Asst. Engr., Freeman, Fox and Partners, Cons. Engrs., London, England. 
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achieved while using high value total stress. Another characteristic how- 
ever, low stiffness, and the deflections may therefore large, including, 
particularly, the change slope near the supports. this respect the live 
loading proportionately more significant, and thus, the stresses induced 
the cable due bending are enhanced interest the designer. 

When the cables are “spun” every effort made ensure that the stress 
each wire identical. Initially each wire behaves independently and, pro- 
vided the cable supported saddles suitable radius, the secondary 
stresses are small. The next stage construction fit the bands act 
attachments for the deck suspenders hangers. These clamp the wires 
together and this extra restraint leads more serious secondary stresses 
the dead load the deck built up. 

Usually, the spun cable wrapped with wire under tension, and although 
the transverse stresses set between the wires are small, the friction 
stresses thus induced deformation the cables are very important. The 
loads causing deformation after wrapping include all the live load and possi- 
bly some part the dead load. 

The most significant feature the deformation the cable respect 
secondary stresses the deflection angle, defined the change slope 
the line resultant cable tension. The deflection angle clearly uniform 
throughout the length between two adjacent hangers, and found that the 
actual slope the cable axis very nearly equal the slope the cable 
tension line except within very short distance from hanger band saddle. 
The angle between the cable axis and the tension line is, therefore, referred 
the local bending angle. convenient consider the corresponding 
local bending stresses separately from other secondary effects, although the 
effect the band clamping and wrapping the local bending angle must 
taken into account. 

The worst conditions secondary stress normally occur the tower, 
probably the side span the backstays are loaded. The total secondary 
stress then found combining the effects the deformation occurring be- 
fore wrapping the cable (band clamping effect) with the stresses due de- 
formation the wrapped cable, and including also the local bending stresses. 

Notation.—The letter symbols adopted for use this paper are defined 
where they first appear, illustrations the text, and are arranged 
alphabetically, for convenience reference Appendix 


LOCAL BENDING STRESSES 


The fundamental behavior tendon carrying large axial tensile stress 
and subjected bending one end shown The angle the 
rotation the tangent (say clamping band the saddle) relative the 
direction the resultant tension. Clearly any section 


from which 
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FIG, 2.—RESTRAINT CAUSED CLAMPING THE CABLE 


Typical band ypical band 


FORM CABLE WITH CLAMPING BANDS 


FIG, 1.—BENDING DEFLECTION SIMPLE TENDON 
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these equations the total tension the cable, refers the radius 
gyration, denotes the dimension measured along the cable and the ten- 


sile stress primary stress). 


Determination the constants from the end conditions shown Fig. 
gives 


where 


which refers the bending stress. will noted that and 
1/g. For circular wire radius 


For circular wire with tons per and 13,000 tons per 
that very small indeed relation the panel length. 

The deflection angle largest the end the span, and the local bending 
stress usually limited supporting the cable curved saddle radius 
say. The stress due individual bending wire the saddle clearly 


r 


will seen that this maximum value does not change with live loading. 
However, there variation stress taken into account when consid- 
ering fatigue effects, due the movement the point where the cable 
tangent the saddle. Normally this movement closely equal 


where the live load deflection angle. Thus, the farthest point touching 
the saddle, the bending stress varies from (r/R) (r/R). 


CABLE BAND CLAMPING EFFECT 


there were cable bands wrapping, there would resistance 
“slip” any wire relative those adjacent, and the deflections and stresses 
all the wires would identical. that case, anycross section the cable 
deflects without rotation. If, however, the cable clamped band that 
prevents slip, the effect the change slope the cable tension would 
exert couple the band shown Fig. 2(b). Some rotation the band, 


and the maximum bending stress (writing the depth the tendon 


therefore, occurs that the couple implicit uniform “primary” stress 
distribution reduced and balanced redistribution the axial stress, 
that is, moments due “secondary” stresses. 

Owing the load applied the cable the band, and (Fig. 2(c) are 
not equal, and include small angle. However, any practical case the dif- 
ference magnitude very small. The “primary” couple may written, 
the two sides the bands. That is, defined the net rotation the 
cable band relative the bisector the lines cable tension the adjacent 
panels. The effective length, the clamping band seen increased 
lengths due local bending. Thus, actual length band plus 

The analysis cable band calmping for unwrapped cable most com- 
monly required for the dead loading. the dead load usually almost uni- 
form throughout the span, the deformation uniform “curvature” with the 
deflection angle changing (n-1) each band, where the 
angular deflection the end panels. convenient write (n-1) 
two found that each step has linear load-deformation relation 
that the total effect may found superposition. 

Firstly, the end cross section the cable the tower made rotate 
tive rotation the two end sections every panel (including the end panels) 
cable tension. 

Secondly, the end the cable returned its true attitude applying 
equal and opposite rotation without further deflection the cable. The be- 
havior this stage shown Fig. where, for clarity, the direction the 
cable tension has been drawn horizontal throughout. The resultant action 
the stress panel represented the tension plus moment 

For the equilibrium band 


From the relative rotation the ends the panel p-1,p 


The equations this type for each panel may solved recurrence rela- 
tion. found practice that the values decrease rapidly away from 
the ends, and the influence one end the bending effect originating from 
the other end the span quite negligible. For example, the effective stiff- 
ness the cable the left any band near the right hand end may, there- 
fore, simply defined coefficient such that from which 


= G 95-1 
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the solution obtained 


For the first panel The maximum secondary strain 
(80 from which the secondary stress, including the sagging defor- 
mation the first step,is 


With the normal order numerical values the “damping” very rapid 
and the effect the first step (the term) small, that this solu- 
tion could also used for live load deformation unwrapped cable. 


SLIP THROUGH HANGER BANDS 


The difference between the stress adjacent panels each wire corres- 
ponding the change moment, supplied friction stresses action 
within the cable band length. assumed that due the clamping pressure 
there limiting friction stress, where P,, the clamping force 
per unit length each side the band and the coefficient friction 
within the band. Appendix that there thus limiting max- 

The coefficient friction, this expression difficult estimate 
because the effect the local constriction the cable increase the 
friction considerably above the normal value applicable when the transverse 
pressure uniform over large area (as for example wrapped length 
cable). The value has been suggested. 

there the band the maximum change moment (1-K)2 
the band nearest the tower. This value may, however, practical case, 


use parameter defining the “effectiveness” the band clampirg, such 
that 


The rotation relative the line tension section the cable the band 
farthest from the tower which slip occurs may written where 
clearly 

Treating the deformation two stages before, there tendency 
slip the first stage. the second step the moment applied rotate the 
end section through angle return its initial inclination. The mo- 
ment thus caused each panel the cable shown graphically Fig. 

The area this diagram corresponds the relative rotation the cable 
cross-sections. Thus, because the tower and the farthest which slip 
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Bands not Bands 


Position along Cable Tower 
Cable bands 


FIG, 4.—CABLE BAND CLAMPING EFFECT 


Whole Section Each Wire 


Stretching Stress 
Section Section Section 


FIG, 5.—BENDING AND STRESS DIAGRAMS FOR 
SIMPLIFIED TWO-WIRE CABLE 


Limit Sli 
Bending Stress 


the above expressions the number bands which slip occurs, 
and thus clearly integer. The range values for each value 
obtained investigating the condition that additional band about slip: 

The maximum moment due the end effect reduced from (1-K) 
from Eq. found almost independent (and thus the band 
clamping effect parameter, A), sufficient approximation for practical 
purposes 


The maximum secondary stress (compare Eq. thus, 


will seen that the reduction stress fairly small. For example, the 
clamping force one quarter that necessary prevent slip, 0.63. 


SECONDARY STRESSES DUE FRICTION BETWEEN WIRES 


clear from the diagram illustrating the cable band clamping effect 
(Figs. and that the deflection parallel wire cable that not wrapped, 
causes rotation any cross section relative the cable axis. Thus, there 
slip between adjacent layers wires. strands are compacted and 
bound together wrapping with wire under tension, transverse stresses are 
set between the wires, and the stresses due friction radically alter the 
stress distribution caused bending the cable. 

The effect wrapping treated causing possible uniform, limiting 
efficient friction 0.20 (see Appendix III). and are the diameter and 
tension forces per strand wrapping wire. The wire used for wrapping 
commonly about 0.15 in. diameter gauge) soft annealed condition. 
such acase the value example for cable radius in. 
tons per in.), the limiting friction shear stress 0.015 tons per in. 

The analysis best approached considering cable which only one 
slip plane allowed. For simplicity the cable taken square made 
two similar rectangular wires each size ag. The deformation due 
the rotation the end section through angle relative the line re- 
sultant tension shown Fig. 5(a). Slips occurs over length only, 
and the upper wire, section adjacent, Ml, and 

The deflected form the wires the length the slip may investi- 
gated similarly. The bending moment each wire distance from 


(9) 


where the deflection, since factors tending cause bending the two 

wires opposing directions must compensated for transverse reactions 

between the wires, because thedeflections are equal. These transverse forces 

are sufficiently small for their effect the friction stress neglected. 


Between and the mean stress each wire varies linearly due the 
friction stress shown Fig. The condition that the limit the 
length over which slip occurs that the strain the middle fiber the up- 


wire (relative the lower wire) corresponds the relative rotation 


whence, substituting for and rearranging terms 


Eliminating from Eqs. and 


This equation may solved for successive approximation. However, 
large number, first retaining the terms only, 


and further 


Thus great simplification possible when larger than approxi- 
mately 100. 


Examining the significance these approximations will seen from 
Eq. that the effect the exponential local bending terms limited 
very short regions the two ends the slipping length. From Eq. the 
slope the end the slip (B) =s/3 whereas for the major part 
the slipping region the slope uniform These angles are thus very 
small compared with the end deflection. Lastly, will seen that the great- 


est variation stress where slip occurs small compared with the 
stress changes along the slip. 


The secondary stress may regarded the sum the “stretching” 


stress (Fig. 5b) 


plus the local bending stress for each wire due rotation determined 
separately previously described. 


MULTI-WIRE WRAPPED CABLE 


These simplifications may now applied cable square section 

typical element this cable from the end shown Fig. 2rp 
the depth the region over which slip not occurring and, thus, between 
and the stress uniform. Since slip taking place the right equal 
and opposite forces act both faces layer and the stress remains 
the same section AB. The stress section CD, and thus layer (p) 
the slip surface between and 

The difference stress between layer (p) and layer thus 


2a, 


The corresponding difference between the extension the middle fibers 
these layers must equivalent the end rotation Thus, 


This equation applies for all values from (q-1). For thecenter layers 


FIG, 6.—BENDING MULTI-WIRE CABLE 
slip has taken 
place are shown 
solid 
FIG, 7.—TYPICAL ELEMENT MULTI-WIRE CABLE 
FIG, 8.—REDUCTION STRESS OMISSION WRAPPING 
FROM THE REGION ADJACENT THE TOWER 
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The solution these equations thus 


p=1 


lg = 
will noted that this the same overall length for the two-layer cable. 
The secondary “stretching” stress obtained the sum the changes 


stress from layer layer, starting from the central layers the cable. 


Thus, the maximum stretching stress, which occurs the outermost layer, 
given 


which 


Therefore, 


The error this approximation the sum required very small, being 
order when large. The “local bending” stress (Eq. 


Eq. allowance made for supporting saddle) must added the 
stretching stress fp7. 

check the accuracy this analysis may made considering the 
deflected form actually implied the stress distribution previously derived. 
The strain energy due bending the wires this deflected form may then 
compared with the total work done the cable due secondary effects. 
Terms not included the above analysis are due 


local bending originating the end the cable, 
(b) smooth curvature the cable throughout the length slip, and 
(c) local bending originating the end the slip. 


Item (a) allowed for separate computation local bending stress that 
then added the stretching stress. Item (b) corresponds strain energy 


— 
tor 
Thus, 


the bending the end the included the two-layer analysis and 
leading stress and corresponds approximately strain energy 


The work done the cable due secondary effects included the 


With typical values the order 1/200, that the accuracy 
the analysis may expected satisfactory. 


CIRCULAR SECTION CABLE 


Results are required, practice, applied the circular section 
cable. The analysis is, however, considerably more difficult because the trace 
the slip surfaces the cross section (more strictly defined the lines 
joining wires between which there slip and which thus form elastic 
layer) are curved, convex away from the horizontal diameter. The geometry 
the circle also more complicated, and because the two slip surfaces 
bounding any layer are not equal width, the stress changes throughout the 
length the layer. 

indication the difference between the rectangular and the circular 
sections can, however, obtained considering regular octagonal sec- 
tion, for which solution may obtained arithmetical integration. must 
assumed that the slip surfaces remain plane but this should not cause seri- 
ous error for the octagon. The results obtained are 


total length slip; 


where the depth the section, across flats. The stress thus very 
nearly the same for square section having the same number slip sur- 
faces (q). The length slip is, however, reduced roughly the proportion 
that simply predicted considering typical slip plane mid-way between 
the center the cable top. therefore suggested that Eqs. and 
may used with sufficient accuracy for the circular cable, taking for the 
size the equivalent square cable Eq. 15, 


INITIAL LENGTH CABLE LEFT UNWRAPPED 


The maximum secondary stress reduced length cable adjacent 
the tower for example) left unwrapped. This case may analyzed 


the same method. The total length slip (including the unwrapped length) 
given ly, 


pressed 


where the factor dependent principally the parameter shown 
Fig. the line drawn represents close approximation for any value 


before. found that the stress may ex- 


EFFECT REVERSAL DEFORMATION 


has been shown that the behavior the wrapped cable governed 


1 
Successive deformations cannot treated superposition and also 
clear that the stretching stresses not return zero removal the load 
causing deformation. the deflection angle decreased from its maximum 
value, slip the reverse direction occurs, starting from the end the cable. 
The change friction stress where the slip reversed and the previ- 
ous analysis may used substituting twice the normal limiting friction 
stress. the deflection varied from the deflection angle for the 
reverse slip and can thus shown that the reverse slip just cov- 
ers the initial slip, and that the bending stress also exactly reversed. 


CABLE BAND EFFECT WRAPPED CABLE 


The rotation the end cross section the cable, relative the cable 
axis, has been considered without allowance for deformation applied other 
points the cable for the clamping effect the bands. many cases 
adequate description the secondary effects the tower thus obtained. 
However, the predicted length slip longer than about two panel lengths, 
modified method analysis must used. 

typical panel subjected bending the ends causing slip the same 
direction shown Fig. For simplicity analysis the stresses, bending 


moments, and deformation angles are written non-dimensional form, that 
is, respectively, 


The boundary the slip zone defined ag, da, shown, that and 
are also non-dimensional. The definition such that the previous 


| 
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Z ag ag 
(very nearly). Thus, for the point 


(24) 


the central region, slip extends throughout the panel, that (for the 


square cable) the stress any wire constant. example is, the limit- 
ing layer passing through 


and, using Eq. 24, 


The stresses outside the central core the left and right hand ends the 
penel respectively (in the upper part the cable) are approximately 


These expressions not apply for very small but the error introduced 
the subsequent moment integrals negligible. 


The total moments corresponding these stress distributions are 


where and are functions only 


These equations completely describe the relation between moments and 


deformation typical panel (p). The relation the next panel the right 
(p+1) 


> oy = pMR + 2 


assuming slip takes place through the band, where the change 
angle between the two panels due the load applied after wrapping. 

Angles and clearly decrease away from the tower. Provided the 
friction sufficiently large, relation the deformation originating 
from the hanger loads (that is, compared with there will sections 
cable over which slip occurs the panels remote from the tower. De- 
noting the panel nearest the tower for which this true the 
bending through angle the end nearer the lower not affected de- 
formation the cable farther from the tower. 

may obtained from Eq. and also from Eq. 29. The value 
for panel found solution Eq. and from Eq. 27. The pro- 
cedure repeated succeeding panels until value obtained which 
exceeds the given end deflection for which the required. The 


Movement 


Boundary 
slip Slipping zone 


Vertical scale 
greatly exaggerated 


FIG. 9.—LENGTH WRAPPED CABLE BETWEEN TWO CLAMPING BANDS 


choice the starting value may then improved obtain the true 
solution. 

The slip regions the bands that are remote from the tower are in- 
creased the deflection angle these bands increased relation the 
friction stress. The method described cannot strictly applied when there 
starting section which the cable entirely unaffected slip. The 
analysis panel subjected secondary bending opposite directions 
the ends becomes very difficult the deformation large that the two 
zones, slipping opposite senses, meet, because the slip boundaries cannot 
then simply defined. found, however, that the relation between stress 
and deformation the tower very little affected conditions three four 
panels distant, would anticipated from the exponential decay” nature 
the disturbance. analysis started from the panel nearest the tower using 
trial value for found that the moments displacements predicted 
working away from the tower successive panels become absurd unless 


Band 


(and thus the maximum stress) correctly chosen within extremely nar- 
row range. is, however, easier towards the tower, and, the prac- 


tical case, any reasonable value can combined with the trial value 


The maximum stress obtained from the previous equation for modi- 
fied apply the top wire layer. Thus, 


which may written 


Thus the increase stress due the cable bands given which 
shown plotted against Fig. for various values the band clamping 


Clamping parameter 


Stress function 


FIG, STRESS WRAPPED CABLE 
DUE BAND CLAMPING EFFECT 


cable bands the dotted line, which shows when the concave- 


upwards deflection angle each band 1/10 the end deflection down- 
wards, 


which 
Pig. 
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NUMERICAL EXAMPLES 


The relative importance the secondary effects may shown refer- 
ence typical long-span bridge with main span 3,000 and suspended 
side spans 1,200 each.. With two in. cables made 0.194-in.- 
diam wire, the quantities governing the secondary stresses may taken as: 
the panel length, 720 in. 

The analysis will applied the side span, for end deflection due 
dead load 0.050 radians (after fixing the clamping bands) plus deflections 
after due live loads and temperature changes 0.020 radians 
downwards 0.015 radians upwards. 

Cable Band Clamping the band length, in., the clamping 
parameter 1.1 and 0.255 (Eq. 5), the maximum stress then 
tons per in. (Eq. 6). The total clamping force necessary 
prevent slip 310 tons, determined placing Eq. (for 
likely that the actual clamping force would this order, 
since the bands are normally designed for conservative value If, 
however, the clamping force reduced half, 0.5 and 0.8, thestress 
reduced tons per in. (Eq. 9). 


Wrapped Cable.—In normal practice the wrapping tension such that the 
limiting friction stress the range 0.005 0.03 tons per sq. in. For 
the live load deflection 0.020 radians the respective stretching stresses 
would the range tons per in. (Eq. 19). 

Protection from corrosion the region where the wires are tangent tothe 
saddle may obtained using independent articulated sheath, for ex- 
ample the Golden Gate bridge. the wrapping omitted for distance 
the stress reduced from tons per in. (Fig. and Eq. 22). 
The reduction naturally less the case 0.005 tons per 
in. 

The stresses may also calculated the cable wrapped before appli- 
cation the dead load; that is, radians. For 0.005 the length 
slip then large, that the effect the cable bands must taken into 
account. From Eq. 23, 29, that with 1.1 and the 
stress (excluding local tons per in. and Eq. 
31), compared with tons per in. previously given. The in- 
crease seen quite small. The final value used for the friction stress, 
must take account the loss tension the wrapping caused the lat- 
eral contraction the cable under load, but the magnitude this reduction 
rather uncertain. 

Local Bending Stresses.—The fundamental equation for local bending gives 
1440 tons per in. (Eq. 2); but the cable bears saddle 
radius ft, the stress limited tons per in. (Eq. 3). will 
noted that 1/g 0.9 in. that the bending stress reduced 
factor each in. along the wire. Local bending may also signifi- 
cant the first cable band. the previous example the rotation the first 
band, slip occurs, 0.013 radians, which, together with the pos- 
sible coexistant local deformation due the load that band, could cause 
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CONCLUSIONS 


Secondary Stresses due Dead Load.—These examples show that the dead 
load stresses are quite large. The stresses would considerably reduced 
over 40% the example given) the band nearest the end the span was 
reléased after the remaining bands were loaded. Some relief may occur 
this way due slip caused the loss clamping force the cable con- 
tracts laterally loading: the bolts clamping the bands are normally re- 
tightened after most the load has been applied. similar reduction would 
obtained the erection procedure was such that the end band was actually 
fitted while the cable was deflected towards its final position. The local bend- 
ing effect the end band may also counteracted relieving radius sim- 
ilar function the saddle radius. 

Secondary Stresses due Live Load.—The live load secondary stresses 
may additional importance, since, due repeated application, their ef- 
fect the ultimate strength the cable may relatively greater than their 
actual magnitude, whereas the converse would normally true the dead 
load stresses. For the example considered, taking typical value for the 
per in. (primary and total secondary stresses respectively). The change 
“stretching” stress due the change deflection from 0.020 radians down- 
wards 0.015 radians upwardsis tons per and the maximum co- 
existant change local bending stress very nearly tons per in. (by the 
saddle). Allowing also for tons per in. change primary stress, the 
drawn wire 100 tons per in., ultimate tensile stress, this corresponds 
factor safety the order 1.5 against failure due fatigue. Alterna- 
tively, considered terms load factor, applying factor 1.1 the dead 
load primary stresses (which are predictable high standard accuracy), 
the factor live loading approaches 2.0 even after application addi- 
tional “factor ignorance” 1.5 applied all the secondary stresses. 
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APPENDIX I.—NOTATION 


radius cable. 
half length side equivalent square section. 
Youngs Modulus. 

tensile stress (primary stress). 


bending stress (secondary stress). 


1960 


flexural stiffness one panel cable. 


Axial force 
ratio net area gross area cable section. 
radius gyration. 


dimension measured along cable. 


panel length: that is, distance along cable between hangers. 
span 


number panels. 


radius individual wires cable. 
radius saddle ends span. 
total tension cable. 


deflection angle. 


deflection angle the end panel. 

difference between the deflection angles two adjacent panels. 
local bending angle. 


deflection angle the end panel due load applied after wrapping. 


rotation cable band relative tension line (see definition text 
description Fig. 3). 


Band Clamping Analysis. 


factor which effect diminishes each panel. 
effective length clamping band (Fig. 3b). 
clamping force per unit length each side band. 
Tle 


coefficient friction within band. 
parameters for effect slip through bands. 


Wrapped Cable Analysis. 


dy, diameter and tension force per strand wrapping wire. 
stress factor for effect clamping bands. 
length slip. 


length slip simple two-layer cable. 
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length cable adjacent tower without wrapping. 
equivalent limiting shear stress. 


stress factor for effect unwrapped length tower. 
coefficient friction wrapped cable. 


APPENDIX 


SLIP THROUGH HANGER BAND 


The difference between the stresses the wires adjacent panels sup- 
plied friction within the cable band length. Considering the forces the 
wires indicated stippling Fig. 11(a), the change tension between the 
ends the band denoted due change the total cable 
tension and the secondary stresses respectively, and equal the fric- 
tion force 

The total change primary tension sin (see Fig. 11(b) where 


the hanger load and sin the inclination the cable the horizontal. 
Thus, 


The change moment and the stress varies linearly 
across the section, 


assumed that clamping the band causes uniform transverse pressure 
that the limiting friction stress may written Thus, for 
slip occur the boundary indicated (that is, with plane lower surface) 


The critical boundary determined the condition =0. normal 
practice the factor safety (C) against slip due primary forces alone 
quite large, and good approximate solution then 


which 


C =27 le ue P,/W sin 


This value slight over-estimate because slip would occur first 
boundary with the convex downwards, but the error small since 


Band axis 


Cable 
(b) 


FIG, FORCES BETWEEN WIRES CAUSED WRAPPING, 


- 
FIG, ACTING PART CABLE WITHIN CLAMPING BAND 
! 
(b) 
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the boundary shown region where the stress change each wire 

the applied moment, greater, some increase can still occur. 
Considering the change moment the force transferred across neutral 
axis (m) times the lever arm, the stress distribution were remain linear 
with the maximum possible force transferred across the axis 


Thus, reasonable approximation the behavior the band given the 
simple relation, that under increasing deflection the cable the moment 
change the band increases linearly with increasing band rotation, toa 
limiting value and thereafter remains constant. 


LIMITING SHEAR STRESS WRAPPED CABLE 


Fig. (a) shows cable wrapped with single layer wire diameter 

ty. The tangential tension force per unit length cable 

thus The ideal wire packing shown Fig. (b). 
Considering the vertical section AA, and resolving horizontally, 


where the number wires the section. 


The limiting longitudinal friction force whence the 
limiting stress 


Similarly for the section BB, 


The cable with ideal packing thus slightly anisotropic this respect. 
the practical cable there are also variations the forces between the wires 
due irregular packing. may noted that the distance between successive 


planes type (on which maximum) larger than between planes 
type AA. 
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VARIOUS INSTABILITY MODES THE FIXED BASE COLUMN 


Sawyer,! J.M. ASCE 


SYNOPSIS 


Equations are presented for the critical buckling load column, fixed 
its base and loaded its upper end, under conditions. The case 
the column loaded its upper end through frictionless roller believed 
presented here for the first time. column loaded its upper end through 
pinned link has appeared rather obscurely the literature. This case 
developed and expanded for comparison with the similar case loading through 
The columns lift-slab operation are loaded hanging links. 
Equations that have possible application this field are derived. 


INTRODUCTION 


The fixed base column exhibits various modes instability depending 
the loading and support conditions the upper end. Two classic examples are 
(1) the column fixed the base and simply supported the upper, loaded end, 
and (2) the column fixed the base and free translate the upper, loaded 
end. Solutions both these cases have been known since Euler’s time. 
well many others, has presented complete derivations 


the controlling differential equations. The first these problems leads 
transcendental equation, 


Note.—Discussion open until December 1960. extend the closing date one 
month, written request must filed with the Executive Secretary, ASCE. This paper 
part the copyrighted Journal the Structural Division, Proceedings, the Ameri- 
can Society Civil Engineers, Vol. 86, No. July, 1960. 
Asst. Prof., Dept. Civ. Engrg., Univ. Florida, Gainesville, Fla. 


“Theory Elastic Stability,” Timoshenko, McGraw-Hill Book Co., Ind., New 
York, 1936. 


which has the smallest root, 4.493. The resulting primary buckling load 
is: 


The second problem has the solution its differential equation the relation- 
ship: 


Por 


ling load is: 


design-office practice, the situation sometimes encountered where the 
use one the other these basic buckling loads ostensibily required. 
One such situation that involving the columns lift-slaboperation. Another 
case might arise design problem where advantageous necessary 
support the expansion end bridge girder, similar structure, ona 
free-standing column. probable that some designers would consider the 
design quite conservative, the columns either condition above were pro- 
portioned using the buckling load given Eq.4 and modified some reason- 
able factor safety, such Actually, the resulting design may either 
grossly conservative dangerously light, depending some rather subtle 
differences the construction details. 

first consideration, the free-standing column supporting the expansion 
end structure will examined. considered that the far end the 
girder structure either simply supported otherwise fixed against trans- 
lation. The expansion detail may simple sliding plate or, more generally, 
detailed sector circle. Conceivably, some designers might consider 
simple link the equivalent the sector circle because many situ- 
ations the action link anda roller are almost identical for small move- 
ments. poorly designed “roller” such that shown Fig. might used 
here. The sliding plate will not considered further and attention will 
focused the roller link designs. 


COLUMN SUPPORTING ROLLER 


Fig. shows fixed-base column supporting roller radius where 
the length the column and decimal number, which usually 
small. The load applied vertically through the center the roller and 
restrained the far-end support conditions the structure that always 
remains the same vertical axis. The significance this that the load 
unable follow the column laterally the buckling load causes deflection. 
Contrast this with the assumptions made the derivation the equations 
leading Eq. The roller will assumed ideal and frictionless that the 
contact force between the roller and the column head will line that 
normal the plane contact, that is, will always radial the roller. 

considering the possible buckling modes, apparent that the column 
cannot buckle single half wave such the classic free-standing column. 


“Old 


such initial deflection occurred, all induced moments would sense 
return the column its straight position. The assumed fundamental buck- 
ling shape that shown Fig. The contact force between the roller and 
the column head canbe broken components—one parallel tothe column 
axis and the other perpendicular it. 

For small deflections, cos and sin approximately. Thus, the 
force parallel the column axis approximated the vertical force 
whereas the force perpendicular the axis approximately The angle 
the rotation the roller during buckling and equal the slope the 
elastic curve the column head, that is, 


Coordinate axes are taken shown Fig. The two components force 
induce moments the column follows: 


Thus, using small deflection theory, the differential equation 


the solution this differential equation can written 


= k2 


The boundary conditions are 


dy/dx=- 9, atx=L..... 


Substitution boundary conditions and and either leads the 
transcendental equation 


The roots this equation vary with the value The critical buckling load 
given 


interesting that the buckling load for this condition varies from 3.17 3.29 
times the buckling load for the classic free-standing column (Eq. range 
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FIG, ROLLER AND LINK-SUPPORTING 
COLUMNS WITH CLASSIC CASES, 


INITIAL ECCENTRIC LOADING ROLLER 


Because expansion contraction the superstructure, the roller 
may move laterally with respect the original position amount and 


(b) 


FIG, 5.—COLUMN SUPPORTING ECCENTRIC ROLLER. 


thereafter remain along this new vertical axis. The conditions then are those 
shown Fig. The moment equation can written 


(15) 
The resulting differential equation 


The solution this equation 

(17) 
The boundary conditions are 


(a) 


FIXED BASE COLUMN 


Substitution these boundary conditions leads the following expression for 
tank 


The resultant eccentricity equal 


e'=e-y-mL@ 
The bending moment equation 


that the maximum bending moment occurs where dM/dx Thus, 


Therefore, the maximum moment occurs where dy/dx This locates the 
point maximum moment the head the column. Substitution the ex- 
pression for evaluated into Eq. 22, permits the calculation e'/e 
for various values This computation has been made and the results 
shown Fig. for several values Note the interesting response this 
type column when the eccentric load increased from zero For 
values approximately 0.33, the eccentric effect the vertical 
component force sufficient cause the column deflect toward the load 
axis somewhat and cause reduction effective eccentricity. For larger val- 
ues the horizontal component force plays more important role 
and finally causes the column forced out laterally away from the load 
axis. 

For design purposes, setof curves similar the familiar secant formula 
curves may developed. For eccentrically loaded column, the maximum 
stress occurs the extreme fiber and may computed the conventional 
short column formula, provided that the maximum eccentricity that occurs 
throughout the column length used. the present case, 


which the distance the extreme fiber and the maximum eccen- 
tricity that occurs determined Eq. located Fig. Eq. 
can rearranged the form 


where the radius gyration the column section. Let take arbi- 


written 


(23) 
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which the critical buckling factor obtained from Fig. for the partic- 


ular value desired, may computed exactly the solution Eq. 
13. Rearrange Eq. the form 


plot the design particular value and initial eccentricity 
arbitrary values are chosen and corresponding values and are 
computed. Then particular value specified, Eq. solved 
for This value (L/r)* may used compute the nominal value 
P/A that caused the stress Substitution Eq. 26. Finally, the nomi- 
nal value P/A found can divided suitable factor safety. 

sample set curves this kind for 32.0 ksi, and the 
factor safety 1.71 are shown Fig. This procedure not advocated 


the best possible way treat column stresses but does reflect the phi- 
losophy existing codes. 


COLUMN LOADED THROUGH LINK 


Although the free-standing column supporting roller radius proved 
considerably stiffer than the classic free-standing column, the substitu- 
tion link length will result catastrophic reduction load- 
carrying capacity. the case the roller, the upper end the link, 
where the load applied, assumed restrained that this point 
must remain along fixed vertical axis. The lower end the link can also 
apply both verticaland horizontal component force. However, this case, 
the contact force not required normal the head the column. The 
resulting components force and the deflected configuration are shown 
Fig. Coordinate axes are taken shown. The moment equation 


(28) 
The resulting differential equation 


The solution this differential equation 


1+m 


The boundary conditions are 


Substitution these boundary conditions Eq. leads the following tran- 
scendental equation 
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Eccentricity 
Extreme 
Fiber Bending 
Gyration 


7.—DESIGN CURVES FOR COLUMN SUPPORTING ECCENTRIC ROLLER 
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FIG, 8.—COLUMN SUPPORTING LINK WITH POSITIVE 
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Note that with the exception the term Eq. identical Eq. for 
the column fixed the base and pinned the top. For the equations 
are identical. Both Eq. and Eq. can solved graphically shown 
both equations represent the same column action since, link zero length 
were buckle completely Eq. indicates, then worst would lie flat 
the column head. Since has zero length and since the top the link 
restrained from sidewise movement the conditions the problem, the 
boundary conditions are converted immediately those the problem de- 
scribed For finite lengths the link, however, real non-trivial roots 
Eq. exist which are quite close zero. For 0.10, for example, the 


smallest non-trivial root 0.518. This fixes the smallest buckling load 


0.268 

This load only about one-ninth the critical buckling load for the classic 
free-standing column about one-thirtieth the critical buckling load for 
the buckling link length 0.10 would allow the top the link rest 
upon the topof the column the case and, thus, next highest 
mode failure. more likely that the buckling the link would force the 
column from beneath the supported structure with resulting catastrophic col- 
lapse. The danger poorly designed roller, such shown Fig. 
obvious. 

almost identical the slope the line tan near the points interest. 
very accurate approximate solution this problem possible, however. The 
horizontal component force causes the top the column deflect laterally 
the same manner simple cantilever beam loaded the 
end concentrated load, Thus, 


Since the column under the simultaneous action the vertical component 
force, approximate correctionfor beam-column action canbe made mul- 
tiplying Eq. the factor 


(37) 


The term the classic Euler buckling load for the free-standing column, 
given Eq. The horizontal force function the vertical load 
Por and the deflection. Thus, for small 

Per 

(38) 


Substitution Eq. into Eq. leads the approximate expression 
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FIG, SOLUTION TRANSCENDENTAL EQUATIONS, 
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Substitution 0.10 into Eq. gives, the critical buckling load 


0.2675 


Compare this with Eq. 35. 


NEGATIVE LINKS 


Neither the foregoing support conditions are compatible with those exist- 
ing for the columns lift-slab operation. These columns are more nearly 
subjected loading condition wherein the load applied through inverted 
links. That is, the loaded ends the links (the links being the jack rods) are 
below the tops the columns. the conditions this problem are used 
set governing differential equation with the link lengths given the 
resulting transcendental equation will the same Eq. except that the 
sign will reversed. This signifies that the inverted link acts the 
sense “negative” link. graphical solution for this equation shown 
Fig. together with the preceding cases. Again, link minus-zero length 
identical the link plus-zero length. Both cases degenerate into the 
conditions expressed Eq. 

noteworthy that the buckling load for the column, with either positive 
negative values approaches that the classic free-standing column 
approaches infinity. For positive values the buckling load 
infinity. That is, the classic buckling load approached from below. For 
negative values the buckling load decreases from 20.19 
the classic buckling load approached from above. Positive values may 
considered defining region instability wherein the potential energy 
maximum. Conversely, negative values define region instability 
wherein the potential energy minimum. link length could 
interesting phenomena could occur. Starting with the link the “normal” po- 
sition sothat considered positive, could applied until the critical 
load was approached. This would rather small load for any moderate 
value 

this critical load, buckling would occur with the column head being thrust 
sideways amount equal the length the link. the link were con- 
structed supposed previously, the column head could then snap back its 
original position, but now the link below that considered negative. 

Thus, the structure has passed from state maximum potential energy 
state minimum potential energy corresponding position where 
this new position, the load may increased further until the higher 
buckling load for negative values reached. The manner which the 
classic buckling load approached from above and below for various values 
and illustration the “snap-through” phenomenon for typical case, 
are shown Fig. 11. 

The case negative values has practical application the field 
lift-slab operations. the beginning lift-slab operation the jack rods ex- 
tend the full length the column that they function link with 


£ 4 


at 
' 
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the column has buckling strength equal that classic Euler pin-ended 
column, namely 


the slab were restrained from sidesway diagonal cables other suitable 
external bracing, the strength the column would increase the lift proceeded 
because the decrease inthe length the link. With the slab the full height 
the column, and the critical buckling load would that given Eq. 

most cases, may not practical brace the slab against sidesway. 
Consequently, the column will act somewhat different manner with the 
column itself acting partial bracing system for part the lift. This 
bootstrap-like action can explained reference Fig. 12. Consider the 
column neutral equilibrium with slightly buckled shape. The top 
the column, which supports the upper end the link, will deflected laterally 
the amount the level the slab where the lower end the link 
loaded the weight the slab, the column deflected amount Be- 
cause the slight difference the lateral deflection the two points, the 


link will inclined small angle that there exists horizontal compo- 
nent equal 


This horizontal component loads the column pair forces separated 
the slab, the vertical component causes moment equal 5), 
that the sum the couple formed the two horizontal components and the 
couple caused the vertical component identically zero. Hence, the part 


the column above the slab level loaded pin-ended column length 
and has its critical buckling load: 


That part the column below the slab level acts the classic free-standing 
column length where and has its critical buckling load: 


Thus, the column loaded this manner will buckle the load given Eqs. 
34, 44, whichever gives the smallest value. Eqs. and are plotted 
Fig. with the previous results. Eq. tangent the curve defined 
Eq. the For Eq. lies above the curve for 
Eq. and, thus, never controls the buckling load for practical values 
Between the values 0.65 and Eq. controls. Thus, the final 
phases lift, the free-standing columns act the classic manner. 

columns having buckling strength equal greater than pin-ended col- 
umns equal length. External bracing onthe provide enough 
stability assure this strength throughout the entire lift. 
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12.—THE LIFT-SLAB COLUMN. 


(a) (b) 
(d) 
(c) 


FIXED BASE COLUMN 
SUMMARY 


These few comments have shown that two basic equations for buckling, 
namely Eq. and Eq. are only special cases general column problem, 
the buckling load for whichis given the roots the transcendental equation, 
Eq. 34. 

Eq. results when takes value zero, whereas Eq. the conse- 
quence letting approachinfinity. problem involving column supporting 
roller radius shown separate and distinct structure from 
the column supporting the link length spite their nominal simi- 
larity. Characteristic curves for all the cases considered are shown 
Figs. and for comparison. 
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DYNAMIC RESPONSE ELASTO-PLASTIC 


Joseph Penzien!, ASCE 


SYNOPSIS 


This paper presents the results analytical investigation involving 
single mass system that elasto-plastic resistance deformation 
relationship and that subjected the ground motion measured the United 
States Coast and Geodetic Survey during the May 1940 Centro earthquake. 
The basic parameters that are varied this investigation are the natural pe- 
riod, ultimate strength, and damping the system. All results are presented 
the form graphs showing the maximum dynamic response the system 
during the simulated earthquake input plotted against natural various 
combinations ultimate strength and damping. 


INTRODUCTION 


The philosophy the design structures withstand strong motion earth- 
quakes has changed appreciably recent years due better understanding 
the true dynamic behavior structures during earthquake. However, 
the past there has been apparent wide gap between analytical results based 
elastic behavior and the results observation. This difference has been 
due the extreme difficulty, general, setting proper mathematical 
model used any dynamic analysis. 


Note.—Discussion open until December 1960. extend the closing date one 
month, written request must filed with the Executive Secretary, ASCE. This paper 
part the copyrighted Journal the Structural Division, Proceedings, the Ameri- 
can Society Civil Engineers, Vol. 86, No. July, 1960. 

Presented the February 1959 ASCE Convention Los Angeles, Calif. 
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Such model should simulate the true resistance-deformation characteristics 
the structure throughout the history deformation imposed the earth- 
quake and should also simulate the true damping the structure. For the con- 
ventional structure, where its stiffness, strength, and damping properties are 
greatly influenced the existence secondary structural elements such 
filler walls, partitions, stairwells, etc., and are also influenced foundation 
conditions, proper mathematical model difficult establish. However, 
considerable progress has been made engineers and researchers recent 
years this direction and, result, the gapbetweenthe theoretical and true 
behaviors structures subjected strong motion earthquakes gradually 
closing. 

While there have been many papers that have contributed better under- 
standing the response structures earthquakes, the writer would 
like mention two that are particularly pertinent. report2 Alford, 
Housner, and Martel has greatly advanced the understanding the 
dynamic response structures subjected earthquakes and has had consid- 
erable influence present day seismic codes. Blume has 
contributed great deal explaining the differences existing between the dy- 
namic behavior structures, predicted the usual methods analysis, 
the true behavior that has been observed. 

the purpose this paper present the results spectrum analysis 
similar those presented Alford, Housner, and Martel except that varying 
degrees plastic allowed inthe structural system. These re- 
sults show quantitatively the importance plastic deformation limiting the 


dynamic response structural systems. This effect was qualita- 
tively Blume. 


METHOD ANALYSIS 


Equation Motion.—The single mass system used this investigation 
represented the simple frame shown Fig. which the quantity repre- 
sents the displacement the mass relative its moving base, represents 
the total displacement the mass relative reference position, rep- 
resents the ground displacement, represents the force offered 


the frame (not including damping), and represents viscous damping coef- 
ficient. 


The equation motion mass can written follows: 


Me + CU + = (1) 
However, 


= Ug + u Ss enue eno 2 6 (2) 
Therefore, Eq. can written 
M i + Cc u Q M uy. 89249 462 .(3) 


This equation motion shows that any analysis determine the frame dis- 
placement, can carried out assuming fixed base. However, forcing 


ner, and Martel, Office Naval Research Contract ONR-244, August, 1951. 

“Structural Dynamics Earthquake-Resistant Design,” John Blume, Pro- 
ceedings, ASCE, Vol. 84, No. July, 1958. 
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function applied the mass equal the negative product the mass times 
the ground acceleration. 

the response the frame determined assuming elastic behav- 
ior, the frame resistance would simply the product the frame elastic 
spring constant times the displacement However, plastic deformation 
allowed, this term becomes considerably more complex. For this investiga- 
tion idealized resistance-deformation relationship assumed shown 
Fig. This relationship limits the total force developed the structure 
some maximum value Asthe mass oscillates back during 
earthquake, the frame may elastic element during certain inter- 
vals time. 


which (critical damping), and the elastic period vibration, 


Using relationships Eqs. and the equation motion (Eq. may writ- 


ten 


Solution this equation will yield the time history the dimensionless ratio 
The three dimensionless parameters and appearing this equation have 
been varied over considerable range this general investigation. 

Integration Equation Motion.—Due the complications imposed al- 
lowing plastic deformation, Eq. has been solved using numerical step-by- 
step integration procedure. this procedure the duration the earthquake 
divided into equal small intervals time (7). Assuming the quantities 
and Uare the beginning any interval, these same quantities are de- 
termined the end that interval using the so-called “mid-acceleration” 
method. 


Let and represent the known values and respectively, 


the beginning any interval. The mid-interval values and and 
respectively) are determined the simple relations. 
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Knowing and the mid-value obtained using Eq. Having 
the magnitudes and the end the interval are obtained using Eqs. 


Substituting the values and just obtained into Eq. the magnitude 
determined. The known end values and are now used the 
beginning values for the next time interval. Thus the calculations proceed 
through all time intervals similar manner. starting the calculations 
the first time interval the initial values and are zero. 

The code for the digital computer used this investigation was written 
that one could obtain either the maximum dynamic response during the earth- 
quake input the complete time history response. The latter, however, 
requires considerable machine time for printing results and, therefore, was 
carried out for special cases only. 

Ground Motion Used Analysis.—The ground acceleration 
nent) measured the United States Coast and Geodetic Survey (USC GS) 
during the May 18, 1940 Centro, Calif. earthquake was used throughout this 
investigation. Fig. shows this ground acceleration plotted against time. 


carrying out the numerical integration Eq. using this ground motion, 
time interval equal 0.01986 sec. was used. 


ANALYTICAL RESULTS 


The maximum dynamic response the single mass system during the peri- 
the earthquake was primary concern this investigation. Therefore, 
all results are presented plotting dimensionless displacement ratio against 
the elastic period vibration the system for various combinations damp- 
ing and yield strength. 

This dimensionless displacement ratio defined the ratio the maximum 
dynamic displacement occurring during the simulated earthquake the static 
displacement produced load the system. assumed, however, 
computing the latter displacement that the system has infinite strength, that is, 
yielding will take place. Therefore, the dimensionless strength parameter, 
effects only the maximum dynamic displacement the above ratio. 

All basic results the investigation are shown Figs. through 12. 
these figures, the dimensionless strength parameter 0.05 1.00, 
the parameter ranges from 0.2 sec 2.4 sec, and the damping 
ratio ranges from 15%. 

The dashed curve which repeated each these figures refers the 
purely elastic case with zero damping and therefore can used basis for 
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comparison. The ordinates this curve correspond the usual acceleration 
spectrum curve based elastic behavior. 


RESULTS COMPUTATIONS 


The maximum dynamic response for the case damping except 
plastic deformation shown Figs. and7. For periods where the 
ordinates curve No. (elastic case) strength parameter 
the response the system purely elastic throughout the earthquake input. 
However, for periods where the ordinates curve No. are greater than 
plastic deformations take place. Therefore, the solid curves these figures 
depart from curve No. points where their ordinates are equal 

Curve No. for 1.00, Fig. shows that plastic deformations take 
place for periods under approximately 0.66 sec. The total dynamic response 
shown this curve considerably less than that shown for the purely elastic 
case. the strength decreased the maximum dynamic response also 
decreases until reaches some optimum value which point further de- 
crease results increase maximum response. Comparison the 
curves Figs. through shows this optimum value approximately 
equal 0.3, 0.2, 0.1, and 0.05 periods 0.3, 0.6, 0.9, and 1.2, respectively. 
For larger periods than 1.2, the optimum value generally about 0.10 but 
curves for additional values would probably show somewhat less 
than 0.10. 

After observation these results and the ground motion shown Fig. 
would seem that the optimum for any period approximately equal the 
amplitude (expressed fraction gravity) the component ground ac- 
celeration having this same period. Therefore, for increasing values 
greater than this optimum value, the maximum dynamic response due more 
and more resonant condition that builds over appreciable increment 
time; whereas, for decreasing values lessthanthis optimum, the maximum 
response seems affected more and more the transient response result- 
ing from separate single ground acceleration pulses. 

should noted that the rapid increase ordinate all curves where 
plastic deformation involved the period approaches zero (curves and 
Fig. does not mean the maximum dynamic displacement increasing 
similar manner. This rapid increase due mainly fact that the static 
elastic displacement for loading approaches zero. Therefore, the period 
approaches zero, the numerator (maximum dynamic displacement) the dis- 
placement ratio approaches some finite value while the denominator (static 
elastic displacement for loading) approaches zero; consequently, the dynamic 
displacement ratio itself goes infinity. 

The dynamic response curves shown Figs. and have been obtained for 
conditions similar those previously discussed for the curves Figs. through 
except that viscous damping has been added the system. comparison 
curve Nos. and Fig. both which represent elastic behavior but 
which represent damping factors and respectively, shows the large 
effect that damping has the maximum dynamic response. The effect low- 
ering the strength parameter generally the same for (Figs. 
was for the previously discussed case involving viscous damping Figs. 


that is, the response reduces with decreasing certain point and 
then increases. 
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Figs. and show similar dynamic response curves for the case 10% 
damping and Fig. shows the response for 15% damping. While dynamic re- 
sponse shows general reduction with increase damping coefficient, this 
reduction becomes less and less for continued increases damping coefficient. 

general interpretation the results Figs. through would indicate 
that stiff structure, that periods approximately 0.8 sec for the 
Centro earthquake, must strengthened considerably over the more flexible 
structures order limit the dynamic response. This usually the case 
the stiffer structures are generally inherently stronger. The usual-type build- 
ing has filler walls, partitions, stairwells, and on, which considerably 
strengthen the structure for small However, since these ele- 
ments also stiffen the structure great extent, one increased 
stiffness interpreting any dynamic response data such those presented 
Figs. through 12. While the resistance-deformation characteristics such 
structures are not similar the idealized case used herein, nevertheless, 
general, the effect forcing the structure deform beyond its ultimate capac- 
ity results appreciable damping which limits the maximum dynamic response. 
This damping would produced partly velocity dependent conditions and 
partly plastic deformations the structural steel elements. damping 
factor (A) 15% would, perhaps, not unreasonable for such cases involving 
large displacements. Therefore, Fig. might used estimate this re- 
sponse. 

For the “modern-type” structure, which the structural frame not as- 
sisted the secondary elements previously mentioned, damping factor 
would more reasonable. Fig. might therefore used this case. The 
reader should keep mind that the dynamic response data Figs. through 
are maximum values during the period the earthquake and not the final values. 
Therefore, the permanent plastic deformations remaining the end the earth- 
quake period are general considerably less than those indicated Figs. 
through 12. These values may considered upper limits only. 

general the damping effect produced plastic deformation appreciable 
under conditions where the maximum dynamic response results primarily from 
resonant condition. show the damping effect plastic deformation for 
true resonant condition, that is, where the natural frequency the system 
equal the frequency the ground forcing function, sinusoidal ground 
acceleration amplitude was imposed the single mass system. 

Fig. shows the resulting steady state dynamic response plotted against 
the strength parameter this strength parameter approaches 
zero, the steady state dynamic response ratio approaches unity, which means 
that the amplitude the dynamic displacement approaches the amplitude 
the ground displacement the strength parameter approaches in- 
finity, the response becomes purely elastic resonant condition which results 
infinite amplitudes displacement. The equivalent damping factor 
plotted Fig. difined the viscous damping coefficient required 
limit the amplitude displacement for purely elastic system those shown 
this same figure for the elasto-plastic case. This damping factor varies 
from 50% for the case zero strength, for the case infinite strength. 


CONCLUSIONS 


Based the results this investigation using the Centro earthquake 
May 1940, the following general conclusions are deduced. 


FRAME RESPONSE 


(1) Buildings designed under present-day code requirements depend great 
deal damping limit their dynamic response during the period strong 
motion earthquake. Because much this damping provided plastic de- 
formations, important that structures designed withstand large dis- 
placements without ultimate collapse. 

(2) The maximum dynamic-elasto-plastic response structure decreases 
with decrease lateral yield strength until this strength reaches some opti- 
mum value which point further decrease strength increase 
dynamic response. 

(3) The optimum lateral strength limiting the maximum dynamic response 
the order 0.1 gravity less for tall flexible buildings that have long 
periods vibration and the order 0.2 gravity for low stiff buildings that 
have short periods vibration. 

(4) The damping obtained from inelastic deformations, along with the in- 
creased strength and stiffness provided secondary components suchas walls, 
partitions, stairwells, and on, are primarily responsible for the apparent 
ability many existing structures withstand strong-motion earthquakes. 


STATIC ELASTIC DISPLACEMENT FOR LOADING = Ma, 


EQUIVALENT DAMPING FACTOR 


FIG, STATE DYNAMIC RESPONSE FOR 
GROUND ACCELERATION AMPLITUDE 


While conclusions and 4are perfectly general, conclusions and based 
the May 1940 Centro earthquake and, therefore, may not necessarily show 
this same behavior for another earthquake somewhat different characteristics. 


ACKNOWLEDGMENTS 


The author wishes express his the University California 
(Berkeley, Calif.) Computer Center for the use their IBM 701 computer 
carrying out the computations for this investigation. also wishes thank 


6 0.6 
j Ug . Ugo sin wt 
0.5 
SZ) 
2 3 4 3° 


July, 1960 


Leonard Housner for his help writing the computer code and Clough 
for his suggestions and comments the preparation this paper. 


ADDITIONAL 


“United States Earthquakes,” Coast and Geodetic Surbey, Govt. Print- 
ing Office, Washington, 


“Computations onthe Response One-Mass Undamped Structures the 


Centro Acceleragraph Record,” Robison, Watson Scientific Com- 
puting Lab., Columbia Univ., New York, 


“Methods Analysis for Structures Subjected Dynamic Loading,” 


Force, Washington, C., December, 1950. 


July, 1960 


Journal the 
STRUCTURAL DIVISION 


Proceedings the American Society Civil Engineers 


MATRIX ANALYSIS PLANE RIGID FRAMES 


Fernando Venancio Filho! 


SYNOPSIS 


The problem plane rigid frames formulated matrix form, using the 
principle virtual work with displacements unknowns. The meth- 
especially suitable for highspeed digital computers. 


INTRODUCTION 


The matrix form the equations plane redundant structures the dis- 
placement and force methods derived following ideas developed 
and making adaptation from the equations deduced that author 
for plane structures composed rigid joints and straight bars with constant 
transverse section. The extension the method structures bars 
and/or variable transverse section obvious and immediate. the matrix 
formulation, one finds the joints equilibrium equations the displacement 
method the elastic compatibility equations the force method. will 
seen, however, that for the structures analyzed, terms the matrix formu- 
lation, the displacement method has great advantage over the force method. 
The equations are solved matrix methods, the final results being the form 
matrix the internal forces the ends the bars. 


Note.—Discussion open until December 1960. extend the closing date one 
month, written request must filed with the Executive Secretary, ASCE. This paper 
part the copyrighted Journal the Structural Division, Proceedings, the Ameri- 
can Society Civil Engineers, Vol. 86, No. July, 1960. 
Asst. Prof., Structures Dept., Instituto Tecnologico Aeronautica Sao Jose dos 
Campos, Sao Paulo, Brasil. 
“Energy Theorems and Structural Analysis,” Argyris, part General 
Theory, Aircraft Engineering, Vol. 27, Nos. and February and March, 1955. 


2547 ST7 


July, 1960 
DISPLACEMENT METHOD 


this method the linear displacements and and the angular displacement 
each joint are chosen redundants. 

Consider any bar the structure, bounded the joints Fig. 
Let the inclination from the bar axis relation with the axis and the 
system orthogonal axis and the axis lying along the bar axes ik. 


Express the displacements the joints and the system terms 
the ones the system xy: 


COS Sin 
COS Sin 


-sin COS 


Assume the whole structure composed n/6 bars, then 


[4] [x] 
which [A] column matrix whose elements are the displacements 
each axis system; [X] column matrix the elements which are 
the unknown displacements each joint; consequently, [a] will 
matrix, which, multiplied [X], will result [A]. The matrix [a] has as- 
pect similar the first matrix the right hand side Eq.1. depends only 
upon the geometry the structure and can found simple inspection. 
Geometrically, its elements are the joint displacements referred the bar 
axes, due unit joint displacements referred the axes. The internal 


Fig. and the external loads. 

fact, and are, respectively, the span, the area and the 
moment inertia the transverse section, and Young’s Modulus the 
material, the parts from these internal forces due the joint displacements 
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will order: 


lik 

+ 


(ik) (ik) (ik) (ik) (ik) 


write under matrix form: 


and 121 
lik 
(ik) (ik) 
\ 
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FIG, 


FIG, 


FIG, 


x (ik) 


x 
| 
(ik) 
(ik) 
(ik) 
(ik) 
= k k 
k 
\ | 
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bar pinned end the equation analogous Eq. is: 


matrix, 


which the elements the principal diagonal are the stiffness matrices 
each bar and the elements outside the principal zero matrices, 
the stiffness matrix the complete structure. 


Referring Eqs. and Eq. can written, for the whole structure. 


which [S] column matrix representing the total internal forces the 

ends the bars. The latter can obtained placing successively, one be- 

low the other, the matrices, which that the left hand side Eq. 

typical. The term [s] also column matrix obtained like manner from 

the column matrix the extreme right hand side Eq. and represents the 
internal forces the ends the bars due the external loads only. 

The joint equilibrium equations will obtained the use the virtual 

work principle, matrix form. Referring force system 


F (ik) (ik) ik 
=E 0 lik 0 3 liz i fy: (3a) 
lik 
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represented the column matrix: 


and the displacements its points application its directions, 
represented the column matrix: 


the work from the forces corresponding the displacements can ex- 
pressed matrix form as, 


which the matrix with (*) the transpose matrix. fact: 


The equilibrium conditions from the joints can expressed then, the 
vanishing the work the forces matrix [S] from Eq. through the dis- 
placements from matrix [a], Eq. thus, 


or, substituting from Eq. the value [S]: 
The matrix [X] the unknowns will found multiplying both members 


Substituting the value [X] from Eq. 11, the problem finding the 
internal forces the ends the bars completely solved: 


which [E] the identity matrix. 
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From the perusal Eq. 12, the following conclusions can reached: (1) 
Matrix [S], giving the internal forces the ends the bars, the solution 
the problem and obtained matrix operations matrices [r], [a] and [s], 
indicated Eq. 12; (2) Matrices [a] and [r] depend the geometric and 
elastic properties the structure, and matrix [s] depends the external 
loads; (3) Once the matrix, which multiplied [s] Eq. 12, ob- 
tained, several cases loading can solved, simply introducing the proper 
matrix; (4) The order the matrix [a]*[r] [a], which the matrix 
inverted, equal the number unknown displacements; (5) This form 
attack the problem particularly suitable for the use electronic digital 
computers. With matrix interpretive program, the matrix operations indi- 
cated Eq. are conveniently performed. 

Consider, simple application example, the structure Fig. Let all 
the bars have the same transverse section moment inertia and the ex- 
ternal loads indicated. The bar extensions will neglected that 
not necessary consider the area the transverse section. The unknowns, 
this problem, will the rotations and joints and and the hori- 
zontal sidesway 

Equation will be, for this particular case: 


Equation will be: 


> (12)) 6 6 
0 


(34) 


(34) 


(34) 
Fy, 


‘34 34 134 * | 
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The equilibrium condition, Eq. 10, will be: 


01000010 


|-& 


01000010 


After making the indicated operations there results: 


and 


12 Wo 
0 0 0 0 0 010 Xp 
34 34 34 
2 
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two joint and one sidesway equilibrium equations which could obtained 
the usual way. 


FORCE METHOD 


The matrix formulation the force methodcan made similar man- 
ner tothe displacement method for structures convenient take 
redundants the bending moments the ends the bars, the continu- 
ous beam problem. The force method, nevertheless, shows its greatest po- 
tentiality when the called group loadings are used and, this case, the re- 
dundants are not, generally, bending moments the ends the bars. 


FIG, 


Consider the continuous beam from with constant transverse section 
area moment inertia and with any vertical external load. 
The basic statically determinate system will obtained placing hinges 


supports and and the redundants will the bending moments the beam 
these supports. Thus, 


\ 
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or, 


which [M] the column matrix six elements representing the bending 
moments the ends the spans; [X] the column matrix formed the two 
redundants and [b] matrix which when multiplied [X] results 
shall noted that the elements [b] are bending moments the ends 
the spans due unit values the redundants. 


For any bar, 12, for instance, the slopes the ends the spans will 


given by: 
The matrix: 


the flexibility matrix the bar 12. The matrix: 


which the elements the principal diagonal are the flexibility matrices 


each bar and the other elements are zero matrices, the flexibility matrix 
the whole structure. Therefore, following Eqs. and 19, 


(12) 


(12) 
(12) 
(23) 
(34) 


which the second the right hand side Eq. 22, matrix formed 


the slopes the ends the spans the basic system, due the external 
loads, Fig. 


Concisely, Eq. can written 
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The elastic compatibility equations the supports can formulated 
the condition the vanishing the work the bending moments matrix 
[b], Eq. 18, through the slopes from Eq. 23: 


The final solution will found substituting Eq. Eq. 18: 


Eq. 24, when written out full, gives the three moment equations the 
given continuous beam. fact, 


011000 
000110 


(34) 


After making the operations, the following moment equations are found: 
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NUMERICAL EXAMPLE 


Consider the Fig. The influence axial forces neg- 
lected, habitually. For computational convenience assume 0.1. 
Matrices [a], [r], and [s] are 


1.5 -0.3 


100 


1.5 


-0.3 


-0.3 0.04 


0.12 0.006 


Matrix [a]*[r] [a] 


0.30 


9.4 0.12 


0.30 


0.12 


0.046 
0.20491 


0.02746 
1.26481 


0.02746 
0.11373 
0.11754 


1.26481 
0.11754 
30.29385 


“Digital Computation for Stiffness Matrix Analysis,” J.S. Archer, Proceedings, 
ASCE, Vol. 84, No. October, 1958. 


and 
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0.07208 7.19105 
0.23529 
0.29998 2.99978 0.29998 
0.21768 
0.01088 0.16766 0.01088 
0.08238 0.35262 


July, 1960 
Finally matrix [E] [r] [a]* 


0.01294 0.01294 
0.11762 0.11762 
0.60000 0.40000 
0.01294 
0.34119 0.34119 


7.19105 
5.29384 
0.83231 
5.29432 
2.99878 
3.35317 
0.83234 
0.35262 


0.07645 
0.11764 
0.01294 
0.11762 
0.40000 
0.25885 
0.01294 
0.65881 


Multiplying this matrix the external load matrix [s], matrix [S] will be: 


343.17 
152.93 
8.07 


152.97 


294.98 
338.56 
8.07 


43.50 
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METHOD FOR ANALYSIS MULTIBEAM BRIDGES 


John Duberg,! Narbey Khachaturian,2 ASCE, 
and Raul A.M. ASCE 


SYNOPSIS 


general method for the analysis multibeam bridges presented. Ex- 
pressions for moment, shear, and torque each beam element, due con- 
centrated load acting any point the multibeam bridge, are derived. The 
method developed applied calculating the vertical moments beam ele- 
ments solid square section constituting multibeam bridge. Influence lines 
for vertical moment multibeam bridges having two, three, and ten beam 
elements solid square section are developed. 


METHOD FOR ANALYSIS MULTIBEAM BRIDGES 


Notation.—The letter symbols adopted for use this paper are defined 


where they first appear the text, and are arranged alphabetically for con- 
venience reference the Appendix. 


Note.—Discussion open until December 1960. extend the closing date one 
month, written request must filed with the Executive Secretary, ASCE. This paper 
part the copyrighted Journal the Structural Division, Proceedings the Ameri- 
can Society Civil Engineers, Vol. 86, No. July, 1960. 
Formerly, Prof. Civ. Engrg., Univ. now with Langley Research Cen- 
Langley Field, Natl. Aeronautics and Space Admin., 
Research Asst., Civ. Engrg., Univ. Urbana, 
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Assumptions.—In the development theory for analysis multibeam 
bridges presented herein, the structure idealized with the following as- 
sumptions. 

assumed that the bridge consists beam elements placed side side 
and connected each other along the span hinges the corners the 
cross section the level top fiber. For convenient reference, the hinge 


defining the point connection any two adjacent beam elements will 
called joint. 


(a) Perspective View the Idealized Structure 


Cross-Section the Idealized Structure 


hinge 


(c) Deflected Shape the Idealized Structure 
Due Downward Loads 


FIG, 1.—IDEALIZED STRUCTURE 


Evidently, this assumption implies that, although the joints move, all the 
components movement and the rotation each element joint are equal 
those the adjacent element the same joint. Furthermore, according 
this assumption, the bridge subjected downward loads, the beam ele- 
ments deflect downward and rotate about the hinges causing separation 
elements all points except the hinge. Fig. 1(a) shows perspective view 
the idealized structure inwhich the elements are shown have solid square 
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sections. The cross-sectional view shown Fig. 1(b) while the lateral de- 
flected shape the cross section shown Fig. 
addition, the following assumptions are made this study. 


The beam elements are made homogeneous, isotropic, and elastic 
material. 


The shear strains are small and can neglected. 
The cross section each element rigid, that is, the relative defor- 


mation the parts forming the beam element small comparison with the 
deformation the beam element whole. 


All beam elements have the same cross-sectional properties and are 
the same material. 
Each beam element simply-supported over the abutments piers. 
The beam elements are prismatic. 
The external forces are concentrated loads and act vertically down. 


Centerline 
support 


Vv 


Centerline 


FIG, 2.—COORDINATE AXES 


The Coordinate Axes.—The coordinate axes are taken shown Fig. 
The x-axis the directionof the span passing through the center gravity 
any beam element. The y-axis the horizontal axis which the trans- 
verse direction the span. The z-axis the vertical axis which, this 
study, parallel the external loads. The origin the coordinate axes 
taken the center end support. 
The positive direction z-axis downward—in the direction acting 
loads. The positive direction and axes are shown Fig. 
Deformations Each Beam Element.—Let designate intermediate 
beam element between elements j-1 and For convenience, the joints 
the positive direction y-axis element will also designated 
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The beam element subjected the external load well 
statically indeterminate distributed forces which act the element Joints 
jandj-1. The lateral position the given force defined e;, which 
the lateral distance the point action from the centerline the 
beam element. Figure 3(a) shows small length element subjected 
the concentrated load which acts distance from the center ele- 
ments. The internal shears are also indicated the same figure. should 
emphasized that the only load acting the element For conveni- 
ence, the internal moments are shown Fig. 3(b). 

The statically indeterminate forces acting beam element joint are 
Qj, rj, and and joint j-1 are rj-1, The distributed forces 
and which are parallel x-axis, are opposite sign and act ele- 
ment joints and j-1, respectively. The distributed forces and are 
Similarly, the distributed forces and act parallel axis, are oppo- 
site sign, and act element joints andj-1. Fig. 4(a) shows small 
length beam element with the statically indeterminate distributed 
forces acting The internal shears are also indicated the figure. The 
internal moments are shown Fig. 4(b). 

According usual engineering theory bending and restrained torsion, 


the differential equations for the deflection element can expressed 
follows: 


which Young’s modulus elasticity psi; Ixy represents the moment 
inertia the cross section beam element about axis in.4; the 
total vertical deflection beam element in.; represents the distributed 


load acting parallel axis lb. per in.; the moment inertia the 
cross section beam element about the axis parallel axis, in.4; 
denotes the total horizontal deflection beam element in.; the tor- 
sion-bending constant, represents the total rotation beam ele- 
ment radians; the shear modulus elasticity psi; the torsion 
constant in.4; denotes the acting torque in.-lb per in.; the gross 
cross-sectional area each beam element in.2; and represents the axial 
deformation centroid beam element in. 

should pointed out that Eqs. through are applicable beam ele- 
ments any cross-sectional shape. beam elements solid hollow 


square section, small zero that the first term the left side 
Eq. can neglected. 


d4w; 
d4y; 
and 
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STATICALLY DETERMINATE FORCES 


STATICALLY DETERMINATE MOMENTS 


FIG, DETERMINATE FORCES AND MOMENTS 
ACTING SMALL ELEMENT 


e 
x 
xy Pad | 
' 
y 
z 
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STATICALLY INDETERMINATE FORCES 


FIG, 4.—STATICALLY INDETERMINATE FORCES AND MOMENTS 
ACTING SMALL ELEMENT, 
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The quantities wj, vj, and can expressed the following form: 


which the deflection beam element parallel axis due tothe 
external 


istributed load in.; the deflection beam element 
parallel axis due the statically indeterminate distributed forces 
parallel axis, due external distributed load in. (Since the only ex- 
ternal loads considered this study are those parallel axis, this quantity 
due the statically indeterminate forces rj-1, and in.; 
the rotation beam element due the external distributed torque 
radians; the rotation beam element due the statically in- 
determinate distributed forces 4j, Tj, and inin., (uj)o 
the axial deformation beam element due external concentric axial 
load (Since the only external loads considered this study are those 
parallel axis, this quantity equal zero); and the axial de- 
formation beam element due statically indeterminate distributed forces 

The quantities vj, and are the total deflections, rotation and de- 
formation beam and have been defined previously. 

For convenient reference the quantities and will 
called the statically determinate deflections, while quantities 
and will referred the statically indeterminate deflections. 

the following paragraphs, the statically determinate and indeterminate 
deflections will discussed separately. 

Statically Determinate Deflections Beam Element j.—As mentioned pre- 
viously, only vertical concentrated loads are taken into account. The only ex- 
ternal load acting beam element the concentrated load which act- 
ing any point along the span. The lateral position defined 
which the distance between the point action and the center the 


element. represented Fourier Series, the differential equations 


for the statically determinate deflections element will 


Wj = (Wj)o + (Wj)i (5) 
. 
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which the length simply-supported span in.; represents the 
distance from the point action the origin in.; and denotes 
number. 

All other terms have been defined previously. 

Integration Eqs. through will result the following expressions: 


and 


The preceding expressions for and satisfy the 


The deflections element the and direction are equal zero 
the supports, that is, for andx 


Rotation due torque zero 


Stresses due bending and due torsion bending are equal zero 
the supports, that is, for and 
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Stresses due axial forces are zero the support, that is, for 
and 


Statically Indeterminate Deflections Beam Element j.—As shown Fig. 
beam element subjected statically indeterminate distributed forces 
which are exerted the adjacent elements. indicated before, the forces 
and act parallel x-axis joints and j-1, respectively. The forces 
and act parallel y-axis, the former joint and the latter joint 
j-1. The forces and joint and j-1, are parallel z-axis. 


shown Fig. 4(a), and have the same direction the positive 
whi 


direction the coordinate axes, the direction and 
oppose the positive direction the coordinate axes. 


Referring Fig. can shown that the following relations are cor- 
rect: 


N 
~"" 


which the distance the joint from the shear center the cross sec- 


tion beam element, and the width the beam element in. 
All other terms have been defined previously. 


The distributed forces and can expressed form series 
follows: 


Similarly 


and 
and 
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the nth Fourier component. 


Substituting the preceding expressions Eqs. through the following 
will obtained: 


-EC 


Integration the preceding equations will result the following expres- 
sions: 


Ixy 
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4- 
1 


Evidently Eqs. also satisfy the boundary conditions mentioned 
the preceding section. 


result the total deflections beam element follows: 


L3 ; 


sin 
2- 


uj = (%j)o + (45)i = 
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Eqs. are for the total deflection beam element Evidently, 
similar expressions can derived for elements j-1 and 

Continuity.—In order satisfy continuity joint the following require- 
ments must met. 


Vertical deflection element joint equal the vertical deflec- 
tion element j+1 joint 

The horizontal deflection element joint equal the horizontal 

The axial deformation element joint equal the axial defor- 


The preceding requirements for continuity can expressed mathematical- 
follows: 


which the axial deformation per unit rotation element joint due 
torsion warping in. per radian per in. 

The quantities Vj+1, and are vertical deflection, horizontal 
deflection, rotation and axial deformation beam element respectively, 
the subscript indicates. 

The coefficient which appears the last term each side Eq. 45, 
generally small for closed box sections. can shown that for such box sec- 
tions the terms containing have the same order magnitude the shear- 
ing deformations associated with bending which have been neglected this 
study. These terms, therefore, can neglected when the beam elements are 
closed box sections. 

beam elements open section, however, better approximation 
include the terms containing even though the shearing deformations asso- 
ciated with bending are ignored. beam elements channel section, the 
terms containing should taken into account. 

Substituting Eqs. through for total deflections element and sim- 
ilar expressions for the total deflections element j+1 Eqs. through 45, 
the continuity relations for each will become 


sin 
Ixy 
Ixy 
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Eqs. 46, 47, and are three simultaneous equations and 
derived the basis the continuity requirements joint ide 
tions similar Eqs. 46, 47, and can derived for each joint the multi- 
beam structure relating the three forces the joint. Each joint contains 
three unknowns and Sp. 

the multibeam structure consists only two beam elements, there will 
only one joint with the three unknowns and this case, how- 
ever, all the terms containing unknown forces joints than will van- 
ish. Generally, the bridge contains joints, there will unknown 
forces each joint, and equations available for their solution. 

Once the unknown forces each joint are determined, the analysis 
multibeam structure becomes possible. Shears and moments well 
stresses each element, due concentrated loads acting various points 
the bridge, can calculated. 

Derivation Expressions for Moments, Torque Element j.— 
order determine the normal and shearing stresses each element, 


necessary calculate the moments, shears, and torque which are carried 
each element. 
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The following expressions are consistent with the ordinary theory bend- 
ing and torsion: 


(Mj)xy 


which denotes the bending moment element about the neutral 
about the neutral plane’zx the shear element corres- 

entiating, for each the following will obtained: 


E Cc n2 m2 + G J L2 [7(j-1)n +1 Tin | 


and 
(Vj)zx = - 5 [9G-1)n + Tin} 
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and are functions distributed indeterminate forces 
order calculate the preceding quantities, necessary calculate 
these indeterminate distributed forces from 47, and 48, the continuity 
equations for joint and similar continuity equations for the other joints. 

The quantity presented Eq. seems the most significant 
quantity this particular study, for through the study this quantity 
that the lateral distribution load multibeam bridge can defined. 

The quantity presented Eq. generally small. However, 
for accurate ation normal stress beam element this quantity 
may necessary. The quantities and T;, although not used 
this study, may used the calculation shear stresses beam ele- 
ment. 

Determination Normal Stresses Each Beam Element.—The normal 


stresses beam element caused and (Mj)zx can expressed 
follows: 


from the xy-plane psi; the distance the point under consideration 
from the xy-plane in.; denotes the normal stress element due 

atadistance from the zx-plane psi; and the distance the 
oint under consideration from the zx-plane in. 


and 
The uniformly distributed normal stress caused axial load 


an 
du; 
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Substituting the differential Eq. terms Eq. 63, the value 
for each will 


addition the normal stresses caused Mj) and axial load, 
element subjected tonormal stresses due torsional warping. 


The normal stresses due this effect can presented the following ex- 
pression: 


which the normal stress element due restrained warping 
caused torsion psi; and the axial deformation the point under con- 
sideration per unit rotation due warping in. per radians per in. 

Substituting the second differential Eq. terms Eq. and 
rearranging the value for each will 


LEQ 


The total normal stress element can obtained the addition 
Eqs. 61, 62, 64, and follows: 


(67) 


dx2 
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Eq. presents accurately the normal stress any point beam element 
From practical point view, the important part the total normal 
and are The quantity may important for 
channel sections. 

order solve Eq. necessary determine the distributed inde- 
46, 47, and for the continuity joint and similar equations for the con- 
tinuity other joints provide sufficient equations for the solution these 
distributed indeterminate forces. 

Shearing Stresses Each Beam Element.—Shearing stresses beam ele- 
ment may calculated from the shear forces and torque given Eqs. 56, 
57, and 58. 

The shearing stresses caused each and values 
are proportional these values, and the constant proportionality can 
derived from the geometry the particular section. 

Since, this study, use has been made shearing stresses, and 
inconvenient present them that they may applicable any shape 
section, their derivation has been omitted. 

Convergence.—The equations derived this section are developed for any 
value order get fairly accurate résults, solutions for three four 
values will necessary. The total results will the summation re- 
sults for each 

should pointed out that Eqs. 46, 47, and for the continuity joint 
and similar equations for other joints not yield convergent results for suc- 
cessive values Consequently, actual solution for the distributed inde- 
terminate forces not possible Eqs. 46, 47, and 48. 

However, from analytical point view, the moments, shears, torque, 
and stresses represented Eqs. through are the quantities generally 
needed. These equations give convergent results for successive values 

this particular study, the value the quantity used for the 
determination the lateral distribution load. For the determination this 


quantity, only the summation three terms corresponding three values 


MULTIBEAM BRIDGES WITH BEAM ELEMENTS HAVING 
SOLID SQUARE CROSS SECTION 


General Expression for Continuity and order show the appli- 
cation the method developed the preceding section calculating 
each element, the simplest case will considered. will assumed that 


the multibeam bridge consists beam elements having solid square cross 
section. 


solid square sections the following expressions are correct: Ixy Izy 


= 12’ 
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Substituting the preceding expressions Eqs. 46, 47, and 48, the continu- 
ity equations joint for each become: 


The expression for for beam elements with solid square section 
from Eq. becomes 


Two-Element type structure consists only two solid 
square section elements connected means hinge. The joint will 
designated the left element will referred and the right element 


= 
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can written follows: 


2(0.296 


simultaneous solution Eqs. 72, 73, and yields the following expres- 
sions for and respectively: 


1 0.857 D 1+ b 


sin 


Since two-element structures are academic interest, and they are dis- 
cussed merely show the application the theory solution actual 

can shown that, this case, L/b has comparatively small effect 
L/b the influence lines (M; drawn for unit load that 
moves transversely from left right the midspan section. 

order obtain the ordinates the influence line for xy, the values 
and Sjn from Eqs. and are substituted assuming 


2 e; 2 e; 
(73) 
Pp; - Pp; 
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The ordinates the influence line have been tabulated coefficients 
for each Table Three values have been used, namely 
Evidently, the even values result zerovalues for and 
Sjn, Since this study assumed that L/2. Summation moments for 
each value results the total moment. The total moments for various load- 
ings are presented the last column from the left Table More accurate 
values (Mj) may obtained taking additional terms the series 


TABLE 


ORDINATES THE INFLUENCE LINE MIDSPAN 


FOR TWO ELEMENT STRUCTURES WHICH L/b 


)/L carried element 
Position Solution 


Point 


Point 


the total Table1, the accuracy may somewhat improved. 
can shown that error the quantities shown Table the neigh- 
borhood 7%. The influence line for (Mj)xy shown Fig. 
Three-Element and Ten-Element 68, and can 
rewritten form allowing their immediate use for the solution any num- 
ber elements simply isolating the coefficients the unknowns. The 
m-element-structure leads the solution the linear matrix equation: 


ST7 
7 
AY 
n=1 0.1027 0.0999 
Point n=3 0.0126 0.0100 
n=5 0.0052 0.0028 
0.1205 0.1127 
n=1 0.1020 0.1006 
n=3 0.0119 0.0107 
n=5 0.0046 0.0034 
0.1185 0.1147 
n=1 0.1013 0.1013 
0.0113 0.0113 
n=5 0.0040 0.0040 
0.1166 0.1166 
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[A] 
[B] 


The matrices [A], and are the stiffness and load matrices, respective- 
ly, for The [x], values obtained are used the evaluation the mo- 
ment carried each element given Eq. 54. 
The matrix [x], has been evaluated the the electronic digital 
computer the University Illinois, using standard library routine. 
The solutions obtained for values three-element structure have 
been summarized Table For convenience, the beam elements have been 
designated Il, and III beginning from the left. The ordinate influence 


matrix order 3m-3 


ELEMENT 


FIG, LINE FOR (M)xy FIG, LINES FOR (M)xy 
ELEMENT WHEN UNIT LOAD ELEMENTS IAND 

MOVES TRANSVERSELY UNIT LOAD MOVES TRANS- 
MIDSPAN, TWO-ELEMENT VERSELY MIDSPAN, THREE- 

STRUCTURE ELEMENT STRUCTURES 


lines are listed terms for each value corresponding particular 
position load. before, the solution based three values namely, 
The sum components moment for each also shown 
Table 

The solution for the three-element structure based upon L/b 20. This 
value L/b average value and has been chosen, since L/b ratio has 
comparatively small effect upon the moments. 


| | 
ELEMENT 
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FIG, 7.—INFLUENCE LINE (M)xy EACH ELEMENT WHEN 
UNIT LOAD MOVES TRANSVERSELY MIDSPAN. 


0.1L 
ELEMENT 
ELEMENT 
0.1L 
0.05L 
ELEMENT 
ELEMENT 
0.1L 
0.05L 
ELEMENT 
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Fig. shows the influence line for Elements The 
ordinates are given terms 


The solutions obtained for (M)xy values ten-element structure are 
listed Tables and for L/b 10, L/b 50, respectively. For convenience 
the elements are designated through The solutions are based three 


TABLE 


THE INFLUENCE LINE MIDSPAN 


FOR THREE ELEMENT STRUCTURE WHICH L/b 


“A 


Solution 


values namely and The results given Tables and are simi- 
lar those Table 
Fig. shows the influence lines elements through for two 


values L/b, and 50. For convenience this case, the ordinates are giv- 


0.7074 0.6690 0.6500 
Point n=3 0.1016 0.0690 0.0546 
0.0460 0.0218 0.0123 
0.8560 0.7598 0.7169 
n=1 0.6978 0.6548 
Point n=3 0.093 0.0736 0.0581 
0.0407 0.0257 0.0147 
0.8320 0.7731 0.7276 
n=1 0.6882 0.6788 0.6595 
Point n=3 0.0852 0.0781 0.0618 
n=5 0.0543 0.0296 0.0171 
0.6738 0.6788 0.6738 
Point 0.0741 0.0781 0.0741 
0.0251 0.0297 0.0251 
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From Tables and can shown that for three values (M) 
values for three-and ten-element structures are roughly within the 
static values. may said that the errors, thiscase, are the same or- 
der magnitude those two-element structures. 

Interpretation Results.—A study influence lines for (M)xy two- 
element and three-element multibeam structures, shown Figs. and 
indicates the load distributed almost uniformly among the elements. 
should pointed out that the influence lines these figures are drawn for 
L/b 20. For L/b values greater than 20, the influence lines may tend be- 
come horizontal. For L/b values less than possible that the loaded 
element may carry somewhat more load than the others. 

From the shape the influence lines and the pattern lateral load dis- 
tribution, one may conclude that two-element and three-element structures 
the lateral distribution uniform. this case, one element 
loaded all elements the bridge deflect almost the same amount. true 
that two-element and three-element structures have practical significance, 
nevertheless, their study contributes the understanding the problem. 

study Fig. indicates that the lateral distribution load ten- 
element structure depends the L/b value. prepared for two values 
L/b, and 50. From the figure may observed that for L/b 50, the 
structure similar two-element three-element structures. For 
L/b 10, however, the loaded element carries greater percentage the 
load, and the structure more rigid. 

should pointed out that the range L/b, selected this study, 
unreasonably high. From practical point view, very seldom that L/b 
exceeds about 30. may concluded, therefore, that, generally, ten- 


element structures greater percentage the load carried the loaded 
element. 


APPENDIX. NOTATION 


gross cross-sectional area beam element, in.; 

torsion-bending constant, in.8, 

distance between the point action concentrated load from 
the support the beam assumed origin, in.; 

Young’s modulus elasticity, psi; 


eccentricity, that is, the transverse distance between the point 
action external load and the center beam element, in.; 


overall depth beam element, in.; 


Ixy moment inertia the cross section beam element about 
the axis parallel axis, in.4; 


moment inertia the cross section beam element about 
the axis parallel axis, in.4; 


torsion constant the cross section beam element, 
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letter designating typical beam element and the joint immedi- 
ately the right the element; 


length simply-supported span, in.; 


total bending moment element about the neutral plane xy, 


bending moment about the neutral plane due the statically de- 
terminate external forces, 


bending moment about the neutral plane due the statically in- 


total bending moment element about the neutral plane zx, 


bending moment about the neutral plane due the statically de- 
terminate external forces, 


bending moment about the neutral plane due the statically in- 


number designating number terms the Fourier Series; 
concentrated load acting parallel axis element 
distributed load acting parallel axis element per in.; 


each the indeterminate distributed loads which are equal mag- 


nitude and opposite sign and which act parallel axis ele- 


amplitude the distributed load corresponding Fourier 
component, lbs per in.; 


each the indeterminate distributed loads which are equal mag- 


nitude and opposite sign and which act parallel axis ele- 
ments and j+1 joint lbs per in.; 


amplitude the distributed load corresponding thenth Fourier 
component, lbs per in.; 


each the indeterminate distributed loads which are equal mag- 
nitude and opposite sign and which act parallel axis ele- 
ments and j+1 joint lbs per in.; 


amplitude the distributed load the nth Fourier 
component, lbs per in.; 


total torque, 
distributed torque acting beam element per in.; 
total axis deformation element in.; 


axial deformation element due statically determinate exter- 
nal axial forces, in.; 


axial deformation element due statically indeterminate dis- 


rj = | | 
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shear corresponding 1b; 


xy) 


shear corresponding 1b; 


total deflection element in.; 


deflection element parallel axis due the statically de- 


terminate external loads, in.; 


deflection element parallel axis due the statically inde- 


total deflection element parallel axis, in.; 


deflection element parallel axis due the statically de- 


terminate external loads, in.; 


deflection element parallel axis due the statically inde- 


distance any point the cross section beam element from 
the plane, in.; 


distance any point the cross section beam element from 
the plane, in.; 


distance the joint from the shear center beam element, 
in.; 


total normal stress element psi; 


plane, psi; 


normal stress element due (My) distance from the 
plane, psi; 


normal stress caused axial load, psi; 


normal stress due torsional warping, psi; 


total rotation element radians; 


rotation element due statically determinate external loads, 
radians; 


rotation element due statically indeterminate distributed 


axial deformation per unit rotation due warping, in. per radi- 
ans per in.; and 


axial deformation per unit rotation due warping level the 
joint, in. per radians per in. 
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FOREWORD 


The papers that form the ASCE-CRC Symposium Metal Compression 
Members are representative all phases Column Research Council ac- 
tivities. Several papers represent summaries both past and current re- 
search, prepared specifically assist the preparation the “Guide De- 
sign Criteria for Metal Compression Members.” Also included are reports 
recently completed research pony truss bridges and plate girders. 

the January 1960 Journal the Structural Division Bruce Johnston 
documented the inception the CRC and the many activities undertaken 
the Council. 

Thomas Kavanagh, the February 1960 Journal the Structural Di- 
vision, reviewed the state knowledge concerning the effective length 
framed columns. 

the April 1960 Journal the Structural Division, Robert Barnoff and 
William Mooney studied the effect floor-system participation pony 
truss bridges. 


Bruce Johnston presented digest the “Guide” the April 1960 Jour- 
nal the Structural Division. 


Note.—Discussion open until December Separate discussions should sub- 
mitted for the individual papers this symposium. extend the closing date one 
month, written request must filed with the Executive Secretary, ASCE. This paper 
part the copyrighted Journal the Structural Division, Proceedings the Ameri- 
can Society Civil Engineers, Vol. 86, No. July, 1960. 

Struct. Metals Div., Fritz Engrg. Lehigh Univ., Bethlehem, Pa. 
Research Assoc., Lehigh Univ., Bethlehem, Pa. 
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summary presented theoretical and experimental investigation 
the strength centrally loaded columns influenced residual stresses 
and variations the yield-stress level. shown that the basic strength 
structural steel columns containing residual stresses may expressed 
terms the tangent modulus. Approximations suitable for design use are 
suggested. Information both the yield stress level and the magnitude and 
distribution residual stress structural members presented. 


research project the “Influence Residual Stress Column Strength 
and the Mechanical Properties Rolled Shapes” has been progress Lehigh 
University under the guidance Research Committee the Column Re- 
search Council. This Committee was assigned the task determining the re- 
lationship between material properties and the strength columns, and the 
first pronouncement the Council (based the recommendation Commit- 
tee was issued memorandum states that the critical ul- 
timate failure load column given 


which the applied average maximum stress column, denotes 
the tangent modulus, the total length pin-end column, represents 
the radius gyration, and L/r the effective slenderness ratio. 

Eq. cannot applied steel columns the stress-strain relationship 
determined from small coupon cut from the section. Early work clearly 
showed this true, and later studied completed prior the time that 
this general investigation was started was shown that residual stresses 
might account for differences column strength much 30%below that 
which would inferred from coupon tests. 

The column curve depends the stress-strain relationship. The latter, 
turn, dependent two important factors; these are: (a) the magnitude and 
distribution residual stresses, and(b) yield stress level the ma- 
terial. Therefore, the objectives the investigation were: (1) determine 
the magnitude and distribution residual stresses columns, and (2) de- 
velop methods predicting the influence these residual stresses 
strength. necessary parallel study, the program included adetermination 
the basic yield stress level the material which columns would fab- 
ricated. 

The program included tests coupons the type performed the mill, 
tests stub columns (short lengths full cross-sectional area), and column 
tests. For the same shapes, residual stresses were measured the section- 
ing technique. Theories were developed for predicting column strength, and 


“The Basic Column Formula,” Column Research Council, Tech. Memorandum No. 
May, 1952. 
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from the measurements made was possible obtain correlation with the 
theory. Maximum strength columns formulas could then written. 

the purpose this paper summarize the findings the investiga- 
tion date and discuss the significance thereof. Reference made, 
throughout, the various progress reports that contain the detailed experi- 
mental and theoretical work. This paper concerned primarily with rolled 
wide-flange shapes ASTM Designation structural steel. brief discus- 
sion built-up columns (welded and riveted) and low-allow high-strength 
steel columns included. The scope limited centrally loaded columns. 


RESIDUAL STRESSES 


Formation.—Residual stresses areformed structural member asa re- 
sult plastic deformations. rolled shapes these deformations always oc- 
cur during the process cooling from the rolling temperature air temper- 
ature; the plastic deformations result from the fact that some parts the 
shape cool much more rapidly than others causing inelastic deformations 
the slower cooling portions. (The flange tips wide-flange shape, for ex- 
ample, would cool more rapidly than the juncture flange and web.) The 
mechanism which residual stresses are formed has been described.4,5,6 
Residual stresses also are formed result fabrication operations. 
The process cold-bending that required inthe straightening operation and 
the process cambering both introduce residual stresses due plastic de- 
Residual stresses are also introduced during the welding oper- 
ation result the localized heat input and resultant plastic deformation.7 
Thus, plastic deformations are necessary order that residual stresses 


formed. hot-rolled welded members, the part cool last usually 
state tensile residual stress. 


Magnitude and Distribution.—Methods for determining cooling residual 
stresses plates and for obtaining qualitative estimate stressesin wide- 
flange shapes are The magnitude and distribution these stresses 
depend the shape cross section, initial temperature, cooling conditions, 
and material properties. 

The measurement residual stresses confirms the trends predicted the- 
oretically. considerable number such measurements have been made and 
they permit good estimate made the magnitude and distribution 
residual stresses likely encountered hot-rolled wide-flange mem- 
bers.5,8 For the purpose making these measurements shapes were selected 
with widely differing geometry. Fig. presents some the measured results, 
showing the magnitude and distribution stresses across the flange and web. 
While the variation considerable, the general pattern the flange simi- 


lar. The residual stresses were determined the “method sectioning” 
described 


“Residual Stress and the Compressive Strength Steel,” Huber and 
Beedle, Final Report Pilot Program, Welding Journal, Vol. 33, No. 12, December, 
1954, 589-s. 

“Residual Stresses Wide-Flange Beams and Columns,” Huber, Lehigh 
Univ., Fritz Lab. Report 220A.25, July, published Proceedings, ASCE). 

“Residual Stress and the Yield Strength Steel Beams,” Yang, 
Beedle, and Johnston, Welding Journal, Vol. 31, No. April, 1952, 205-s. 

Column Strength,” Fujita, Dissertation, Lehigh Univ., August, 1956. 

Magnitude and Distribution Residual Stress,” Fujita, Lehigh Univ., 
Fritz Lab. Report 220A.20, June, 1955, 
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FIG. 1.—RESIDUAL STRESS DISTRIBUTION WIDE-FLANGE SHAPES 
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The magnitude and distribution residual stresses inthe flanges also may 
estimated from the results stub-column test.4,9 

Table contains listing all shapes studied thus far. also contains 
summary all the measurements made the program. Table presents 
the average value the residual stress different positions the cross 
section and gives the variation well. Insofar columns are concerned, 
the most important the stresses are those the flange tips and from 
this program tests the average compressive stress 12.8 ksi witha 
maximum 18.7 ksi and minimum 7.7 ksi. 

The magnitude the residual stress the flange tips may also deter- 
mined indirectly from stub-column test, and the difference between the 
yield stress level and the proportional limit. This value has been found 

approximately 12.0 ksi which agrees well with the measured value 12.8 ksi 
noted 

addition average values, Table gives the maximum and 
minimum values measured residual stresses. Fig. 2(a) gives the frequency 
distribution flange-tip stress determined actual sectioning.5 Fig. 2(b) 
shows the same information determined indirectly from the proportional 
limit obtained for the same group “stub-column” tests.5 

The variation residual stress within material from one ingotis relative- 


but larger variations may exist between material from different 
lots. 


ksi 


FREQUENCY 


FREQUENCY 


Or = (Oy (Compression) 


STUB-COLUMN TESTS 


(a) FLANGE TIPS 


FIG, 2.—MEASURED RESIDUAL STRESS 


Theoretical studies show that cooling residual stresses are constant along 
the member except for distance approximately equal the larger cross- 
sectional dimension the ends. This due the uniformity the shape and 
its manufacture. Measurements are agreement with the theory shown 
Fig. 3.2 While atrend may exist for the variation residual stress func- 
tion the geometry the cross-sectional shape, has not been possible, 
yet, show precise relationship. 

Influence the Apparent Stress-Strain Relationship.—As implied previ- 
ously, residual stresses affect the average stress-strain relationship the 
complete cross section. Fig. presented show this influence diagram- 

matically. Fig. 4(a) shows short length wide-flange shape with sim- 
plified linear distribution residual stress the flanges. The cutting ofa 


«The Influence Residual Stress the Instability Axially Loaded Columns,” 
Huber, Tall, and Beedle, Lehigh Univ., Fritz Lab. Report 220A.26 (in 
preparation). 
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coupon from the flange the member would relieve the residual stresses that 
were present the coupon prior sectioning, and the stress-strain relation- 
ship determined from this coupon would shown the dashed line Fig. 
4(b). If, now, the load considered being applied the entire cross sec- 
tion containing its residual stresses, evident that when the applied stress 
becomes equal the difference between and then yielding will com- 
mence the flange tips. Thus, 


which the yield stress the proportional limit and denotes the 
yield stress level. The superposition stresses the “proportional limit” 


TABLE 2.—RESIDUAL STRESSES SHAPES DUE COOLING 


Mini- 
mum 
Columns 


Web Center, 


Note: Residual stress measured gage 
length at the center of the stub column 


Average stress at flange centers 
(Ten) 
Average stress at flange centers ksi 


20 Average stress at flange edges 


Average stress at flange edges hei 


LENGTH STUB COLUMN -FT 


LENGTH STUB COLUMN 


LENGTH STUB COLUMN 


(a) ALONG BEAM (b) WITH STUB-COLUMN LENGTH 


FIG, 3.—VARIATION RESIDUAL STRESS 


shown Fig. 4(c). Yielding occurs when the flange tip residual stress 
plus the applied stress equal 

When more load applied, the average stress and average strain are 
longer proportional one another and non-linear stress-strain relationship 
results for the section whole. Fig. 4(d) shows wide-flange shape with 
flanges partially yielded (shown shaded) and the corresponding stress distribu- 


tion for Above the proportional limit and below the yield stress 
the strain given 


Flange Edge, Flange Center, 
Maxi- Aver- Mini- 
mum age mum 
8.8 -21.8 -41.0 
Region containing cold-bending yield lines > 
(Ten) 
1°) 3 6 9 
20 
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Stress 


Average Strain 


(b) 


Portion 


Plastic Portions 


FIG, 4.—INFLUENCE RESIDUAL STRESS 
STRESS-STRAIN CURVE 
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FIG, 5.—STRESS-STRAIN CURVE 
FOR ASTM 
UMNS BASED MEA- 
SURED RESIDUAL 
STRESSES 
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which Young’s modulus elasticity and defined Fig. 4(c).4 The 


stress distribution completely yielded section shown Fig. 4(e). 
soon yielding has progressed across the entire section, there further 
influence residual stress. The average stress versus average strain curve 
for the entire cross section shown Fig. 4(b). The yield stress level 
unaffected the residual stresses. 

Thus, the effect residual stress the stress-strain relationship 
lower the proportional limit and cause the stress-strain diagram non- 
linear beyond that point and the yield stress level. The proportional 
limit may computed from Eq. For typical wide-flange column 
3.0) with -34 ksi and with -13 ksi the theoretical stress-strain 


FIG, 6.—TYPICAL STUB-COLUMN TEST 


shape and the web area. 

obtain experimental correlation witii the foregoing predictions, tests and 
measurements stub columns have made. typical stub column test 
shown Fig. Typical results are shown Fig. which also shows the con- 
trast with the curve based coupon tests. Fig. the stress-strain curve 
annealed stub column (the residual stresses have been removed); the 
comparison with Fig. clearly verifies the influence residual stresses 
the average stress-strain curve. 

interest the fact that the average stress-strain curve, Fig. near- 
identical with the average stress-strain from stubcolumn 
tests (Fig. 18, presented subsequently). 

Cold-Bending Residual Stresses.—Residual stresses due cold bending 
can predicted with reasonable assuming certain initial cooling 
stresses and the extent deformation. Some measurements are shown 
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Fig. general, the maximum and minimum stresses are the same or- 
der magnitude the cooling residual stresses. While the influence the 
stress-strain curve might predicted, this has not been done studies 


(described subsequently) show unnecessary for the purposes this 
present study. 


MECHANICAL PROPERTIES 


Proportional Limit.—The general effect residual stresses the stress- 
strain relationship was discussed under the heading “Residual Stresses: In- 
fluence the Apparent Stress-Strain Relationship.” The proportional limit 
reduced below that obtained coupon tests, and may computed either from 
Eq. may measured stub column test. Fig. clearly shows this 
effect, and further comparisons have been 


COLUMN STRESS- 
STRAIN CURVE FOR AS- 
DELIVERED MATERIAL 


FIG, COLUMN STRESS- 
STRAIN CURVE FOR AN- 
NEALED MATERIAL 


The average value the proportional limit 21.7 ksi determined indi- 
rectly residual-stress measurements this program tests. The aver- 
age value determined the offset method from stub column tests 20.7 ksi, 
and the frequency distribution curve for these specimens shown Fig. 
10. Since the offset method was used determine the proportional limit, the 
actual value even lower than 20.7 ksi. expected that the stub 
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FIG, 9.—COLD-BENDING RESIDUAL 
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column proportional limit would lower than the value determined indirectly 
from Eq. the first place, the residual stresses are probably greater than 
the measured value small amount. Secondly, the flange-tip values Orc 
are averages measurements corners and two sides, 
and deviations from these averages will reflected lowering the pro- 
portional limit. Small inaccuracies alinement also influence the observed 
proportional limit. 

The yield may bedetermined from the results ASTM accept- 
ance test may determined from laboratory coupon test from the 
results laboratory stub column test. Factors such upper yield point, 
strain rate, and the web strength compared strength cause the yield 
stress level full cross section (stub column) tested the laboratory 
compression markedly less than the tensile yield pointdetermined the 
routing ASTM acceptance test.4,10 


TEST (Web) 


STRESS 


Zero 
Strain Rate 


WEIGHTED COUPON AVERAGE 
STUB COLUMN 


FIG, 11.—INFLUENCE SEVERAL VARIABLES YIELD STRESS LEVEL 


Fig. illustrates the different stress-strain curves that may obtained 


depending the type test that performed and shows diagrammatically 
the influence the various factors follows: 


(1) Starting with curve the yield value reported mill-type accept- 
ance test usually (though not always) the upper yield point. Occasionally 
specimens will not exhibit upper yield point (dashed line, curve B), which 
case the yield stress usually recorded strain 0.5%. The upper yield 
point from 0%to than the yield stress level. 

(2) Comparing curves and Fig. 11, the ASTM acceptance test usually 
will show higher yield stress level primarily because the effect strain 
rate. The “static” yield level (the value measured zero strain rate)10,11 


“Qn the Yield Properties Structural Steel Shapes,” Tall and Ketter, 
Lehigh Univ., Fritz Lab. Report 220A.33, November, 1958. 

“Material Properties, Residual Stresses and Column Strength,” Gozum 
and Huber, Lehigh Univ., Fritz Lab. Report 220A.14, May, 1955. 
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web coupon from 10% 15% lower than the strength obtained the ac- 
ceptance 

(3) There difference strength between the various elements 
rolled shape, the web being stronger than the flange. Thus the average yield 
level stub column lowered still further comparison with the mill ac- 
ceptance test because the latter made from the web (curves and D). 

(4) Fig. seen that stub column(tested “zero strain 
rate”) reflects the effect residual stress upon the stress-strain curve and 


averages out the differences between web and flange material, shown 
curve 
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FIG, 12.—STRESS-STRAIN CURVE SHOWING 
INFLUENCE STRAIN RATE 
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FIG, 13.—INFLUENCE STRAIN RATE 
YIELD STRESS LEVEL 


Tests have shown that even “very slow” laboratory strain rate used 
testing coupons(an elastic strain rate micro-in. per in. per sec) can raise 
the apparent yield stress level much 5%.11 The effect strain rate 
the yield stress level typical tension coupon test shown Fig. 12. 
From the results such tests the relationship shown Fig. may ob- 
tained. The measured yield level compared with the static yield level shown 
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function the plastic strain rate. Since the plastic strain rate mill 
test about 1,000 micro-in. per in. per sec, increases the yield stress 
above the static value may expected high 16%. 

interest note that the yield stress level obtained from stub- 
column test agrees well with tests tension and compression coupons the 
results the latter are averaged according their respective areas the 
cross section, and both the coupons and the stub column are tested the 
zero strain rate.10,11 This shown diagrammatically Fig. (curves 
and F). 

comparison the results acceptance tests and stub-column tests 
importantin establishing yield stress use abasic column 
formula. this relationship can established for fairly large sample, then 
could applied with confidence the larger body acceptance test data 
available the mills give reasonable estimate the actual strength 


Point 


No. of Mill Tests 3124 


Minimum Specified Yield Strength 


Yield Point Low 31,090 psi 
Strength 


Standard Deviation 


Strength = 39,360 psi Coefficient 78 % 
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FIG, 14.—CUMULATIVE FREQUENCY DISTRIBUTION 
YIELD STRESS (MILL TESTS) 


structural steel columns. Figs. 14,15 and show the distribution the yield 
stress determined number methods follows:4,10,11 


ASTM acceptance tests the mill: Fig. (3,124 Specimens), Fig. 
(3,010 Specimens), and Fig. 16(a) (35 Specimens). 

Simulated ASTM tests: (Fig. 16(b) (35 Specimens). 

Stub-column tests: (Fig. 16(c) (35 Specimens). 

Stub-column tests: (Fig. 16(d) (47 Specimens). 


Items and 3(Fig. 16) are for the same control group. Comparing Fig. 16(a) 
with Fig. 16(b) shows that the mill test results may approximated within 
about laboratory test that simulates mill test procedures. The aver- 
age value the control group mill tests (Fig. 16(a)) was 42.9 ksi with 
standard deviation 4.4 ksi. Further comparisons this type are avail- 

The average value the yield stress level determined simulated mill 
tests was 41.2 ksi(Fig. 16(b)), while the average strength the corresponding 
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stub columns (Fig. 16(c)) was 34.0 ksi. The average value the yield stress 
level for all stub columns tested the program was 34.5 ksi with variations 
from 24.6 ksi 43.0 ksi (Fig. 16(d)). 

The probable ratio thebasic compressive static yield stress 
level) the “acceptance test” strength was found 0.80. Fig. plot- 
ted the distribution the individual ratios stub column yield level mill 
test “yield point.” The average 0.80 with minimum 0.62 and maxi- 
mum 0.92. From this information would concluded that the static yield 
stress level wide-flange column averages about 20% less than the value ob- 
tained the ASTM-type acceptance test. pointed out, this difference due 
mainly the strain-rate effect, but alsois the higher yield stress 
the web and the difference between upper and lower yield points. 

Applying this average ratio (0.80) the average value the mill tests 
shown Figs. and (41.7 Ksi) there obtained probable compressive 
strength for this material 33.4 ksi. The average value obtained from all 
stub columns tested the program was 34.5 ksi suggesting that the sample 


-——All specimens (3010 Specimens) 
— (1699 Specimens) 
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FIG, 15.—MILL TEST YIELD POINT FREQUENCY 
DISTRIBUTION CURVES 


was fairly representative. The corresponding average stress-strain curve for 
stub column shown Fig. 18. 


CENTRALLY LOADED COLUMNS 


Influence Residual Stress Column Strength.—For mate- 
rial exhibiting idealized stress-strain cruve and free from residual stresses, 
the buckling strength defined the Euler buckling curve and the yield 
stress the material. The existence residual stresses the cross sec- 
tion reduces the buckling strength, since there early localized yielding 
certain portions the cross section. This reduction greatest when the 
slenderness ratio between and 90. 

This effect residual stress may illustrated example. Consider 
column rectangular cross section containing residual stresses, and bend- 
ing about the weak axis. construct the column curve, first the tangent 
modulus determined for various stress levels using stress-strain dia- 
gram such that shown Fig. Then curve would drawn 
shown Fig. 19(b). Applying the tangent-modulus formula, Eq. the re- 
sulting column curve stress-slenderness ratio would obtained, Fig. 19(c). 
The influence residual stresses when buckling occurs the inelastic range 
thus seen comparing the solid line with the dashed line obtained for 
member without residual stresses. 


600 Average = 44) ksi-———> 
NUMBER OF 
SPECIMENS 
200 
30 4 50 60 


COLUMN STRENGTH 


~Average 429 ksi 
(s=44 


Oy ksi 


(a) Mill Test (Web Coupon) 


Average 41.2 ksi 
ksi) 


—— Average 34.0 ksi 
ksi) 


FREQUENCY Specimens 


(c) Stub-Column Test 
(Static Strain Rate 


Average 345 ksi 


FREQUENCY 


(d) Stub-Column Test 


16.—YIELD STRESS LEVEL DETER- 
MINED VARIOUS METHODS 


155 
FREQUENCY Specimens 
20 » 40 50 60 
FREQUENCY Specimens 
20 3% 40 50 60 
Oy ksi 
(b) Simulated Mill Test 
20 30 40 50 60 
ksi 
40 
ksi) 
20 30 40 50 60 


July, 1960 


Eq. valid only for the special case rectangle bent about the weak 
axis. Proceeding now toadiscussion the problem determining the strength 


actual column, the basic equation for the critical strength column 
containing residual stresses given by® 


which the moment inertia that portion the cross section which 
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207 ksi 


FIG, 18.—AVERAGE STRESS-STRAIN 
CURVE FOR STUB COLUMNS 


remains elastic (Fig. 4(d), for example). terms the average critical 
stress, Eq. may written 


Eq. the basic equation for column containing axially symmetric cooling 
residual stresses. Since the flanges contribute significantly the flexural re- 
sistance, evident that residual stresses the flanges are more 
pronounced influence column strengththan are residual stresses inthe web. 
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Effect Flexure Axis.—There pronounced difference the behavior 
as-delivered columns that dependent the axis about which the member 
bends. column will usually carry more load about one axis than about an- 
other normal it, these directions being called the strong and weak axes, re- 
spectively. 

This difference behavior may shown follows: rectangular 


bent about the weak axis the quantity Eq. may obtained 
from 


and Eq. would reduce also very nearly true for wide- 
flange section bent about its strong axis, since the web contributes only asmall 
portion the moment inertia and thus the action similar that two 
rectangles. However, for the rectangular section bent about its strong axis 


(a) (b) (c) 


Member free Residual Stress 
Containing Residual Stress 


FIG, 19.—INFLUENCE RESIDUAL STRESS COLUMN CURVE 


and for H-section bent about its weak axis, the term will consider- 
ably less than Thus the buckling strength will less than the value pre- 
dicted Eq. and would computed according Eq. 

Fig. illustrates this difference for idealized case (parabolic residual 
stress pattern, residual stress atflange edges equal -20 ksi, residual stress 
flange centers equal +10 ksi, yield point stress equal ksi, effect 
web neglected.) The lower curve for flexure about the weak axis 
section, while the upper curve for flexure about the strong axis. 

Effect Relationship.—Columns material without definite 
yield level, and with continuously curving stress-strain diagram, with 
aluminum alloys, for example, show influence residual stress for the 
whole range L/r. Although such materials can have much higher buckling 
strengths for shorter column lengths (that is, L/r <40), for the medium 
long columns the effect different stress-strain relationships quite small, 
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FIG, 20.—INFLUENCE RESIDUAL STRESS 
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FIG, 21.—STRESS-STRAIN RELATIONSHIP AND 
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particularly comparison with other effects such variations yield level 
and residual shows column curves calculated for log- 
arithmic, parabolic, and idealized elastic-plastic stress-strain relationships 
for individual fibers. The column cross section rectangular with para- 
bolic residual stress distribution. 

Cold-Bending Residual Stresses.—The study the effect cold-bending 
residual stresses axial column strength indicates that for shortand medium 
length columns, these stresses are more critical than are cooling 
This means that findings based members with cooling re- 
sidual stress patterns will conservative when applied straight members 
whose cooling patterns have been modified cold bending (Fig. 22). 

The Column Curve.—For wide-flange shapes with axially symmetric cool- 
ing residual stresses the solution Eq. (the column curve) may obtained 


Cold Bent 
As-Delivered (Cooling) 


Weak Axis Column Tests 


FIG, 22.—COLUMN CURVES REFLECTING INFLUENCE 
ROLLING AND COOLING AFTER ROLLING 


either from residual stress data, from the average stress-strain curve ob- 
tained stub column Both methods are essentially the same, the 
solution Eq. requiring the function relating and the geometry the 
non-yielded portion the cross section that particular stress, which makes 
possible solution for L/r. The solution will obviously depend the distri- 
bution residual stress. Analytical expressions for this function for certain 
residual stress distributions are available.4,9,11 

The use stub-column obtain the column curve somewhat sim- 
pler approach. the usual situation assumed for loaded wide-flange 
column, namely, that the web does not yield (or that when does the flanges 

“The Influence Residual Stress the Instability Columns,” Huber, 
Dissertation, Lehigh Univ., May, 1956. 


Residual Stresses High Strength Steel,” Feder and Lee, Lehigh 
Univ., Fritz Lab. Report 269.2, April, 1959. 
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have completely yielded), then the following equations will hold true:4 


Since for most wide-flange columns the ratio about Eqs. re- 
duce 


The magnitude isdetermined conducting stub-column test wide- 
flange shape containing residual stresses, plotting the average stress-strain 
curve, and then determining the tangent various stress levels. 

Fig. was drawn using residual stress data (except that the effect the 
web has been neglected). Fig. has been drawn using the stub column stress- 
strain curve Fig. 18. The two methods have been and were 
found very good agreement for steel. Further, the results tests 
also correlate well with them. 

Column Curve Approximations.—A simplification Eqs. might 
desirable for use basis for arriving design formulas. will noted 
from Fig. 23(b) that the curve for buckling the strong direction approxi- 
mately parabolic shape and that for buckling the weak direction 
may approximated straight line. Thus the two curves Fig. 23(b) 
(exact solution), could replaced the following relationships, the first 
which the same form suggested Friedrich Bleich:14 


for L/r 


“Buckling Strength Metal Structures,” Friedrich Bleich, McGraw-Hill Book 
Co., New York, 1952. 
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for and 


for L/r< Eqs. isthe yield level and the proportional 

limit. The latter value could determined either from stub-column tests 

from residual stress measurements using Eq. parabolic column curve, 


tangential the Euler curve, was originally proposed Johnson15 
1893, the basis test results. 


for use Eq. The variation with and flexure axis shown 


Fig. 23(a) andis the information used the “exact solution” curves 
Fig. 23(b). 


Fig. shows column curves for number specimens determined 
the “stub column” method and also the straight line approximation (Eq. 10) us- 
ing average values obtained for this group. (In Fig. 24, many that follow, 
the information presented non-dimensional form order adjust for 
variations and E.) The approximation Eq. also shown which 
all the data Table has been included. While the scatter considerable, 
the straight-line approximation evidently good one. 

Fig. shows the results weak-axis column tests comparison with the 
same straight-line approximation arrived from data from all the stub col- 
umns tested. The circles show the maximum load the columns carried, not 
the point first bending. The results annealed columns are also shown; 
their strength clearly above that the as-delivered members. 

Fig. shows the results column tests and the parabolic approximation 
Eq. using, again, the average properties from all the stub columns test- 
ed. Again there good correlation between theory and test. 

evident therefore that Eqs. 10, based residual stress and tangent 
modulus considerations, could form the basis for design curves. 

this point cognizance has been taken the fact that structural 
steel strain hardens after the plastic strain has reached about ten times the 
elastic limit value, shown the first portion Fig. 27. This results 
diagram that differs from Fig. 19. Instead the yield stress, 


Et= 900 ksi, about shown Fig. the column curve affected 


for low (L/r)-values. Below L/r column will carry greater average 
critical stress than the full yield value. 


and Structures,” Salmon, Vol. Longmans, Green and Co., 
London, 1948, 
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10° 


for L/r 122. 


orderto simplify thedesign procedure eliminating the consideration 
flexure axis, has suggested that 0.5 used with 


parabolic curve, using the single curve for both strong and weak axis bending. 
The resulting “CRC curve” shown Fig. 


general terms 


the past was common explain the reduction column strength 
region L/r 100 asdue accidental eccentricities and initial curvature. 
Such accidental eccentricities, fact, were estimated and appear certain 
design formulas. Although the secant type formula derived the basis 
idealized stress-strain curve with yield point the elastic limit (33,000 
psi), assumes certain value for accidental eccentricity. This latter value 
was arrived analytically correlation with study column tests. Since 
any correlating column tests must have included as-delivered specimens that 
contained residual stresses, the magnitude the accidental eccentricity 
initial curvature must necessarily have been arbitrary since considerable 
portion the reduction column strength now known tobe due the pres- 
ence residual stresses rather than eccentricities. 

Thus, design curves for column strength based the tangent modulus 
method modified the presence residual stress, reflect actual conditions 
rather than reliance assumed irregularities. 


BUILT-UP MEMBERS 


pilot investigation into the influence residual stresses the behavior 
built-up columns has been carried Certain the findings are sum- 
marized herewith. 

Residual Stresses.—Welded columns will have high residual stresses; this 
particularly true the tensile residual stresses. For H-shaped members 
the compressive residual stresses may also high. Fig. shows residual 
stresses welded H-shaped member. The tensile stresses approach the 
yield value and the compressive stress the flange tips averages about 
ksi. These stresses are compared with those rolled shapes and Univer- 
sal plates prior welding, Fig. 30. Fig. shows that compressive 


“Guide Design Criteria for Metal Compression Members,” Column Research 
Council (in preparation). 
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Fig. drawn taking strain hardening into account. Assuming L/r 
the onset strain hardening the column curve, Eqs. become 


for 
Eqs. have been replotted Fig. Eqs. 12, and 14, together with 
column test results. All information presented dimensional basis us- 
mental data has been adjusted proportionately these figures, which repre- 
sent the values most appropriate result this investigation. (Using the 
value 0.70 has been used for these curves. This corresponds 
Orc which the peak(mode) the frequency distribution, Fig. 2.)16 
(Actually, the yield stress level, experimentally, was 33.4 ksi, but this 
close the specification minimum 33.0 ksi that the latter was used the 
average value—particularly the factor safety would account for any devi- 
ation from this average.) 

The test points for Fig. have been plotted from the data contained 
Figs. and 26. For the foregoing average values, take the form 


xx 


for 20< 122, 


(13) 


for 122, and 


“Influence Residual Stress the Strength Structural Members,” 
Ketter, Welding Research Council, Bulletin No. 44, November, 1958. 
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FIG, 27.—INFLUENCE STRAIN HARDENING COLUMN CURVE 
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FIG, 28.—COLUMN CURVES AND TEST RESULTS 
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stresses welded H-shaped members may higher than those rolled 
wide-flange shapes. Compressive residual stresses universal plates 
ksi ksi have been observed flange tips (Fig. 30).7 

Since the magnitude anddistribution welded residual stresses marked- 
influenced the geometry, further workis withcross- 
sectional shapes other than the H-section. would expected that the use 
welded H-shape columns would replaced more frequently the future 
the economical “box” section. 

Riveted built-up columns have considerable variation residual stresses 
that afunction the geometry the component parts. Figs. show 
some measurements that have been made. 

can predicted the same techniques were used for rolled shapes with 
symmetrical cooling residual stresses. The results, correlated with actual 
column tests, are shown Fig. 32. The tangent-modulus concept not rea- 
listic for the prediction these welded columns, but ultimate- 
strength gives good correlation.7 


RIVETED 
COLUMN 
COLUMN 
Welded Specimen 
Specimen 
Tangent Modulus 
Ultimate Strength 
FIG, 31.—SIMPLIFIED 32.—TEST RESULTS AND COLUMNCURVES 
STRESS PATTERN FOR FOR WELDED AND RIVETED MEM- 
RIVETED COLUMN BERS—WEAK AXIS 


General conclusions regarding the column strength built-up members 
(particularly the welded ones) cannot made until further studies are com- 
pleted. Even though the strength the welded H-columns was proportionately 
less than that the riveted as-rolled columns, very important not 
draw the conclusion that welded columns usually will “weaker” than corres- 
ponding rolled shapes. The effect would quite different the cross section 


were box form. Studies the effect cross-sectional form are necessary 
and are underway. 


LOW ALLOY HIGH STRENGTH STEELS 


program tests high strength low alloy steel (ASTM A242) has indi- 
cated that the column strength such steel can predicted inthe same man- 
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ner for steel (Fig. The residual-stress distribution rolled 
shapes the same for steel, the formation residual stress being more 
dependent shape than variation material properties. 


Column Tests (A242) 


FIG, 33.—COLUMN CURVES AND TEST RESULTS 


FOR A242 STEEL COMPARED WITH 
STEEL 


Because the high yield-stress level, the influence residual stress 


columns high strength steel not pronounced steel shown 
Fig. 33. 


SUMMARY 


The strength centrally loaded steel columns may expressed 
the tangent modulus (Fig. 19). This modulus depends the state re- 
sidual stress the member. 

tic deformations that occur during cooling after rolling, after welding, 
cold-straightening operations. For rolled welded members, the part cool 
last usually tensile state residual stress (Figs. and 30). 

Insofar H-columns are concerned, the most important the stresses 
are those the flange tips; for rolled shapes the average compressive stress 
there about ksi (Table 2). 

The effect residual stress onthe stress-strain relationship tolower 
the proportional limit and cause the stress-strain diagram non-linear 
beyond that point and the yield stress level (Figs. and 18). 


The pro- 


portional limit for the shapes studied was about ksi (62%of yield stress 
level). 


Residual stresses reduce buckling strength because early localized 
This reduction 


yield that occurs certain portions the cross section. 
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greatest when the slenderness ratio between and and amounts about 
25% (Figs. and 26). The critical stress the inelastic range 
function the moment inertia the elastic portionof the cross (Eq. 
and Fig. 4(d)). 

There difference behavior as-delivered columns that de- 
pendent the axis about which the member bends. Columns given slen- 
derness ratio bent about the weak axis will carry lessload than columns the 
same L/r bent about the strong axis (Fig. 20). This because the highest 
compressive residual stresses are found the tips the flanges. 

Factors such upper yield point, strain rate, and difference strength 
between web and flange cause the yield stress level full cross section 
(stub column) about than the yield point determined routine 
ASTM acceptance test From “mill test” average all data avail- 
able date about ksi, there obtained probable compressive strength 
about ksi for steel. 

Approximate column formulas which may adapted fordesign use have 
been developed which agree well with test results (Fig. 28) and which also re- 
flect the influence strain-hardening (Fig. 27). desired take into 
account the effect flexure axis, then two formulas would necessary (Eqs. 
and 13). Otherwise single curve could used (Eq. 15) which provides 
good approximation the theoretical and experimental results. 

Cold-bending residual stresses are more critical than are cooling 
stresses. Thus, findings based members with cooling patterns may ap- 
plied rolled wide-flange columns (Fig. 22). 

10. Columns built-up welding will contain tensile residual stresses close 
the yield point(Fig. 29). Compressive stresses may higher lower than 
those that form due cooling, depending onthe geometry the cross section. 
Although tests H-shaped welded members exhibit strength that com- 
paratively less than that corresponding wide-flange shape, might ex- 
pected that welded “box” columns would somewhat stronger than corres- 
ponding rolled member. 

11. Low-alloy, high-strength steels contain residual stresses whose magni- 
tude and distribution are similar those steel. Since the formation 
residual stress does not depend the same extent the magnitude the 
yield stress level does the geometry, the influence residual stress 
the strength columns higher strength steels notas pronounced 
columns steel (Fig. 33). 

12. Column strength affected the presence residual stress, and 
dependent its magnitude and distribution. Columns riveted and high- 
strength steels, for example, with proportionately smaller magnitudes com- 
pressive residual stress, have greater strengths than columns welded 
shape steel (Fig. 34). 

13. Residual stresses are the major factor affecting the strength central- 
loaded columns the intermediate ranges (40 L/4 120), and design 
curves for column strength based the tangent modulus concept, modified 


the presence residual stresses, reflect actual conditions more closely than 
has previously been possible. 
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FIG, 34.—DIAGRAMMATIC COLUMN CURVES FOR 
A242 AND ROLLED SHAPES AND FOR 
RIVETED AND WELDED COLUMNS 
H-SHAPE—WEAK AXIS 


NOMENCLATURE 


cross-sectional area: 
web area; 
width flange wide-flange shape; 
depth wide-flange shape; 
Young’s modulus elasticity: 
tangent modulus; 
moment inertia: 
moment inertia the unyielded (elastic) part; 
total length pin-end column; 
ratio; 
L/r effective slenderness ratio; 


standard deviation, statistical measure the scattering obser- 
vations; 
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flange thickness; 
directions for flexural axis; 


unit strain; 


strain rate the plastic range; 
stress: 


applied average maximum stress column; 
stress proportional limit; 

yield stress level; 

static yield stress level (zero strain rate); 


“dynamic” yield stress level strain rate other than zero); 
residual stress flange tips; 


residual stress web center; and 


residual stress flange center. 


DEFINITIONS 


the process for any structure part structure 
pass from one deflection pattern into another without change load. 

Critical Load.—The maximum load column will carry. not coinci- 
dent with the buckling load for axially loaded column. 

Plastic Strain Rate.—The strain rate the plastic range. For coupon 
test, the plastic strain rate, unlike the elastic strain rate, independent 
the elasticity the testing machine. 

Stub Column.—A short compression test specimen, sufficiently long for use 
measuring the stress-strain relationship for the complete cross section, 
but short enough avoid buckling the elastic and plastic ranges. 

Yield Point.—The first stress material, less than the maximum attain- 
able stress, which increase instrain occurs without stress. 

Yield Stress.—The stress which material exhibits specified limiting 
deviation from the proportionality stress strain. 

Yield Stress Level.—The average stress during actual yielding the plas- 


tic range. remains fairly constant for structural steel provided the strain 
rate remains constant. 
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LATERAL BUCKLING BEAMS 


FOREWORD 


The papers that form the ASCE-CRC Symposium Metal Compression 
Members are representative ofall phases Column Research Council activi- 
ties. Several papers represent summaries both past and current research, 
prepared specifically the preparation the “Guide Design Cri- 
Metal Compression Members.” Also included are reports recent- 
completed research pony truss bridges and plate girders. 

the January 1960 Journal the Structural Division Bruce Johnston 
documented the inception the CRC and the many activities undertaken the 
Council. 

Thomas Kavanagh, the February 1960 Journal the Structural Di- 
vision, reviewed the state knowledge concerning the effective length fram- 
columns. 

the April 1960 Journal the Structural Division, Robert Barnoff and 
William Mooney studied the effect floor-system participation pony 
truss bridges. 


Bruce Johnston presented digest the “Guide” the Jour- 
nal the Structural Division. 

the July 1960 Journalof the Structural Division Lynn Beedle and Lam- 
bert Tall provided summary investigations the strength centrally 


loaded columns influenced residual stresses and variations the yield- 
stress level. 


Note.—Discussion open until December 1960. Separate discussions should sub- 
mitted for the individual papers this symposium. extend the closing date one 
month, written request must filed with the Executive Secretary, ASCE. This paper 
part the copyrighted Journal the Structural Division, Proceedings the Ameri- 
can Society Civil Engineers, Vol. 86, No. July, 1960. 

Asst. Chf., Engrg. Design Div., Alcoa Research Labs., New Kensington, Pa. 
Chf., Engrg. Design Div., Research Lab., New Kensington, Pa. 
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SYNOPSIS 


This paper intended provide brief summary information that may 
used background material writing specifications for the design 
beams and girders whose strength controlled lateral buckling. 

general formula developed for elastic buckling strength beams, with 
coefficients depending principally the conditions loading and support. 
Values the coefficients this formula are tabulated for wide variety 
conditions loading and support. The relationship between theoretical buck- 
ling strength and experimental results discussed. 

Various design formulas that have been proposed for beams are compared 
with each other and with experimental and theoretical beam strengths. 


INTRODUCTION 


This paper intended provide brief summary information that may 
used background material writing specifications for the design 
beams and girders whose strength controlled lateral buckling. The pa- 
per not intended recommend specific design procedures. not possi- 
ble cover all the available material this subject paper reasonable 
length. For the sake brevity, detailed information the form charts and 
tables presented only cases where the data are not readily available 
published form. Much the material this paper taken from report 
which the authors the request Research Committee Col- 
umn Research Council (40). 

Notation.—The letter symbols adopted for use this paper are defined 


where they first appear and arranged alphabetically, for convenience re- 
ference, Appendix III. 


BUCKLING STRENGTH BEAMS AND GIRDERS 


Elastic Buckling.— for the elastic buckling strength beams can 
used basis point departure the preparation design rules. 
The results many the theoretical investigations that have been made 
the elastic buckling strength beams through 10) can expressed the 
following equation, which derived Appendix 


which the are coefficients, refers the critical stress for lateral 
buckling, the modulus elasticity, denotes the moment inertia 
the beam cross section about the y-axis, coefficient, Lis the distance 
between points support against lateral bending and twisting, denotes the 
torsion warping constant, refers the distance from shear center point 
application transverse when load below shear center and 


Numbers parenthesis thus, (1), refer corresponding items Appendix II.— 
Bibliography. 
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negative otherwise), the modulus rigidity, the torsion constant, and 
denotes the section modulus for stress compression flange. The factor 


The coefficients Cg, and depend mainly the conditions load- 
ing andsupport for the beam. Values these coefficients for number dif- 
ferent cases are listedin Table These coefficient values are limited 
prismatic beams bent about principal axis. cases where transverse loads 
are applied, assumed that the lines action the loads pass through the 
shear center and the centroid and that the loads are attached the beam 
such way that their lines action remain parallel their initial directions 
the beam deflects. also assumed that the shear center lies prin- 
cipal axis through the centroid. 

Elastic buckling theory has been confirmed tests for several types 
loading and shapes cross section. Many these tests have been performed 
aluminum alloy members. The aluminum beams tested include symmetri- 
cal I-beams under uniform bending moment (12), unequal end moments (13), 
and under concentrated load the center the span (14); and deep rectan- 
gular beams (43), channels (15), Z’s (15), and unsymmetrical I-sections (16) 
under uniform bending moment. Tests structural carbon steel beams sub- 
jected two concentrated loads the top flange the quarter points have 
also demonstrated the validity elastic buckling theory (17). 

Inelastic Buckling.—Bleich (18) has pointed out that possible obtain 
lower limit the theoretical buckling stress the inelastic range sub- 
stituting the tangent modulus corresponding the maximum stress the 
beam for the elastic modulus the elastic buckling formula. Tests alu- 
minum alloy beams have shown that this method gives close approximation 
the experimental buckling stress when the bending constant along 
the length (12)(13). Tests aluminum alloy beams subjected unequal end 
moments, with the ratio the end moments varying from 1.0 -1.0, resulted 
experimental critical stresses below 39% above the val- 
ues computed the method suggested Bleich (18). 

showed the tangent modulus curve for structural steel suggested Bleich (18) 
somewhat unconservative comparison with some ofthese tests results. 
However, probable that tangent modulus curve based compressive 
tests the full cross section the members used the beam tests would 
have been better agreement with the test results, Results Austin’s study 
are shown Fig. 

second method estimating the effect plastic action onthe strength 
beams and girders assume that the relationship between the elastic and 
inelastic buckling strength the same for beams for columns. Inelas- 
tic buckling strength beams may then estimated from column curve. 
This procedure applicable aluminum alloy members since the tangent 
modulus buckling curve agrees well with tests both columns (23). 
This method has also been recommended basis for design steel beams 
(8)(39). 

third method taking into account the effect inelastic behavior 
determine the relationship between elastic buckling strength dir- 
ectly from tests beams. This procedure takes into account the effect any 
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TABLE 1.—VALUES OF COEFFICIENTS IN FORMULA FOR ELASTIC BUCKLING STRENGTH 
OF BEAMS 


Condition of 
Restraint 
Against Rota- 


tion about Values of Coefficients 
Bending Moment Vertical Axis 
Loading Diagram 


1.30-1.32% 


Fixed Salvadori (4) 


Simple Support Salvadori 
Fixed 


Fixed 

Fixed Austin, et al(5) 


Simple Support s Austin, Yegian 
and 
Fixed 


Simple Support Pettersson(6) 
Fixed Austin, et al) 


Simple Support A Austin, Yegian 


and Tun 
Fixed 


Simple Support J Winter, 
Schrader(10) 


Cantilever Beams 


Ti moshenko, (1) 
Kerensky, 
Flint and 
Brown(®) 


Poley!?) 


@ Minimum value applies to beams with negligible resistance to warping and can be used to obtain conserva- 
tive buckling stresses for all beams. 


b Average value giving buckling stresses within about 10 per cent of values computed by Pettersson in all 
cases. 


© For cantilever beams, no single combination of values of K and C, applies to all proportions of beams. In 
Case No. 11, the variation of C; corresponding to a single value of K is approximately a minimum for k = 1.0. 
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No. Data 7 
Beams Restrained Against Lateral Displacement at Both Ends of Span 
om 
3. = 
w 
Supported End 
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BEAMS 
affected the test results. 


beam tests the University Washington shown Fig. 
Buckling Beams Under Various Special Conditions. 


Oblique Loads.—The case ofa beam subjected uniform bending moment 
that does not lie one the principal planes the cross section has been 


The latter showed 


discussed Goodier (2) and McCalley (20). 
that the equation for the critical moment takes the form Eq. with 


and denote principal axes, and the x-axis the axis normal the plane 
bending. 


@ 101 25.4 BEAM 
© 27 BEAM 
0 18 50 BEAM 
4 106 '5 BEAM 


25 


20 


APPROX. CURVE 
OF LOWER BOUND 


BUCKLING STRESS BY TEST, KS! 


10 
THEORETICAL ELASTIC CRITICAL STRESS, KSI 


FIG, BENDING TESTS ROLLED 
STRUCTURAL CARBON STEEL 


Nonuniform Cross Section.—The effectonthe buckling strength variations 
the flange thickness along the length I-section girder has been investi- 
gated Kerensky, Flint and Brown (8), who have suggested 
basis design for such members. Gatewood (21) has shown that for can- 
tilever beams with negligible warping rigidity, variations the cross section 
have the same the buckling strength cantilever columns. 
concise solution for the lateral buckling strength tapered narrow rec- 
tangular beam, subjected constant bending moment and simply supported 
the ends, has been presented Lee (41). Tapered cantilever I-beams 
were investigated experimentally Krefeld, Butler and An- 
derson (42). 

Special End Conditions.—Austin, Yegian, and Tung (5) gave solu- 
tions for elastic buckling strength I-beams under transverse load (either 
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other factors such residual stress and initial crookedness, which may have 


The curve for inelastic buckling strength struc- 
tural carbon steel beams that has been suggested Austin the basis the 


BLEICH TANGENT 
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uniform load concentrated load the center) with various degrees re- 
straint against rotation the beam about either vertical transverse 
horizontal axis the ends. Each type restraint was considered vary be- 
tween zero and complete fixity. all cases, the beams were considered 
fixed against rotation about longitudinal axis the ends. Tests Hechtman, 
Styer and Tiedemann (17) demonstrated that web-angle connections 
the ends beams prevented rotation about longitudinal axis and restrained 
rotation about transverse horizontal axis, but did not provide much restraint 
against rotation the flanges about vertical axis. Top and seat angle con- 
nections resultedin higher beam buckling loads than did web angle connections 
since the former provided some restraint against rotation the flanges about 
vertical axis, well restraining rotation about the other two axis. 

Frequently cantilever beam simply the overhanging end beam that 
extends over two more supports. this case, the supported endof the can- 
tilever beam not fixed against lateral bending the beam flanges but some 
restraint supplied continuity the support. such cases, conserva- 
tive estimate the buckling strength can made considering the warping 
constant, Cw, zero the buckling formula. 

beam supported both ends continuous one both supports, the 
end conditions for any one span are generally between the cases complete 
fixity andsimple support covered Table The effect continuity has been 
discussed Salvadori (22) and Nylander (3). 

Large for which the ratio relatively large 
may develop appreciable curvature the plane the applied bending moment 
before lateral buckling takes place. some cases, this curvature can have 
appreciable effect the buckling strength. Pettersson (6) has devel- 
oped accurate expression for the increase buckling strength that results 
fromthe deflection simply-supported I-beam subjectedto uniform bend- 
ing moment. However, Clark and Knoll (24) have shown that, when 
the ends the beam are restrained against lateral rotation, the buckling 
strength not appreciably increased the deflection prior buckling. Ini- 
tial camber beam may, course, cancel any effect deflection the 
buckling strength. 

Reduced Torsional Stiffness Plate Girders.—The torsional stiffness 
plate girders may reduced because warping the cross sections and al- 
result slipping between the various parts built-up sections. 
Chang and Johnston (25) have shown that computing the torsional stiff- 
ness riveted plate girders, accurate value the torsional constant, 
can computed considering the parts the flange between the outside 
lines rivets single solid piece, and the parts extending beyond the out- 
side lines rivets individual elements. The tests Chang and Johnston 
(25) showed that, the rivets are sufficiently close together, this method 
computing causes maximum error about 20%, which would result 
error not more than 10% the computed elastic buckling stress. Ny- 
lander (26) has demonstrated that distortion the cross section not very 
significant for beams proportions. For very deep plate girders, 
the tests Chang and Johnston (25) indicated that the shape the cross sec- 
tion maintained partly through the effect the stiffeners. Winter has 
recommended (9) that for beams with very thin webs and for latticed trusses, 
the torsion constant should takenas zero computing the buckling strength, 


since this assumtion does not cause mucherror the buckling stress for these 
members. 
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Eccentric Loads and Crookedness.—The behavior beams under eccentric 
loads has been investigated Pettersson (6) and the effects initial bow and 
initial twist have been studied Kerensky, Flint and Brown (8). general, 


these factors affect the buckling strength the same way that they for 
columns. 


NOMOGRAPH FOR BUCKLING STRENGTH BEAMS 


Eq. appears first glance too involved for use design tool. 
However, all the quantities except are tabulated for many structural 
sections and may computed for other sections (27)(28). The quantity can 
approximated the quantity for which formulas are readily available 
The arithmetic involved calculating critical stresses from Eq. (or 
values “equivalent slenderness ratio” the equation put form that 
can used with column buckling curve) can, greatly reduced with the 
aid nomograph, such that shown Fig. Fig. gives the equivalent 
slenderness ratio L/r) corresponding Eq. for unsymmetrical 
beam under uniform bending moment. For symmetrical sections, the values 
L/r from Fig. can multiplied make them applicable 
beams with other distributions bending moment. For unsymmetrical sec- 
tions under nonuniform bending moment, the quantity should multiplied 
take into account the effect applying transverse loads above be- 
low the shear center, the quantity Fig. should replaced 
+C3 ly. Again, the resulting value KL/r must multiplied 
make Fig. applicable beams whose section properties are outside the 
range coveredinthe chart, the values Cw, for the 

section can multiplied constant such 0.1 0.01. Use Fig. with 
these modified section properties (the true values must used) 
will give the correct value L/r for the beam. For any particular materi- 


al, the answer obtained from Fig. could put the form allowable 
stress rather than equivalent slenderness ratio. 


SIMPLIFIED FORMULAS FOR BUCKLING STRENGTH AND 
ALLOWABLE STRESS 


There are many that result simplifications Eq. and al- 
many approximations that canbe used reduce the complexity the com- 
putations. 

Symmetrical Sections.— For sections that are symmetrical about the hori- 
zontal axis about point, the quantity Eq. zero. Further 
simplication results there transverse load, the transverse loads 
are applied the neutral axis, which case The expression for elas- 
tic buckling strength can then written 


fer 


Large copies this chart will furnished upon request Clark, Alcoa 
Research Laboratories, Box 772, New Kensington, 
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Eq. can put the form 


which function the conditions loading and restraint, the modulus 
ofelasticity and Poisson’s ratio ofthe material andofthe parameter 


EFFECTIVE LENGTH, KL, IN. 
20 40 60; 100 200 


TORSION WARPING 


TORSION CONSTANT, 


00 59 MOMENT OF 
INERTIA ABOUT 


AND COMPRESSION 
FLANGE 


- noe 


EFFECTIVE LENGTH, 
KL, IN. 


- Noe 
o 


SECTION MODULUS, 


EQUIVALENT SLENDERNESS RATIO 
FOR BEAM UNDER UNIFORM 
BENDING MOMENT, KL/ 


NOTE: DOTTED LINES SHOW EXAMPLE IN WHICH KL «72 IN., J*0.124 Cy +3047 INS 
ly INS 4.1 1N., So#125 INS 


3, (** et ) 


2.—NOMOGRAPH FOR EQUIVALENT SLENDERNESS RATIO, KL/r, 
BEAMS UNDER UNIFORM BENDING MOMENT 


Eq. can used basis for design suitable factor safety intro- 
duced and allowance made for inelastic action. The equation can supple- 
mented series ofcurves giving values for number conditions 
loading and restraint plotted against Cw/(L2 J). Values the cross section 
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properties and are tabulated for standard aluminum I-, Z-, and channel 


sections (27). Values and Cw/G for rolled steel sections have al- 
been published (28). 


Symmetrical I-Sections. 
General Formula.—For doubly symmetrical I-beams, 


which the distance between centroids beam flanges (taken depth 
beam approximate formulas). Eq. can written 


Values the coefficient can plotted against the parameter (L2 J), 
for different conditions loading and restraint. approach this kind 
used German specifications (29). additional modification the coeffi- 
cient Eq. can introduced take approximate account the effects 
transverse loads being applied the top bottom flange rather than the 
neutral axis (23). 

Parameters L/ry and t/ry d.—Eq. can further simplified intro- 
ducing approximate expressions for the section properties the beam. 
this procedure, the effect the web the section properties generally ig- 
nored, and the flanges are considered rectangular sections. Expressions 
this kind were derived Winter and other discussers Vries’ paper 
equation this type, which has been proposed basis for new 
British specifications (8), the following: 


which the thickness the compression flange and denotes the radius 
gyration the beam cross section about the y-axis. The proposed value 
the coefficient for structural carbon steel 381,000,000 psi (8). 

Parameter t.—An additional simplification results cases which 
the term t/ry large comparison unity, that the first term 
under the radical Eq. can neglected. When expressed terms 
the flange width the resulting expression for the critical stress 


Eq. was developed Vries (30), who selected value 20,000,000 psi 
for the coefficient C7. Vries has shown that within the elastic range, Eq. 
reasonably accurate andgenerally conservative the buck- 
ling strength rolled steel I-beams usual proportions under number 
loading conditions (30). This has been further substantiated Hechtman’s 
tests (17). Eq. was used basis for the present American Institute 

Steel Construction (AISC) specification (31). 
rough idea the degree approximation involved Eq. can ob- 
tained comparing with Eq. which quite accurate for symmetrical 
beams, demonstrated subsequently the results Figs. 8for val- 
Ac/(Ac At) 0.5. For beam sections given material, the maxi- 
mum difference between buckling stresses corresponding Eqs. and oc- 
curs the transition point between the elastic and plastic buckling curves. 
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For many materials, this transition point occurs stress the vicinity 
two-thirds the compressive yield strength, fy. Therefore, find the ap- 
proximate maximum difference between Eqs. and for given material, one 
can set the critical stress Eq. equal 2/3 fy, find the proportions 
beam having this critical stress, and determine the critical stress computed 
for this beam from This procedure results the following equation: 


The following relationships were used obtaining Eq. 13.1 (8); 
0.73 (or 21,200,000 psi for steel) and 4.2 (8). 

The ratio stresses given Eq. plotted against d/t for various ma- 
terials Fig. Fig. shows that the formula gives good accuracy 
for A-7 carbon steel and I-sections, but can considerably error for 
structural shapes other materials. For plate girder sections all mate- 
rials, Eq. can give appreciable errors. 

Parameter L/ry relatively deep beams, the quantity 
may small that Eq. can approximated the expression 


for 


formula ofthe above been presented the American Iron and Steel 
Institute (AISI)(35) and Priest and Gilligen (32). For simple 
structural shapes, the value approximately proportional the flange 
width, that can expressed terms the parameter L/b instead 
L/ry. Formulas involving L/b are usedin the American Assn. State High- 
way Officials (AASHO)(33) and American Railroad Engineers Assn. (AREA)(34) 
specifications for steel bridges. 

his discussion Vries’ paper 1947 (30), Winter mentioned that two 
sets design formulas could used, one set involving the parameter L/r 
and the other set involving the parameter Ld/bt. questionable cases, al- 
lowable stresses calculated from both sets offormulas and the higher 
allowable stress would used. Recently, Austin has made similar sugges- 
tion (36), incorporating the parameters L/b and The theoretical basis 
for this proposal illustrated Fig. graph originally presented Aus- 
tin (36). The curves labeled and Fig. are the theoretical buckling 
stress curves for beams under uniform bending moment top flange 
loading, respectively. The diagonal line labeled Fig. can approxi- 
mated Eq. and the horizontal line labeled can approximated Eq. 


BEAMS 
The specific design formulas suggested Austin are follows: 


But not less than 


12,000,000 
884, 


which the allowable stress. 

the formulas involving the parameter Ld/bt, the quantity would 
replaced for riveted plate girders and for welded plate 
girders, has been suggested Vries (30). the formula involving L/b, 
girders. avoid the necessity computing allowable stresses both sets 
formulas, could specified that Eq. would used except when the val- 
(L/b)2is less than the value Ld/bt. This limitation can derived 
approximately from Fig. 

Unsymmetrical I-beams symmetrical about the vertical 
axis but unsymmetrical about the horizontal axis, subjected uniform bending 


moment, Hill (16) derived the following approximate equation for the elas- 
tic buckling stress: 


Eq. the same Eq. for this case, except that the quantity approx- 
imated the distance from the shear center the centroid, 


Winter (9) derived the following approximate equation for the same case: 


Eq. can derived from Eq. replacing the quantity the expression 
The foregoing equation has been further simplified the proposal for the new 
British specifications (8) the following pair equations: 


for) 
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ALUMINUM T-I ALLOY A-8 NICKEL A-7 CARBON 
606!-T6 STEEL STEEL STEEL 


AND d/t (Eq.6) BUCKLING STRESS 
BASED Ld/bt (Eq.7) 


APPROX. RANGE OF d/t VALUES FOR 
ROLLED W AND I SECTIONS 


MAX. RATIO BEAM BUCKLING STRESS BASED 


RATIO, d/t 


3.—COMPARISON Ld/bt BEAM FORMULA WITH 
FORMULA BASED L/ry AND d/t 


N 
xi > 


UNIFORM LOAD, 
TOP FLANGE 


STRESS PARAMETER 


RATIO 


FIG, 4.—CRITICAL STRESS CURVES FOR I-BEAMS 
SUPPORTED THE ENDS 
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which and are the distances from the centroid the cross section 
the centroid the compression and tension flanges, respectively. 

different approach the design unsymmetrical sections has been sug- 
gested Vries (38), who the formulas involving the par- 
ameter d/b used, with the quantity modified follows: 


For Riveted 
Plate Girders: 
For Welded 
Plate Girders: 


Calculations elastic buckling stress, using the foregoing formulas, have 
been made for two basic girder sections, considering several different sets 
flange proportions for each girder. These basic cross sections are the same 
those used the recent study Kerensky, Flint, and Brown (8). How- 
ever, the values length used the present computations are somewhat dif- 
ferent from those used the British study. Two different length values were 
assumed for each girder. For one set computations, the two flanges were 
assumed the same thickness but different widths. another set com- 
putations, the two flanges were assumed the same width but different 
thicknesses. The results these calculations are shown Figs. 8.5 

material considered structural carbon steel, critical stresses above 
the level 33,000 psi are not practical interest. Higher values may 
interest, however, for high-strength steels sluminum. the charts that 
apply the girders having flanges equal widths but different thicknesses, 
vertical line drawn value the ratio Ac/(At Ac) corresponding 
width-to-thickness ratio for the outstanding compression flange. Gird- 
ers with values Ac/(At Ac) less than the values corresponding these 
vertical lines would not permitted AISC specifications (31). 

Although the comparisons made are limitedto specific cross 
sections and lengths, fairly wide range proportions covered. should 
noted that allthe formulas that are compared Figs. are intended 
apply beam subjected uniform bending moment, with the exception 
Eqs. and 14. The constant 20,000,000 used Eq. was recommended 
Vries for beams uniformly loaded the top flange (30). However, several 
discussers Vries’ paper showed that theoretical considerations the 
uniform bending caseled formula with coefficient the vicinity 
20,000,000. 

somewhat higher value for uniform bending can derived fol- 
lows. Vries recommended value 24,000,000 for beams subjected 
uniform transverse load the centroid. The theoretical critical stress this 
case 13% greater than for the case uniform bending, shown the coef- 
ficients Table Therefore, reasonable value the constant for uniform 
bending might 24,000,000/1.13 21,200,000. this value had been used, the 
curves Figs. for Eqs. and would have been about higher. 

Kerensky, Flint, and Brown (8) compared critical stresses computed from 
Eq.6 with the “exact” solution for some built-up sections incorporating angles 


was assumed 29,4000,000 psi and Poisson’s ratio 0.25. 
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and channels the flanges well sections the type considered Figs. 

and applies channel and zee sections subjected 
end moments. Although the accurate expression for for these shapes 
somewhat more complicated than Eq. Hill has pointed out (15) that Eq. can 
used computing the buckling strength standard channels aid 
Eq. the resulting errors critical stress not exceeding about 6%. Eq. 
can used approximation for channels and zees subjected lateral 
loads, provided that the loads pass through the shear center (in back the web 
the case the channel) that they not tend rotate the section. 

Box and Multiweb Sections.— applies box sections and other shapes 
with more than one web well the more conventional beam shapes. For 
box section beams, the torsion constant usually large that lateral buck- 
ling need not designingthe beam. Hat sections andother open, 
double-web shapes frequently have sufficiently large values and that 
there little tendency fail lateral buckling. 


CONC LUSIONS 


Elastic buckling strength beams can computed for wide variety 
cross sections and conditions loading and support. Some this informa- 
tion summarized Table which gives values coefficients used 
Eq. 

Elastic buckling theory has been verified number experimental 
investigations aluminum and steel beams. 

Inelastic buckling strength can predicted from the elastic buckling 
strength with the aid relationship determined from beam tests from 
column curve. 

Computation beam buckling strength can simplified means 
nomograph such Fig. which can applied unsymmetrical cross sec- 
tions and variety conditions loading and support, including transverse 
loads applied points above below the shear center. 

number simplified formulas have been developed for designing 
beams the basis their buckling strength. Some these methods and the 
degree approximation involved them have been related this paper. 
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APPENDIX 


GENERAL EQUATION FOR LATERAL BUCKLING STRENGTH BEAMS 


Consider prismatic beam bent about principal axis. transverse loads 
are applied, assume that their lines contains aprin- 
cipal axis through the shear center and the centroid. 

The potential energy such beam (11) buckles 


which the displacement parallel the initial position the x-axis and 
Bis the rotation the cross section 
Equilibrium considerations require that 


Ely u"=MB. 


The load can expressed terms the moment, 


Use Eqs. and eliminate and Also, make the follow- 
the bending moment the beam and function determine the 
buckling load, set Eq. becomes: 
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and, 


and solving Eq. for one obtains the equation 


The values Cj, Cg, and depend the moment function 
The function that describes the boundary conditions and but 
usually not affected greatly the beam cross section and length. Hence, 
the coefficients Cg, and depend the boundary conditions and 
the loading, but are relatively independent the proportions the beam. This 
demonstrated the fact that numerical solutions for the buckling strength 
beams under variety loading conditions can expressed the form 
Eq. 20, which values the coefficients depend only onthe conditions load- 
ing and support (Table 1). 
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APPENDIX --NOTATION 


The symbols used this paper are defined follows: 

flange width symmetrical I-section channel 

width compression flange 

width outstanding portion compression flange 

width tension flange 

coefficients 

torsion warping constant 


distance between centroids beam flanges (taken depth beam 
approximate formulas) 


distance from shear center shear center lies be- 
tween centroid and compression flange, negative otherwise) 


modulus elasticity 


tangent modulus 


a 
C1, 
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fer stress for lateral buckling 

allowable stress 

compressive yield strength 


distance from shear center point application transverse load 
(positive when load below shear center; otherwise negative) 


modulus rigidity 

moment inertia compression flange about axis parallel web 
moment inertia tension flange about axis parallel web 
moment inertia beam cross section about x-axis 

moment inertia beam cross section about y-axis 

major moment inertia (when different from Ix) 

minor moment inertia (when different from ly) 


torsion constant 


coefficient 


distance between points support against lateral bending and twisting 
bending moment any point beam 
maximum value for given beam 
transverse load unit length beam 
radius gyration beam cross section about y-axis 
section modulus for stress compression flange 
thickness compression flange 
displacement parallel initial position x-axis 


potential energy (increase strain energy during buckling minus work 
done external forces) 


coordinate through centroid normal axis beam and normal plane 
bending 


coordinate through centroid normalto axis beam and parallel plane 
bending, positive downward 


distance from centroid ofcross sectionto centroid compression flange 


distance from centroid cross section centroid tension flange 


coordinate along centroidal axis beam 
=2/L 


rotation cross section 
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DISCUSSION 
WIND FORCES STRUCTURES: NATURE THE 


Closure Robert Sherlock 


ROBERT SHERLOCK, ASCE.—Appreciation expressed those 
who have presented discussions this paper. Several the discussions have 
dealt with the variation wind velocity with height coastal areas during hur- 
ricanes, subject whichadequately controlled observations are indeed mea- 
gre and where code recommendations must, therefore, based inferences 
what little available. result these discussions, has be- 
come clear that separate tables must presented for design velocities 
coastalandinlandareas. The one-seventh- power the variation 
wind velocity with height inland stations has been retained, was done 
Table for basic velocities (fastest mile wind level) mph, 
but the assumption, based paper Mr. Bunting (6), that the variation 
wind velocity with height could neglected mature hurricanes, has been 
discarded. Curves and Fig. (which was entitled Fig. must, 
therefore, discarded, must the design velocities Table for basic 
winds above mph. The reasons for this drastic change will appear subse- 
quently. 

Another drastic change that Fig. has been superseded two maps, 
Figs. and which first appeared July, 1959 manuscript the United 
States Weather Bureau and have since been published Proceedings, 
ASCE, April 1960 (9).“ They were prepared the Office Climatology 
the USWB under the direction Thom. The maps show contour 
lines, each which gives the fastest mile wind above the ground 
the points through which the contour line passes. The contour lines were 
based upon average about USWB records each 141 airport 
stations throughout the United States. means statisticalanalysis the rec- 
ords were extrapolated velocities having 50-yr probable period recur- 
rence Fig. and 100-yr probable period Fig. should noted that 


Note.—This paper part the copyrighted Journal the Structural Division, 
Proceedings the American Society Civil Engineers, Vol. 86, No.ST March, 1960. 
July, 1958, Robert Sherlock. 


Numerals parentheses, thus (1), refer tocorresponding items the Bibliography. 


197 


2563 
= 


emjen(pe ‘speeds para 

3 SUC; 

‘peqeaete st eansodxe 
242 JI GZ jo 
& 403 BsOzTUN AT OATIETOS 
8} 
eunsodxe 10; speeds 


\ 

\ 


198 July, 1960 
\ 
\ = 
Ye \ og, 
4 : \ ™ No | 
\ - 4 L . 
LS \ | 
Qi ‘ \ \ - | 
| / ua 
\ 
| 4 ! 
4 
4] / 
A ‘oN Vo S | 


os 
es dea 03 epee 
eq ‘speeds puys 


40 

242 GZ jo 
e@ansodxe 40; speeds 


199 
* o+ 
\ 


200 July, 1960 


these are “probable” values and that, over long period time there may 
occur interval which may include two more basic winds having 
“probable period recurrence” 100 yr. Nevertheless, this method gives 
the designing engineer “conservative” anda “more conservative” basis for 
choosing the wind velocity which appropriate the design each type 
structure, taking intoaccount the wind-sensitivity the structure, and the haz- 
sensitivity depends the importance wind loading compared the im- 
portance allother kinds their effect the strength and stabil- 
ity the structure. 


FIG, 7.—SIMPLE WIND ROSE 


Mr. Gilbert the typographicalerrors whichunfortunately occurred 
the original printing this paper. The corrections were printed the Sept- 
ember, 1958, issue the Journal the Structural Division, pp. 1787-71. 

Mr. Gilbert asks the writer comment the design usefulness “wind 
roses,” since they are sometimes available engineers working outside the 
United States where detailed weather data are not readily available. The an- 
swer that wind roses may may not useful the design engineer who 
wishes select the fastest-mile wind, the fastest-1-min velocity, for use 
his design. depends the type wind rose and the amount supple- 
mentary information available him. 
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Fig. very simple wind rose given the form bar diagram and 
indicating the percentage time when the wind has blown from the different 
points the compass some previous sample wind history. For example, 
this diagram indicates that, for the 5-yr sample shown, the winds Grand 
Rapids blew from the southwest for about 8.6% all the hours. Wind velocities 
are not given. 

Fig. type wind rose which eachof the bars has been sub-divided 
into the percentage time when the wind velocity lay between certain limits. 


0-3 4-7 13-18 
Wind Speed MP.H. (Hourly Average) 


FIG, 8.—MODIFIED WIND ROSE 


Usually, the diagram basedon the hourly average obtained from the USWB 
some local installation instruments. This wind rose shows that about 
0.3% all winds blew from the southwest with hourly average velocity be- 
tween mph This little use the designer unless able 
convert the maximum hourly average wind velocity into distribution from 
which can obtain the fastest-mile wind some other subdivision which 
needs. For this purpose may approximate the answer using published 
information from other locations which gives the distribution the average 
wind velocity for shorter periods time recorded longer period (10). 


LEGEND 
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The writer has assumed that Mr. Gilbert’s question referred tothe simpler 
forms the many types wind roses which are available. type which in- 
volves more complicated statistical concepts and which was developed for use 
planning aircraft flights described Crutcher (11). 

Mr. Scruton would like see the subject the maximum gust period, which 
significant various design problems, more detail. Mr. Daven- 
port raises the same question regard sharp edgedstructures. Un- 
fortunately, the writer knows progress which has been made this sub- 
ject beyond that quoted previous paper (12) under the heading “Minimum 
Effective Gusts.” 

Mr. Perrin points out several problems, connection with the choice 
wind velocities and gust factors for the design guyed tower, which the 
judgement the engineer may supplement modify the specified Code 
Manual. The writer has limited his presentation material widespread 
usefulness structural design the assumption that the more highly spe- 
cialized structural problems will receive the attention engineers with spe- 
cialized competence make direct such modifications. has had 
the problem limitation, since the number special problems very great. 
example such special problem the selection appropriate wind ve- 
locities used for the design the bowl, supporting structures, and ma- 
chinery large radio telescope. 

Mr. Feld presents list eighteen severe storms which have previously 
been reported the appropriate literature and for which has reduced the 
recorded wind velocities allow for instrumental errors. This im- 
portant since the long-time records wind velocity reported the USWB are 
not all the use the same instruments, the same methods 
reporting. has indicated later paragraph his discussion. This 
list may prove valuable those who feel that particular structure deserves 
closer study probable basic wind velocities than can obtained from Figs. 
also cites the failure radio towers which subsequent stress 
analysis failed check recorded wind velocities. mentions the possibility 
that the higher kinetic energy driving rain may have been responsible. 
This question has been under study recent years but the writer feels that 
the results are not yet ready for inclusion recommendations the Commit- 
tee Wind Forces. The writer knows acase where subsequent stress anal- 
ysis has shown that the failure include the “moving load” effect the design, 
Fig. was sufficient account for the failure guyed tower. 

The writer indeed grateful Mr. Davenport for his extended discussion, 
especially regard the effect surface roughness onthe variation wind 
velocity with height. His remarks contrasting surface roughness with the sta- 
bility the atmosphere controlling factor the variation wind veloc- 
ity with height are pertinent the later discussion the Brookhaven records 
this closure. This not because the degree surface roughness Long 
Island stands alone determining the variation wind velocity with height, 
because these records illustrate the stage where the transition from conditions 
roughness over water those over land major influence. 

not practicable toseparate the influences surface roughness and at- 
mospheric stability the usual cases, since both influences are variable. 
For example, even uniformly wooded area the coefficient roughness 
will decrease with increasing wind velocity, due the streamlining the leaves 
summer and the twigs winter they bend the wind. Nevertheless, 
may useful todiscuss, briefly, the question atmospheric stability and sur- 
face roughness influences the variation wind velocity with height. 
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The widely accepted exponential relation given 


x 


where u(Z) the average velocity height and denotes the average ve- 
locity height Sutton (13) gives the equation the form 


n/2-n 
u(Z) = uy E ee (2) 


where parameter which varies with the conditions stability. Sutton 
states that, “No reliable values are yet available for large lapse rates and large 


Brookhaven Storms (Coastal Area) 


fastest peaks 
Carol 1954 Extremes 


Edna 1954 Extremes 


Height Above Ground (Feet) 


FIG, 9.—VERTICAL WIND PROFILE VS, STABILITY 


inversions.” Hewson (14) states that “has the value 1/4 under average 
conditions lapse rate, witha range from 1/2 for large inversion 1/5 for 
large lapse.” 

The curves showing the variation wind velocity with height for each 
these three values are shown The data for the three Brookhaven 
storms (15) are also shown and will seen that they lie the right the 
line which indicates average lapse rate 0.25), that is, they lie the area 


800 
600 
400 
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moderate stability, where relatively rapidincrease velocity with height 
may expected because the small amount vertical interchange present 
under conditions moderate stability. 

The vertical temperature gradients are available (15) for the two hurricanes 
Carol and Edna (1954), but not for the violent extra-tropical storm Novem- 
ber 25, 1950. The temperature gradients which are reported for Carol and 
Edna are not inconsistent with the rapid variation wind velocity with height 
shown Fig. 

will shown subsequently that the wind passes from the water area 
the landarea must undergo atransition from vertical profile that com- 
patible with the surface drag over the water one that compatible with the 
surface drag over land. This quite different however, from where 
condition surface roughness. Consequently, the explanations the 
crease velocity with height must sought some combination surface 
roughness and atmospheric stability. 

may that some future Committee Wind Forces will see fit include 
additional tables based several different conditions surface roughness, 
but this time not recommended. 

The writer not certain that agrees with Mr. Davenport’s statement, 
given his item (2), that there “no marked tendency for violent thermal 
interchange” large scale, mature, extra-tropical storms. the zone 
transition, following the passage the cold fronts, that the most violent and 
prolonged gusts occur extra-tropical storms and these are caused thermal 
interchange between the higher and lower strata the storms (16). the 
statement restricted tothe last stages maturity where decay impending; 
that is, where occluded situation has set impending, then true 
that thermal interchange may small. 

There considerable duplication the discussion concerning the variation 
wind velocity with height hurricanes, and therefore, this closure will deal 
with the subject general rather than the discussion separately. Some excep- 
tions this procedure will made where seems desirable for clarity. 

has been assumed the writer that the “first gust” thunderstorms 
and the “plow winds” hurricanes, mentioned Mr. Bunting, have both been 
included the statistical analysis upon which Figs. and are based. 

The writer does not understand the statement Mr. Bunting that “the com- 
puted maximum gradient winds were obtained from averaging over 
minute intervals that would improper compare observed minute 
maximum winds with computed minute maximum gradient winds.” The gra- 
dient wind equation based the assumption steady-state conditions and 
contains time factor (27) (28): 


which the angular velocity the earth’s rotation, denotes the lati- 
tude, the radius curvature the wind path, density, and 
the gradient wind. 

The Ekman spiral fitted the storm data taken Ann Arbor (17) shows 
that the time required establish steady state conditions appropriate for the 
spiral and gradient wind theories lay somewhere between min and min, 
while the Ekman spiral Curve Fig. 2A, fitted the peaks the storm 
November 25, 1950 Brookhaven, shows that the required time lay somewhere 
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between min and the period the peaks. course, these two storms were 
extra-tropical character while those quoted Mr. Bunting were definitely 
tropical storms. 

states that the work quoted Mr. Bunting the “gradient winds” 
were really estimated cyclostrophic winds obvious, from the 
discussions, that these “gradient winds” cannot used the conventional 
meaning the term, was done Fig. and Table and, consequently, 
curves and Fig. must discarded together with the design 
velocities above mph Table 

Mr. Feld has also included abstracts from material reported Smith and 
Singer (15). based observations made the 410 tower Brook- 
haven, I., during the severe storm November 25, 1950, and the two hur- 
ricanes 1954, Carol and Edna. The following comments will draw onthe re- 
port Smith and Singer, rather than the abstracts quoted Mr. Feld: 

The records from the storm November 25, 1950, were discussed under 
the heading “Extra-Tropical Cyclones.” The velocities were reported the 
fastest 6-min means and the fastest peaks each the tower stations. The 
data were plotted Curves and Fig. and Ekman spiral was fitted 
the peak values and extrapolated the gradient level. Since these are ex- 
treme values they may used directly the Committee Wind Forces. 
However, the records from the two hurricanes were reported “Means and 
Standard Deviations 1-Min Mean Speeds, 2-Sec Lulls, 2-Sec Peaks,” and for 
the purposes the Committee Wind Forces some approximationto extreme 
values must made statistical extrapolation. 

Since only the averages deviations are given, will assumed 
that the distribution curves for the data are not skewed and that the symmet- 
rical normal law may used acceptable approximation choosing ex- 
treme values. This assumption supported many observations made during 
storms Ann Arbor, example which shown Fig. from the storm 
January 19, 1933. Fig. 10(a) shows the variation skewness with height 
for the distribution 1/2 sec average velocities within min intervals, for 
total hr. The positive skewness the bottom stations the tower 
indicates preponderance low velocities with less frequent high velocities, 
while the negative skewness the upper stations indicates preponderance 
high velocities withless frequent low velocities. height about 125 ft, 
the high and low velocities have equal chance occurrence, that is, the 
distribution curve has zero skewness. assumed that acceptable ap- 
proximation extreme values may obtained adding 1.5 times the stand- 
ard deviation Fig. 10(b) shows that least 93% 
all cases will included regardless which skewness used. The data for 
curves and Fig. 2A, for Carol, 1954, were derived this way, were 
the data Table 

The storm winds observed Brookhaven are typical the transition zone 
which present where the wind crosses coast line and blows over land 
surface after having blown over extensive water surface such the ocean, 
Long Island Soundor the Great Lakes. The coast lines this case are irreg- 
ular and impossible define exactly. The winds blew from the two quadrants 
the compass from the east through the north the west, and the coast line 
varied between approximately miles and miles. 

Storm winds approaching coast line have velocity profile which com- 
patible with the pressure gradient, the curvature the wind path, the rough- 
ness the surface, and other things which affect the air movements. 
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FIG, SKEWNESS WIND DATA VS, HEIGHT ABOVE GROUND. 
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FIG, 11.—RAPID INCREASE WIND VELOCITY WITH 
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the wind passes the coast line there aretwo main influences which are re- 
sponsible for slowing down its velocity. First, increased fric- 
tional drag which generates increase turbulence and this, turn, diffuses 
the drag upwards. Second, the supply water vapor longer available 
from the ocean and consequently can longer provide source energy 
wherein water vapor condensed into cloud droplets with the accompanying 
release latent heat. The greater the distance which the wind moves inland 
from the coast line, the greater will the height which the increased tur- 
bulence will effective modifying the wind profile. 

Fig. diagrammatic representation the foregoing idea regarding 
the role frictional drag. Brunt (18) gives approximate equation based 


TABLE 2.—EXTREME VALUES WIND VELOCITY (MPH) 
BASED BROOKHAVEN DATA 


Height, 


Minimum mean speeds. Second peaks. 


previous work Taylor, for the “height which turbulence effective 


time interval the eddy diffusion heat. may used also for the 
eddy transfer velocity. 


which the height (in centimeters) which turbulence significantly 
effective time interval (in seconds), the eddy diffusivity 10° 
the distance inland from the coast line (in centimeters) and 
the average wind velocity (cm height time then 


Because the approximations used the derivation Eq. should 
considered giving only the order magnitude The degree approx- 
imation also affected the loss water vapor source energy when 
the away from the coast causing the upper winds todecrease more 
than the amount caused surface drag. The moderate thermal stability, pre- 
viously discussed, alsoignored Eq. The value 10°, recommended 
Brunt, has been successfully used Hewson (19) for investigation 
which was approximately mph for period about hr. 


2 
Carol Edna 
47.9 
150 69.2 
355 104.6 
| 
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FIG, 12.—HEIGHT WHICH TURBULENCE EFFECTIVE TRANSITION ZONE, 
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Fig. shows curves fitted tosome for different val- 
ues and From consideration these curves and the foregoing dis- 
cussion, concluded that coastalareas, such Brookhaven, Long Island, 
the rapid increase wind velocities with height need not continued beyond 
600 ft, above which constant velocity may used. 

The Brookhaven data for Edna, given Table are plotted 
Fig. together with the curves for the 6-min winds and the peaks the 
storm November 25, 1950. Three 0.3-power curves are shown for compar- 
ison. They seem indicate that the exponent 0.3 gives conservative esti- 
mate the variationof windvelocity withheight during these three Brookhaven 
storms. 

Fig. assigns basic wind velocity mph mph the Long Island 
area (fastest mile level, 50-yr period return). This far 
greater than the extreme 1-min winds shown Table for Carol and Edna 
the level and even exceeds the extreme 2-sec peaks for this level. This 
indicates that, this coastal area, even stronger winds must expected 


the structural designer than those observed during Carol and Edna. Two ques- 
tions arise: 


(1) Should the very rapid increase wind velocity with height 0.3), 
shown these two hurricanes their decaying stages and the one intense 
extra-tropical storm, used for future guidance this coastalarea when the 
specified basic wind (75 mph mph) twice great the basic veloc- 
ities accompanying the observation (Carol 43.8 

(2) Should used those coastal areas which may visited hur- 
ricanes during their mature stages? 


Mr. Gentry quotes the fastest mile wind recorded the USWB Miami 
for each five hurricanes. compares the velocities the Airport Station 
(50 ft) and the City Office (229 ft). The ratios are, respectively, 
1.01, 1.14, 1.26, 1.30, and 1.53, and the corresponding values (the expon- 
ent) are 0.0, 0.082, 0.083, 0.085 and 0.28. Only the largest the exponents 
compares all closely with those obtained Brookhaven, shown the 
appropriate curve for the corresponding basic wind Fig. 13. Mr. 
Gentry points out that “these data unfortunately cannot used satisfactorily 
for determining the variation the wind with height because the different 
exposure the two stations and the distance between them (approximately 
miles).” might also added that each pair the readings are not coinci- 
dental, that the two stations are different parts the transition zone for 
winds blowing from the ocean, and that the anemometer the City Office 
mounted atall office building which affects the wind differently for different 
wind directions. All this emphasizes the difficulty obtaining reliable 
quantative information regarding the variation wind velocity with height 
any situation except where anemometers are mountedon exposed tower. 
these five hurricanes, the increase maximum with 179 increase 
height varied from 53%, and only one storm five gave results com- 
parable the storms Brookhaven. Although both Miamiand 
transition zones coastal areas one likely experience mature hurri- 
canes, while the other, the two recorded hurricanes were decaying 
stage. 

Riehl states that “no one knows how high wind speeds can rise trop- 
ical storms.” Also, that wind speed 200 mph may attained “over the 
sea and few hundred feet above the ground” very small area. also 
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quotes Hughes (21) the effect that average knots (103 mph) 
was observed near the center 28storms during 84reconnaissance flights 
the United States Navy near the 1,000 level. This average value should 
increased two factors least 1.3 each, one obtain the extreme storm 
and the other obtain the fastest mile wind within that storm. This would 
give “fastest mile wind” least 175 mph the 1,000 level. 

Fig. does not specify basic velocities excess 120 mph but even this 
value would require the fastest mile wind 600 295 mph 
the 0.3 exponent were used. Also, for the benefit engineers who may feel 
the need still higher basic winds unusual situations, basic winds high 
130 mph are included Table 3(b). The accompanying wind speed 600 
would need 320 mph 0.3 were used exponent. 

The velocities 295 and the 600-ft level are believed ex- 
cessive, but the writer knows well-controlled observations mature 


TABLE 3.—FASTEST MILE WIND (MPH) FOR VARIOUS 
HEIGHT ZONES ABOVE GROUND (SEE FIGS. AND 


Height Zone, Basic Wind-Velocit 


(a) For Inland Areas 
100 115 120 
100 110 125 
105 115 130 145 
(b) For Coastal Areas 
100 105 115 
125 130 135 145 150 
140 150 160 165 170 175 180 185 190 


400 
700 1,000 
1,000 1,500 


600 1,500 


hurricanes, such those which were made onthe tower Brookhaven for the 
two decaying hurricanes. Consequently, prudent compromise between reck- 
lessness and extravagant conservatism requires careful consideration any 
relevant material from which reasonable inferences may drawn. be- 
lieved that itis the duty the Committee Wind Forces make recommen- 
dations for the guidance structuralengineers evenin this areato which field 
research has not yet penetrated. 


The answers the two questions posed ina previous paragraph will 
adopted follows: 


(1) coastal areas where the basic wind mph, the variation wind 


velocity with height will adpoted for heights the 


600-ft level, above which constant velocity will used. 


. 
150 
400 
130 
150 150 
400 180 
400 600 195 
200 


ST7 


DISCUSSION 


211 


(2) coastal areas where the basic wind mph, variable exponents 
will adopted such that never greater than 200 mph 600 Table 


3(b): 
V30 


varies from 0.3 when mph 0.143 when 130 mph. Curves 
are shown Fig. for some the intermediate values 


the present state information seems undesirable for the Committee 
Wind Forces totry between the wind velocities which, coastal 
areas, are appropriate for structures the coast-line and those located 
farther inland. recommended that coastal areas considered extend 
least miles inland (22) from well-defined coast line. 

For all situations which, the judgement the engineer, deserve the clas- 
sification “coastal area,” the wind velocities shown Table 3(b) should 
used together with the “Basic velocities” shown Figs. Special anal- 
ysis should made for promontories and areas where the coast line 
precipitous. 

For all inland areas Table 3(a) should used together with Figs. 
These values are minimum increased where higher basic veloc- 
ities have been recorded locally, are likely occur, due mountains, 
gorges, for other unusual reasons. 

There considerable decrease the density the air locations high 
above sea level and consequently there areduction the dynamic pressure 
(velocity pressure) compared the same wind speed and temperature 
sea level. If, however, the fastest mile wind, for which the contours Figs. 
and are drawn, should occur temperature greatly below that the 
standard atmosphere (59° F), the loss density due elevation may com- 
pletely compensated the increase density due lower temperature. For 
example, Brekke (23) gives table which shows that, 4,000 elevation and 
standard temperature, there decrease 14% the density compared 
sea-level elevation, but the temperature were reduced to-10° there would 
overall reduction density the atmosphere. Consequently, rec- 
ommended that the engineer make reduction wind pressure based re- 
duced density high elevations unless has investigated the probability that 
the fastest mile windused the design may coincident with significant 
reduction temperature. 
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SUMMARY RECOMMENDATIONS 


The scope these recommendations has been limited those general con- 
siderations which are common most problems the selection appropri- 
ate wind velocities for the design bridges, buildings, towers, and similar 
structures. has been assumed that highly specialized problems will receive 
the attention engineers with competence make direct the necessary 
modifications additions the recommendations. 

Several important changes have been made the original recommendations 
this paper. from additional information furnish- 
the writers the discussions. Consequently, complete summary 
revised recommendations given here: 


All recommendations herein are based conditions open, level coun- 
try standard reference, and the influence shielding, deflections 
and channeling due unusual topography orto large obstructions shall eval- 
uated the engineer each individual 

The effective ground (that is, the height open, level country its 
quivalent) which minor obstructions are unimportant, and which approxi- 
mates the height anemometers now use the USWB the first-order 
stations airports, shall taken ft. 

Basic windvelocity shall definedas the fastest-mile-of wind, reduced 
height ft, shown Figs. and for the various geographical 
reas the United States. Fig. shows the basic winds which have statisti- 
cally derived probable period recurrence and Fig. probable pe- 
riod recurrence 100 yr. 

wind velocity with a50-yr period recurrence (Fig. may 
used for all permanent structures except those which, the judgment the 
engineer, have high degree wind sensitivity and unusually high degree 
hazard life and property case failure. the latter case, 100-yr 
recurrence should used 6). the case temporary struc- 
tures, such those used during construction only, basic winds having veloc- 
ities equal three-quarters the velocities shown Fig. may used 
the discretion the engineer. 

special locations, such promontories, mountains, gorges, 
for other unusual reasons, where records experience indicate that the 
basic velocities given Table are inadequate, higher basic veloc- 
ities shall used the discretion the engineer. 

Table (a) shall used inland areas obtain the variation wind 
velocity with height after Figs. has been used obtain the minimum 
basic wind velocity for that area. For each basic wind velocity stepped ap- 
proximation the one-seventh-power law has been used Table (a) toa 
height 1,000 above which constant velocity has been used. 

Table used coastal areas obtain the variation wind 
velocity with height after Figs. has been used obtain the minimum 
basic wind velocity for that area. For each basic wind velocity stepped ap- 
proximation power law has been used height 600 ft, follows: 
(a) where the basic windvelocity equal less than mph, exponent 
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0.3; (b) where the basic wind velocity greater than mph, exponent 
varying from 0.3 0.143 130 mph. Above 600 feet constant velocity 
used each case. Coastal areas shall considered extend least 
miles inland from well-defined coast line along the oceans, great lakes, 
other large bodies water. 

Appropriate gust factors shall used for structures which are small 
enough responsive gusts involving less than one mile passing wind. 
The gust-factor (gust velocity velocity the fastest mile wind) should 
bear some relation the minimum size gust necessary envelope the 
structure and its accompanying pattern flow (24). gust-factor 1.3 will 
allow for gusts approximately one-second duration (25) which, mph 
basic wind, would have downwind about 130 ft. This will adequate 
for sign boards andsmall gust-factor 1.1 willallow for gusts 
approximately 10-sec duration which, mph basic wind, would have 
length downwind about 1,300 ft, and which adequate for structures having 
horizontal dimension transverse the wind about 125 ft. 

Guyed towers shall for the moving load effects aerial jets 
(26) and descending gust fronts, indicated Fig. 

10. those geographical areas storms may followed strong 


winds, least 1/2 in. thickness ice shall assumed all structural 
members, 
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CONTINUOUS GIRDER BRIDGE WITH VARIABLE MOMENT INERTIA* 


forces continuous beams with variable sections has been developed Mr. 
Sami paper. simplify this method, Mr. Sami took advantage the 
symmetry the spans. recent years, many simplified analytical and graph- 
ical methods have been developed. 

solving the same problem with various computation methods, the advan- 
tages the methods With this purpose the illustrative prob- 
lem shown Mr. Sami solved bythe graphical method moment distribu- 
tion (Fig. 18). The distance the central point determined from 


The central-point location was explained elsewhere.® Line the closing 
line the fixed end bending moments which may determined any method. 


Point 


FIG, 18.—BENDING MOMENT DIAGRAM 


January, 1960, Sabri Sami. 
Assoc. Bridge Engr., Calif. State Highways, Sacramento, Calif. 
Discussion “Arches and Bents,” Transactions, ASCE, 1950, 630. 
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With graphical construction athree line polygon and cross lines, the bending 
moment diagram drawn. 

All maybe The addi- 
tional graphical constructions for the varied properties beams may com- 
pleted fraction this time. 

interesting note that the redundant forces computed Fig. 18, differ 
from those computed the author amount which less than 3%. The 
graphical construction the bending moments may further simplified 
drawing the bending moment closing line directly through the central point 
(vertical GH). 

conclusion, may stated that, selecting the computation method 
solving continuous beams, important consideration should given the de- 
signer’s personal conception expressing the redundant forces either nu- 
merical factors, graphical ordinates. 

performed tremendous algebraic work simplify the expressions 
for the redundant forces. practical application the method only the final 
equations diagrams are needed. Yet, the designer always has prepared 
trace the effect various features continuous beams. This fact may 
found quite involved. gratifying that the author performed such tremendous 


work illustrating quite unique approach solving continuous beam with 
redundant forces. 


STABILITY CONSIDERATIONS THE DESIGN 
STEEL PLATE 


EREMIN,2 this paper Mr. Massonnet has presented some 
interesting tests steel plate girder stability. The tests have proveu that the 
stresses inthe stiffeners, web plate, and girder flanges are interrelated. From 
the tests may concluded that tight fitting the ends intermediate ver- 
tical stiffeners does not contribute uniform distribution stresses the 
web plate. Likewise, the stresses the horizontal stiffeners (when they are 
placed short distance from the flanges) may affected the stresses 
the flanges. 

present, there tendency for continuous steel plate girders de- 
signed with haunched web plates the intermediate support. would inter- 
esting extend the tests the steel plate girders with variable width web 
plate. Furthermore, the stiffners are often combined with intermediate dia- 


phragms wind cross frames. Obviously, the stresses the cross frames 
the distribution stresses inthe stiffeners. This may especially 
true relation the horizontal stiffeners. 


January, 1960, Charles Massonnet. 
Assoc. Bridge Engr., Calif. State Highways, Sacramento, Calif. 
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PROPERTIES STEEL AND CONCRETE 


AND THE BEHAVIOR 


x 


Discussion Lester Gabriel 


LESTER ASCE.—The author quite rightly points the 
fact that the material testing procedures commonuse and the interpretations 
placed thereon are attempt make explicit some the more significant 
variables design model thought apply. 

design model (method), the consequences the design are, neces- 
function not only the explicit factors material composition, stress, 
strain, modulus elasticity, yield point, etc., but also many (perhaps hun- 


dreds even thousands) neglected factors. good design method will sat- 
isfy the following: 


(1) The variables the properties the material that have significant 
effect will well understood and properly evaluated that the explicit fac- 
tors will properly influence the design procedures. 

(2) All the major factors influence will have been made explicit and 
properly evaluated that those many factors which are considered negligible 
and omitted are fact negligible and act chance variable. negligible 
factors influence act something other than chance variable then the to- 
tal effect the negligible factors may quite significant. 

(3) The design method itself will subject test and the consequences 
thereof will reflect the value the method. the results the tests fail 
satisfy the constraints the design method, then the entire method itself 
must scrutinized for erroneous assumptions; that the explicit factors, 
the negligible factors that are assumed act chance variable, and the 
test instruments themselves (the maintained hypotheses). 

those initiated the area reinforced concrete design the elastic 
analysis method, immediately apparent just how important the modulus 
elasticity concrete fact and the ratio the modulii steel 
concrete, underly the very substance the method the transformed sec- 
tion for elastic design. The ACI Building calls for which 
turn, gives value 1,000 This the well known secant modulus 

and the design methods present and past accept this estimate one very 
significant factor the many explicit factors the design method. 


February, 1960, George Winter. 
Asst. Prof. Civ. Engrg., Sacramento State College, Sacramento, Calf. 


“Building Code Requirements for Reinforced Concrete,” American Concrete In- 
stitute, ACI, 318-56. 
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The fact that known that the secant modulus 1,000 
low for moderate strength concrete and high for strength concrete,8 
coupled with the fact that any kind variability ignored the clas- 
sical design methods, puts this variability the group neglected fac- 
tors. The fact that the variability does not present itself chance 
variable but rather afunction itself, leads abias the design method. 
How then does this one type variation (which not subject chance alone) 
effect the strength the structure? How does effect the apparent factor 
safety How variations other factors the design model, both explicit 
and neglected, effect the strength the structure the apparent factor 
safety 
The writer wishes lend support the author’s statements 


“The attempt reduce the amount judgment and increase the fac- 
tual and quantitative basis which safety and load factorg are deter- 
mined. While unlikely that will ever possible deal with safe- 
entirely statistical-probabilistic basis, ------- (information) 
can incorporated the safety factor explicitly and numerically.” 


This writer feels very strongly that now—when new methods design and 
analysis are bursting forth, when older methods are being re-examined, and 
when the philosophy the action structures under load (s) and all the de- 
sign methods related thereto are being explored the time attempt 
include probabilistic estimates variability parameter the methods 
themselves. This would lead (1) more significant testing and test proce- 
dures, (2) realization the part the designer that variability does exist 
and must considered the design, and (3) more sophisticated concept 
safety factor load factor terms probability failure and the related 
consequences. 

Once design method developed and accepted “modus operendii” 
the design profession, becomes exceedingly difficult effect changes. 
with anxiety that this writer has followed the development plastic design 
that the margin safety design “is not less than that provided against fail- 
ure simply supported beam according usual past practice.”9 

This writer believes that the structural engineering profession has arrived 
that degree sophistication whereby notions factor safety margin 
safety that describe vague condition expected strength without describ- 
ing some estimate likelihood failure are grievously wrong. Variations 
all kinds must considered and part the thinking the developers 
design methods. develop new methods they are safe old methods 
step backwards because the time lag generates the eventual defini- 


“Reinforced Concrete Fundamentals,” Phil Ferguson, John Wiley and Sons, 
Inc., New York, 1958, 10. 


“Plastic Design Steel Frames,” Beedle, John Wiley and Sons, Inc., New 
York, 1958, 


DIGEST THE GUIDE DESIGN CRITERIA 
FOR METAL COMPRESSION 


Discussion Richard Pian 


RICHARD PIAN,3 ASCE.—One the objections the initial-yield- 
criterion approach beam-column design, outlined the author, fol- 
lows: “The effect end restraint cannot rationally introduced into the se- 
cant formula analysis since the procedure basedon the analysis unre- 
strained member.” 

The writer has developed generalapproach for the analysis metal col- 
rational basis similar the secant formula. non-dimensional 
formula for column design with different degrees imperfections, different 
end conditions, and with materials having different yielding stresses can 
written the following form. 

(1) For columns with initial deflection 


(2) For columns with initial eccentricity 


which the axial stress, denotes the yielding stress, 


set curves may plotted from Eqs. for various values 
and these curves can used for any Oy. number other curves are in- 
cluded the four-quadrant chart 8). The third quadrant curves, which 
are first used for computation, give the relation between L/r and 
Here the effect end conditions compression members can taken into 


Prof. Engrg., Arizona State University, Tempe, Ariz. 


“Column Design Simplified,” Richard Pian, Proceedings, ASCE, May, 
1957. 
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The term definedas the effective length the column. Next, the straight 
lines the fourth and second quadrants are drawn for different values 


The curves the first those plotted from values obtained eval- 
uating Eq. Eq. 36. 
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